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ABSTRACT 


An analysis considering nonuniform radiosity and assuming diffusely emitting and re- 
flecting surfaces is presented. The shields are planar and circular, while the sources 
may be surfaces of revolution. The effects of source geometry, number of shields, 
spacing, and surface properties are discussed. Shield conductivity and selective coatings 
are investigated. Specularly and diffusely reflecting surfaces are discussed under the 
simplifying assumption of uniform radiosity. Also, the effects of thermal radiation on 
the temperatures of a conducting structural member are analyzed. Results are presented 
showing the effects of significant parameters on the temperature profiles and heat- 
transfer rates. The computer program for the shield analysis is listed. 



THERMAL ANALYSIS OF SHADOW SHIELDS AND STRUCTURAL 
MEMBERS IN A VACUUM 
by Robert J. Boyle and Richard H. Knoll 

Lewis Research Center 
Cleveland, Ohio 

SUMMARY 

An analysis is presented which predicts heat-transfer rates as a function of the num- 
ber of shadow shields used. The system analyzed consists of two sources and the shields. 
One of the sources may be solar radiation. It is shown that shields form an effective 
means of reducing radiant heat transfer. A possible application of these shields is the 
reduction in the heat transfer to a cryogenic propellant. 

The shields are flat disks while the sources may be surfaces of revolution. The 
sources and shields have a common axis of revolution. The prime analysis used is for 
diffusely emitting, diffusely reflecting surfaces assuming nonuniform radiosity. The anal- 
ysis considers thermal radiation and radial shield conduction. The effects of the simpli- 
fying assumption of uniform radiosity are considered. It is concluded that for many cases 
it may be necessary to use the more complicated assumption of nonuniform radiosity in 
order to have accurate heat-transfer predictions. Specularly reflecting surfaces are ex- 
amined under the simplifying assumption of uniform radiosity. As the spacing between 
surfaces increases, the heat transfer for specular surfaces may significantly exceed that 
for diffusely reflecting surfaces. 

In addition to the analysis for the shadow shields, an analysis is presented to predict 
the temperature distribution along a cylindrical conducting strut. This analysis for the 
strut considers both internal and external radiation. The struts are needed in order to 
attach the shields to the rest of the vehicle. The thermal analysis for the strut concludes 
that increasing the external emissitivity is effective in lowering the conducted heat trans- 
fer out of the strut. The interaction of the struts and shields was not considered in the 
analyses. 

For both the shield and strut analyses, predicted temperatures are compared with 
experimental data . Heat-transfer rates for the shields are also compared with experi- 
mental data. 



INTRODUCTION 


The design of space vehicles for long-term missions often requires thermal protec- 
tion systems. This is especially true when cryogenic propellants are involved. Shadow 
shields are a means of reducing the radiant energy absorbed by a propellant tank. This 
results in reduced boiloff . Shadow shields can also be useful in protecting a relatively 
warm body, such as a payload, from receiving excessive amounts of solar energy. The 
primary area of application for shadow shields is in a vacuum. Otherwise, convection 
may be the dominant mode of heat transfer. The basic purpose of this report is to present 
an analysis for multishield configurations in a space environment. 

Previous studies considered shadow shields for thermal protection of cryogenic pro- 
pellants on long-term missions, and the protection of payloads passing close to the sun. 
These studies considered flat disk, spherical, and conical shields (refs. 1 to 5). The 
analyses in these references were generally for shadow shields of uniform temperature. 
Reference 3 analyzed flat disks exposed to solar radiation for both zero and infinite radial 
conductance . Reference 6 determined the local and overall heat-transfer rates for two 
flat circular sources. In references 3 and 6, the surfaces were treated as having nonuni- 
form radiosity. 

In this report, the effects of finite shield conductance and variable surface properties 
are considered. The analysis presented is applicable for any number of shields, and data 
are given for up to three shields. Assumptions made concerning the radiosity and reflec- 
tivity of the surfaces can affect the predicted heat-transfer rates. The prime analysis 
assumes diffusely emitting and diffusely reflecting surfaces of nonuniform radiosity. The 
consequences of making the simplifying assumption of uniform radiosity are examined. 
Specularly reflecting surfaces are also analyzed under this simplifying assumption. 

The shields stand off from the surface which is being protected. It is necessary to 
support the shields, and an analysis for the appropriate structural members is also pre- 
sented. However, the interaction between the shields and supports is not analyzed. 
Therefore, the thermal analysis for the shield support is applicable to cylindrical mem- 
bers in space. The analysis for the support gives the axial temperature distribution along 
a tube or rod. This analysis considers both conductive and radiant heat transfer. 

Comparisons for both the shield and strut analyses are made with experimental data 
reported in reference 7. A computer code for the shield analysis is listed and discussed. 
The code is for diffusely emitting and reflecting surfaces of nonuniform radiosity. 
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METHOD OF ANALYSIS 
Shield System 

The steady-state temperature distributions in a shield system and consequently the 
heat -transfer rates are determined by talcing heat balances on elemental areas of each 
shield. These heat balances account for the energy emitted, absorbed and reflected as 
well as conducted to or from the element under consideration. All of the elements on the 
shield as well as the surroundings and the elements of the adjacent surfaces enter into the 
radiant heat terms. Adjacent elements of the same shield also enter the conduction terms 
of the balance. Appendix A gives the symbol list used in the analyses. Appendix B pre- 
sents the equations for a system consisting of shields and two sources. One of the 
sources may be solar radiation. Figure 1 gives a schematic of two possible shadow- 
shield configurations. From an analytic standpoint, sources are similar to shields ex- 
cept that they have temperature distributions which are known a priori. The sources can 
be at any temperature; therefore, they may act as sinks. 

In order for the analysis to be applicable to a system of shadow shields, the system 
must satisfy certain requirements. These conditions are: 

(1) There can be no more than two sources and these may be surfaces of revolution 
but they cannot see themselves. 

(2) The shields have to be circular and planar. 

(3) The shields and the sources have a common axis of revolution. 

(4) The size of the shields and sources are such that a surface of either a source or 
shield sees only one surface. When one of the two sources is solar radiation, the outer- 
most surface, and only this surface, has solar radiation incident upon it. 

The derivations of the governing equations in appendix B are made under the follow- 
ing assumptions: 

(1) The environment between each pair of surfaces has a single temperature and con- 
stant reflectivity. The environment is that space which completes the enclosure between 
two surfaces. 

(2) There is no thermal gradient in a shield in either the circumferential or axial 
directions. 

(3) Heat is transmitted to or from a shield only by radiation. Under this assumption, 
each shield is isolated from its support structure. 

(4) While each shield has a finite thickness, it is assumed that there is no radiation to 
or from the edge of the shield. 

(5) The equations are derived for steady-state conditions. 

(6) All surfaces are opaque so that a + p = 1 at every point on the surfaces. 
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(7) E missivity and absorptivity can each be separate functions of radial position and 
temperature, but are constant in the circumferential direction. These properties are 
taken to be only a function of the temperature of their surface. Emissivity and absorp- 
tivity may be functions of wavelength, but this is not accounted for in the analysis. 

(8) Each shield is homogeneous so that the thermal conductivity of the shield is a 
function only of temperature. 

(9) Each of the surfaces is diffusely emitting and diffusely reflecting. 

A comparison is made in a later section of the report between the heat -transfer rates 
for specularly reflecting surfaces and diffusely reflecting surfaces. The equations used 
for determining the heat -transfer rate between specularly reflecting surfaces of uniform 
radiosity are presented in appendix C. In the analysis for specular surfaces, it is as- 
sumed that the reflectivity is independent of the angle of incidence. Also, polarization 
effects were neglected. 

The equations presented in appendix B for predicting the performance of a shadow- 
shield system result from each of the surfaces being divided into a series of concentric 
annuli and heat balances being made on each annulus. Heat enters an element due to solar 
radiation absorbed on the surface or from thermal energy radiated from an adjacent sur- 
face. Energy leaves the annulus due to the emissive power of the element. Additional 
heat enters or leaves the element due to conduction from the two adjacent elements on the 
same shield. As the number of annuli increases, each of the annuli approach a differen- 
tial area in size. When the annuli represent differential areas, the results obtained are 
for surfaces of nonuniform radiosity. When only a single annulus is used for each sur- 
face, the results obtained are those for uniform radiosity surfaces. Unless otherwise 
noted, the results in this report are for diffusely reflecting surfaces of nonuniform radi- 
osity. In connection with the discussion on specularly reflecting surfaces, a discussion 
is included on the errors which result from assuming surfaces to have uniform radiosity. 


Thermal Distribution Along Structural Member 

The analysis in appendix D yields the steady-state, one -dimensional temperature 
distribution in the axial direction along a structural member which may be either a rod or 
a tube. Figure 2 gives a representation of a strut. In order for the analysis to predict 
the temperature distribution for a rod accurately, there should be low thermal resistance 
in the radial and circumferential directions. The environment is taken to have a uniform 
temperature and reflective properties. As in the derivation of shadow-shield equations, 
it is assumed that the surfaces are diffusely emitting and diffusely reflecting. It is nec- 
essary to divide the strut into several segments due to the axial temperature variation; 
and, therefore, the results obtained for the internal radiant heat transfer approach the 
case of nonuniform radiosity. 
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The governing equations are determined by taking a heat balance on the elements of 
the strut. Heat enters an element of the strut by conduction from an adjacent warmer 
element and may be conducted away to an adjacent colder element. Additional heat enters 
or leaves the element by thermal radiation from the inside surface of the tube. If the 
ends of the tube are closed, the end surfaces are taken to be at the same temperature as 
the ends of the tube. These two end temperatures are the known boundary temperatures. 
If the ends of the tube are open, the end disks have the properties of the surroundings, 
which are assumed black, but not necessarily at zero temperature. 

The equations presented in appendix B for the shadow shields, and in appendix D for 
the strut are relatively complex. The solution of these equations is discussed in appen- 
dix E, and appendix F presents a computer code for solving the shadow -shield equations 
of appendix B. 


COMPARISON OF RESULTS 

This section of the report discusses the comparison of experimental and analytic re- 
sults. In addition, mention is made of previous analyses with which the work presented 
herein was compared. 


Shadow Shields 

The analysis presented in appendix B yields results consistent with the assumptions 
of uniform radiosity surfaces when each surface is taken to have a single annulus. Under 
this condition, the results are comparable with those presented in reference 1. 

When the analysis is performed using a large number of annuli for each surface, the 
results obtained approach those for nonuniform radiosity surfaces. In reference 6, ana- 
lytic results are presented for the heat -transfer rates between two surfaces wherein it is 
assumed that the surfaces have nonuniform radiosity. Solving the equations of appen- 
dix B for the case of nonuniform radiosity gives results which agree with those of refer- 
ence 6. 

Experimental data for shadow-shield systems are given in reference 7 and some of 
these data are reproduced in the following figures. Figure 3 gives a comparison of ana- 
lytic and experimental heat -transfer rates as a function of spacing for a one- and three- 
shield system. The shields were made of copper so that their radial temperature distri- 
bution was nearly uniform. The shields were evenly spaced between two plane sources 
with the warmer source being maintained at 444 K (800° R) and the colder one at 77. 8 K 
(140° R). This figure contains data for both a one- and three -shield system at two dif- 
ferent shield emissivities. In one case, the shields were painted to give a nominal 
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emissivity of 0.94, the same as the sources. While in the other case, the shields were 
glass blasted to give a nominal emissivity of 0.27. (The surface properties of both the 
painted surfaces and glass-blasted shields were determined to be a function of tempera- 
ture so that there was a variation in emissivity from one surface to the next.) As is ex- 
pected, the heat transfer decreases with the addition of more shields and increases with 
increasing emissivity. 

As the emissivity approaches unity, the difference in the heat-transfer rates between 
assuming uniform and nonuniform radiosity goes to zero. Therefore, the analytic predic- 
tions for the high-emissivity cases are nearly the same using either assumption. Only a 
single curve is shown for each of these cases. The single glass-blasted shield (e « 0.27) 
was located between two high-emissivity sources. For this case also, both heat -transfer 
rates are nearly the same. Therefore, only one heat-transfer rate curve is given. 

There is a significant difference between the two heat-transfer rates for the three 
glass-blasted shields (e ~ 0.27). As can be seen, the assumption of nonuniform radiosity 
gives better agreement with the experimental data. 

In figure 4, the experimental and analytic temperatures are plotted as a function of 
radial position for a three-shield configuration. The shields were made of Mylar which 
is a very low conducting material; and, consequently, a large radial temperature gradient 
is present in each of the shields. All of the surfaces were painted with a high-emissivity 
coating (e » 0. 94) in order to give a high heat-transfer rate to the colder source. A high 
heat-transfer rate was desirable in order to improve the accuracy of its measurement. 
The analytic radial temperature distribution closely follows the experimental one for each 
of the shields, even to the extent of predicting the temperatures at the outer edge of each 
of the shields. 

In figure 5, the radial temperature distribution is plotted for a two-shield configura- 
tion. In the experiment, all surfaces except one were coated with a high-emissivity paint 
(e « 0.94) . The side of the warmer shield facing the colder shield was left unpainted ex- 
cept for a rim which covered 0.16 of the radial distance. The purpose of this test was to 
see if the calculations would accurately predict the behavior of a shadow- shield system 
when the emissivity of a surface was not constant. The shield with a targeted rim was 
made of Mylar with a thin layer of aluminum laminated to the Mylar on each side. This 
resulted in an emissivity of about 0.03 for the unpainted center portion. The colder of 
the two shields was a plain Mylar shield painted on both sides with a high-emissivity 
paint. Even though the aluminum increased the conductivity of the warmer shield, the 
thickness of the aluminum was so small that the overall conductivity was still low. If the 
rim of the shield were not coated, one would expect a temperature distribution similar to 
that presented in the previous figure. However, coating the outside surface of the shield 
causes much more energy to be emitted from the surface. This energy causes a local in- 
crease in the temperature of the adjacent colder shield. This effect is illustrated in this 
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figure. The analytic radial temperature distribution closely follows the experimental dis- 
tribution for both the warmer and colder shields. 


Structural Member 

In reference 8, analytic data were presented for the temperature distribution and 
heat-transfer rates for a pipe penetration with an adiabatic external surface . Setting the 
external emissivity equal to zero simulated this condition when using the analysis pre- 
sented herein. For similar cases, results calculated using the equations presented in 
appendix D, agreed to within the accuracy of the plot reported in reference 8. 

In reference 7, experimental temperature distributions were reported for tubes de- 
signed as scaled-down representations of support members. Figure 6 gives two of these 
temperature distributions and an analytic comparison for each distribution. One of the 
struts was left unpainted while the other was partially painted in the circumferential di- 
rection to give an average emissivity of 0.6 in the axial direction. The manner in which 
the second strut was painted is shown in figure 6. The measured emissivity of the bare 
metal over the temperature range to which the struts were subjected is closer to 0.25 
than the handbook value of 0. 3 which was reported in reference 7. The analysis herein 
uses an emissivity of 0.25 for the bare metal and accounts for internal emissivity. The 
analytic curves in reference 7 were based on an emissivity of 0. 3 and did not account for 
the effects of internal emissivity. Both struts were made of stainless steel which has a 
relatively low thermal conductivity. When considering only conduction with a constant 
thermal conductivity, a straight line would connect the end points in figure 6. It can be 
seen from this figure that at the higher emissivity the temperature distribution departs 
more from the straight-line distribution. 

One of the possible reasons for the discrepancy between the analytic and experimen- 
tal results in figure 6 is that the analysis assumes the strut to be a cylinder while the ex- 
perimental struts formed a cone at each end. Also, the calculations were carried out as- 
suming that the bare metal inside surfaces of the tube were diffusely reflecting. In ref- 
erence 9, the effects of having specular surfaces on the heat transfer through a noncon- 
ducting tube are examined. In this reference, it is shown that for an 1/ d^ of 20 and an 
internal emissivity of 0.25, which is representative of the struts reported in reference 7, 
over four times as much energy can be transported through the tube when the surfaces 
are specularly reflecting as opposed to when they are diffusely reflecting. Since much of 
this energy would be reflected at the end of the tube, the temperature profile at the cold 
end would be altered if the surfaces were specularly reflecting. If, in addition, the ener- 
gy was polarized so that the component of reflectivity was large in the direction parallel 
to the axis of the tube, larger heat-transfer rates would result. 
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Consequences of Assuming Specularly Reflecting Surfaces 
and Uniform Radiosity 


The primary analysis used in this report (appendix B) assumes that all surfaces are 
diffusely emitting and reflecting and have nonuniform radiosity. In this section of the re- 
port, the change in the heat-transfer rate which results from surfaces being diffusely 
emitting but specularly reflecting is examined. Also, the effects on the heat-transfer 
rate due to the simplifying assumption that the surfaces have uniform radiosity are dis- 
cussed. It is simpler to calculate heat -transfer rates assuming surfaces of uniform radi- 
osity. However, this assumption does not yield temperature profiles for each shield and 
may significantly underestimate the real heat-transfer rates. 

Specular surfaces . - In appendix C, the method used to calculate the heat-transfer 
rates and temperature of each shield for specularly reflecting surfaces is presented. 

This method is for diffusely emitting surfaces under the assumption of uniform radiosity. 

An indication of the effects of assuming specular rather than diffuse reflections is 
given in figure 7 for various surface emissivities and spacing ratios. The ordinate gives 
the ratio of the heat-transfer rate to a source at zero temperature for specular surfaces 
of uniform radiosity to the heat-transfer rate for diffusely reflecting surfaces of uniform 
radiosity. Both surfaces have the same emissivity and curves of constant emissivity are 
plotted. It can be seen that as the spacing ratio increases at constant emissivity, the 
ratio of the heat -transfer rates also increases. In addition, as the emissivity decreases 
the ratio also increases. 

In figure 8, the same parameters as in figure 7 have been plotted except that a shield 
has been placed midway between the two sources. The surfaces of the shield are specu- 
larly reflecting like the sources and have the same emissivity as the sources. The re- 
sults of placing the shield between the sources is to magnify the ratio of the heat -transfer 
rates. In figure 8, the shield temperature is uniform due to assuming uniform radiosity 
and emittance. Because of the assumption of uniform radiosity, these two figures only 
give an indication of the effects of specular reflectivity. It would be more accurate to 
consider the effect of nonuniform radiosity on the heat-transfer rates. However, this 
was beyond the scope of this report. 

Comparison of uniform and nonuniform radiosity assumptions . - The two previous 
figures presented a comparison of the heat-transfer rates for specularly reflecting with 
diffusely reflecting surfaces under the assumption of uniform radiosity. In figures 9 
and 10, comparisons are made in the heat -transfer rates between uniform radiosity and 
nonuniform radiosity under the assumption of diffusely reflecting surfaces. In figure 9, 
the ratio of the heat -transfer rates of nonuniform to uniform radiosity is given as a func- 
tion of spacing ratio for constant emissivity. Both surfaces have the same area; and, 
again, the ratio is for a cold-source temperature of zero. It can be seen from figure 9 
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that as the spacing increases, the ratio of the heat-transfer rates reaches a maximum and 
then decreases. However, the ratio of the heat-transfer rates at any spacing decreases 
with increasing emissivity. 

In figures 10(a) and (b), the same parameters as in figure 9 are plotted except that a 
shield has been placed midway between the two sources. Figure 10(a) is for a noncon- 
ducting shield while figure 10(b) is for a uniform temperature shield. A comparison of 
these two figures shows that the effects of the shield conductivity assumptions on the heat- 
transfer ratio are small in relation to the effects of the radiosity assumptions. In each of 
these figures, the effect of placing a shield between the source is to magnify the ratio of 
the heat-transfer rates. 

In the previous four figures, the heat-transfer ratio was given as the dependent varia- 
ble. In the next two figures, the actual heat-transfer rate is given as a function of emis- 
sivity. In figure 11, the nondimensional heat-transfer rate is plotted as a function of 
emissivity for two different spacing ratios and the various types of surfaces which have 
been discussed. In figure 12, the same parameters have been plotted with the addition of 
a single shield midway between the two sources. 

At an L^/R of 0. 1 in either figure 11 or 12, the difference between the heat-transfer 
rates for diffuse surfaces of uniform radiosity and specular surfaces of uniform radiosity 
is small. The difference between diffuse surfaces of uniform radiosity and those of non- 
uniform radiosity is larger. At an Lrp/R of 1.0, there is a significant difference in the 
heat-transfer rates between specular and diffuse surfaces of uniform radiosity. How- 
ever, an L T /R of 1.0 represents a relatively large spacing for shadow-shield applica- 
tions. 

By comparing figures 11 and 12, it can be seen that as the number of surfaces in- 
volved increases, the differences in the heat-transfer rates caused by the various as- 
sumptions generally increase. The curves presented in this report for diffusely reflect- 
ing surfaces of nonuniform radiosity might be significantly different from the results ob- 
tained for specularly reflecting surfaces of nonuniform radiosity if the number of shields 
were large and the emissivity low with a large spacing ratio. 

In the previous six figures, the cold-source temperature was taken to be zero. When 
the temperature of the colder source is increased, the change in the heat-transfer rate 
caused by the different analyses becomes greater. This is illustrated in figure 13 where 
the increase in the heat-transfer rate is given for various source temperature ratios. 

The increase in the heat-transfer rate is due to assuming specularly reflecting surfaces 
over diffusely reflecting surfaces. The increase in the nondimensional heat transfer to 
the colder source in terms of the warm-source temperature ^Qj^aA^T^ is plotted as 
a function of spacing ratio for constant source temperature ratios and a system consisting 
of two plane sources with no shields. It is necessary to plot the increase in the heat- 
transfer rate rather than the ratio of the heat-transfer rates because for some cold- 
source .temperature ratios, the heat-transfer rate for diffusely reflecting surfaces goes 
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to zero. It can be seen from figure 13 that the cold-source temperature has to become a 
fairly large fraction of the warm-source temperature before the increased heat transfer 
due to specular reflectivity becomes much greater than the increased value for a cold- 
source temperature of zero. 

It is simpler to calculate heat-transfer relations assuming uniform radiosity sur- 
faces . If one does not require the temperature profile across the shield and the emissiv- 
ities of the surfaces involved are not too small, the heat-transfer relations presented in 
reference 1 for diffusely reflecting surfaces of uniform radiosity may suffice. The infor- 
mation presented in the previous figures should aid one in judging whether the assumption 
of uniform radiosity is satisfactory. 


PARAMETRIC DATA 
Shadow Shields 

4 

The equations used to determine the heat-transfer rates are linear in T when there 
is no conduction or when the shields have uniform temperatures. An additional require- 
ment for linearity is that surface properties be independent of temperature. In appen- 
dix B, it is shown that when the equations are linear, the heat-transfer rate for any pair 
of source temperatures can be found using only a pair of nondimensional heat-transfer 
rates. Each of these rates gives the contribution of one of the sources. The equation for 
the addition of the heat-transfer rates is: 



• * / 4 

Qj/Aj is the net heat-transfer rate per unit area to source 1. Qj/aA^Tg is the 
nondimensional heat-transfer rate to source 1 due solely to a temperature on source 2. 
Qj/ctAjTj is the nondimensional heat-transfer rate to 1 resulting from the temperature 
of source 1. T^ and Tg are the source temperatures. The application of this equation 
is illustrated in a subsequent portion of the report. 

One of the applications of shadow shields would be the reduction of the heat transfer 
to a liquid hydrogen tank. The emissive power of a source at liquid hydrogen tempera- 
ture is low. Therefore, much of the data in this report is for the heat-transfer rate to a 
source at 0°. This is done to avoid the necessity of giving two heat-transfer rates for 
each case. Nonconducting shields yield higher heat-transfer rates than conducting 
shields. For conservative results, heat-transfer rates are for nonconducting shields un- 
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less otherwise stated. For simplicity, the data given in this report are for shields of the 
same size. Planar sources have the same surface area as each side of the shields, and 
nonplanar sources have the same lateral radius as the shields. Also, for convenience, 
data are given for black surroundings at zero temperature. 

Source ge ometry . - The heat-transfer rate for two plane sources with and without a 
single shield placed between them has been given in figures 11 and 12. In figure 14, the 
heat-transfer rate to a source at zero temperature is given as a function of emissivity for 
two nonplanar sources and various spacings. Each source is an oblate spheroid whose 
major radius is |/iT times the minor radius. It can be seen from this figure that the 
heat-transfer rate is a strong function of emissivity. As the emissivity decreases, the 
effects of the spacing ratio become somewhat more pronounced. 

In figure 15, the parameters given in figure 14 are plotted for a single planar shield 
placed midway between the two nonplanar surfaces. With the addition of a single shield, 
the effects on the heat-transfer rate due to the overall spacing become more pronounced 
especially at lower emissivities. 

Number of s hield s. - The next two sets of figures are used to determine the net heat- 
transfer rate to a plane source when both sources have nonzero temperatures. Figures 
16(a) to (d) give the heat-transfer rate to a source of zero temperature, and figures 
17(a) to (d) give the heat-transfer rate when the other source has a zero temperature. 

In figures 16(a) to (d), the heat-transfer rate to source 1 is given as a function of the 
number of shields placed between the two sources. For each of the figures, the overall 
spacing ratio (Lrp/R) is held constant and curves are given in which all of the surfaces 
have the same emissivity. It should be noted that figures 16(a) and (b) have one ordinate 
scale while figures 16(c) and (d) have another scale. Since a fractional number of shields 
has no physical significance, only integer values of the abscissa have meaning. By com- 
paring particular points on each of these figures, it can be seen that increasing the spac- 
ing ratio increases the efficiency of shadow shields. 

Figures 17(a) to (d) are the complement of figures 16(a) to (d) in that the heat- 
transfer rate to source 1 due to the presence of a temperature on source 1 is given. 

Again, note that there is a scale change in the ordinate of these figures. Figures 17(a) 
and (b) have one scale while figures 17(c) and (d) have another scale. There is a net heat 
transfer out of the source for all cases; therefore, the ordinate of figures 17(a) to (d) is 
negative. It can be seen from these figures that as the spacing between the source and 
nearest surface becomes large, the heat-transfer rate out of the source approaches the 
limiting value of Qj/crAjTj = e. This quantity is the amount of heat which would leave 
the source if the nearest surface was at an infinite distance from the source. 

In order to show how the data in figures 16 and 17 may be used to determine the net 
heat transfer to a source, an illustrative example is given. The example is for two plane 
sources with two equally spaced shields. A spacing ratio (L,p/R) of 0. 3 is chosen and the 
emissivity of each surface is 0.1. The warmer source has a temperature of 300 K 
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(540° R). The colder source has a temperature of 100 K (180° R). 

The first nondimensional heat-transfer rate in equation (1) (Qj / aA jTj) is found from 

figure 16(c). The second nondimensional heat -transfer rate (Qi/^jT 4 ) is found from 
figure 17(c). Substituting the appropriate values in equation (1) gives the following result: 


Q 

A 


1= 5.67X10 8 [(2.1X10 3 )(300) 4 + (-0. 065)100 4 ] 


V 0 .595 ^ = M89_Btu 

A 1 m 2 (hr) (ft 2 ) 

Figure 18 illustrates the effects of varying the temperature of the colder source on the 
heat-transfer rate. In this figure, the nondimensional heat-transfer rate in terms of the 
temperature of the warmer source is plotted as a function of emissivity for various tem- 
perature ratios. It can be seen from this figure that the colder temperature has to be a 
fairly large fraction of the warmer temperature before the heat-transfer rate changes a 
great deal with respect to temperature. It can also be seen from this figure that for some 
temperature ratios the minimum heat-transfer rate has a negative value. This means 
that the absolute minimum heat-transfer rate may occur at an emissivity other than zero . 
When all surfaces have the same emissivity, the minimum heat-transfer rate is not a 
strong function of emissivity. If the emissivity were allowed to vary from one surface to 
another, the minimum could be lower. This can be seen from the dashed curve in fig- 
ure 18. For this curve, the colder source has an emissivity of 0.9 while all other sur- 
faces have the emissivity given by the abscissa of figure 18. The temperature ratio 
(Tj/Tg) is 0.75. 

Up to this point, only heat-transfer rates have been discussed. Shield temperatures 
are also of interest since the temperature gradients in the shields may affect the design 
of the shield system. The following figures give temperature profiles for some represen- 
tative nonconducting shields. Figures 19(a) to (c) present the shield temperature against 
radial position for a one-, two-, and three-shield case. For each of the three figures, 
the overall spacing between the sources is held constant. It can be seen that the shield 
temperatures drop off sharply toward the outer edge of each of the shields. The temper- 
ature difference across the shield become slightly greater as the temperature of the 
shield increases. The short horizontal line crossing each of the curves indicates the 
temperature of an infinitely conducting shield under the same conditions. 

Targeting . - When the emissivity of the surfaces of a shadow-shield system are uni- 
form in the radial direction, reducing the emissivity generally reduces the heat -transfer 
rate to the colder source. It is sometimes possible to gain a further reduction in the 
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heat-transfer rate for low-emissivity surfaces by using a high-emissivity coating on the 
outer edges. Shields for which this is done take on the appearance of targets due to cir- 
cumferential symmetry. It is generally necessary that there be more than one shield in 
the system for targeting to produce a reduction in the heat-transfer rate. 

Figures 20(a) and (b) illustrate the effects of targeting for uniform temperature 
shields and nonconducting shields for a system consisting of two plane sources separated 
by two shields. In order to obtain a decreased heat-transfer rate, it is necessary to tar- 
get the surface of the warmer shield facing the colder shield. The percentage change in 
the heat-transfer rate due to targeting is plotted as a function of the percentage area of 
the shield which is targeted. The results for two overall spacing ratios are shown. From 
figures 20(a) and (b), it can be seen that for this case there is some decrease in the heat- 
transfer rate due to targeting uniform temperature shields, but no improvement is notice- 
able for nonconducting shields . The targeted areas of the uniform temperature shields 
act like fins in that they are able to take energy from the shield and dissipate it to space. 
The nonconducting shields are less efficient at doing this. The changing of the shield 
temperatures due to targeting is illustrated in figures 21(a) to (c) for nonconducting 
shields. In these figures, the temperature profiles are plotted as a function of radial po- 
sition for an untargeted system and two combinations of targeted shields. The upper 
curve gives the temperature of the warmer shield in each figure. In each case, the non- 
dimensional heat-transfer rate is also given. Even if targeting does not provide an im- 
provement in the heat transfer for nonconducting shields, targeting can drastically alter 
the temperature profile of the shield. 

The heat-transfer rate varies little while the temperatures vary substantially in each 
of the three cases presented. In figure 21(b), the temperature of the outer edge is 
greatly lowered. This could prove useful in attaching shields to their support members. 
Shields could be thermally bonded to the support structure. If the rim temperature was 
lower than that of the adjacent structure, the shields could act as fins to dissipate heat 
from the strut. In this way, the local temperature of the strut would be lowered and this, 
in turn, might reduce the conducted heat out of the strut. The interaction of shields and 
their support structure is not analyzed in this report. In figure 21(c), the radial temper- 
ature difference has been lessened for each of the shields by means of targeting. 

Shield co nductivity . - Figure 22 illustrates the effect of shield conductivity on the 
heat-transfer rate. As the conductivity of a shield increases, the temperature of the 
shield becomes more uniform and the heat-transfer rate decreases. Shields which have 
the same values of the scaling parameter (kt/oTgRp where T2 is a reference tempera- 
ture) have the same radial temperature profiles and heat-transfer rates. In figure 22, 
heat-transfer ratios for partially conducting shields are given as a function of this param- 
eter. The ordinate of figure 22 is the ratio of the heat-transfer rate to the heat-transfer 
rate for uniform temperature shields and data are given for one-, two-, and three -shield 
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configurations. Only when the conductivity parameter has a value of zero or infinity are 
the equations for the heat-transfer rates linear. It should be noted that at these two ex- 
tremes the abscissa has been made discontinuous in order to show the effects of the two 
assumptions which yield linear equations. From this figure, it can be seen that when 
there are several shields involved, there can be a significant change in the heat-transfer 
rate due to thermal conductivity. These curves are given for an emissivity of 0.1 and 
although it is not shown on the figures, the ratio of the heat-transfer rates increases with 
decreasing emissivity. 

When the thermal conductivity has a finite value other than zero, it is necessary to 
solve nonlinear equations in order to determine the temperature distributions. In appen- 
dix B, a solution is given in which the conducted heat transfer is treated as a known quan- 
tity during each iteration cycle. As the conductivity of the shields increases, it becomes 
less desirable to do this. In appendix E, an alternate solution is given in which the ra- 
diant heat transfer is treated as a known quantity during each iteration, and the tempera- 
ture distributions are found by solving the conduction equations. 

Shield posit ion. - In the previous figures, the shields were placed so that there were 
equal spacings between surfaces. Figures 23(a) to (e) illustrate the effects of varying the 
spacings between surfaces for two emissivities. In each of the figures, the variable giv- 
ing the relative spacing (L/Lrp) is measured from the warmer source. 

The ordinate in each figure is the ratio of the heat-transfer rate to the heat-transfer 
rate for the same number of evenly spaced shields. In figure 23(a), there is a single 
shield present and in figures 23(b) and (c), there are two and three shields, respectively, 
with the spacing between shields remaining constant. For these three cases, the spacing 
for the minimum heat-transfer rate is halfway between the two sources for both emissivi- 
ties. As the emissivity decreases, the heat-transfer rate is more sensitive to spacing. 
The heat -transfer rate decreases with increasing number of shields, so that the ordinate 
is different for each of these figures. The actual heat -transfer rate is lowest for three 
shields. 

Figures 23(d) and (e) are for configurations in which the shields have the same spac- 
ing relative to the sources. In figure 23(d), both shields come closer to the center as the 
abscissa is increased; and in figure 23(e), the center shield remains fixed while the two 
outer shields approach it for increasing values of the abscissa. In both of these figures, 
the spacing for minimum heat transfer is a function of emissivity. 

Solar radiat ion. - When one of the sources is solar radiation, this radiation is inci- 
dent only on the outermost surface. If the equations are linear, the net heat-transfer 
rate to a source can be found by using two generalized heat -transfer rates. The heat- 
transfer rate due to solar radiation is proportional to the solar flux (cp) times the solar 
absorptivity ( a ) . The solar absorptivity is taken to be independent of the emissivity of 
the surface. In this way, increasing the solar absorptivity is equivalent to increasing the 
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solar flux. Therefore, these variables can be combined with others to form nondimen- 
sional heat- transfer rates. 

Figure 24 gives the nondimentional heat-transfer rate due to solar energy as a func- 
tion of the spacing ratio for two nonconducting shields and a plane source. All surfaces 
have the same emissivity, and lines of constant emissivity are plotted. It can be seen 
that as the spacing ratio decreases, the heat- transfer rate increases with the greatest 
change coming for low-emissivity surfaces. As the spacing ratio reaches a very small 
value, the heat-transfer rates become less dependent on the emissivity of the surfaces. 

In figure 25, the nondimensional heat-transfer rate is plotted as a function of spacing 
ratio for surfaces of constant emissivity. In this figure, the solar flux is taken to be 
zero. The inner shield is placed midway between the source and the outer shield. The 
heat-transfer rate to the plane source is negative since there is no solar flux. From this 
figure, it can be seen that as the emissivity decreases and spacing ratio increases, the 
heat-transfer rate approaches the limiting value of QjyoAjTj = e. 

The information presented in figures 24 and 25 can be utilized in order to determine 
the net heat-transfer rate into a source when the system is exposed to solar radiation. 
The solar flux varies with the square of the distance from the sun. At one AU, the solar 
flux is about 1390 W/m 2 (442 Btu/(hr)(ft 2 )) . If it is desired to approach a one-tenth AU 
distance from the sun, the solar flux would be 1.390x10^ W/m 2 (4.42x10^ Btu/(hr)(ft 2 )) . 
The sun is not really a point source so that a conical array of shields might be needed. 
However, at this distance from the sun, the half angle of the cone is less than 3° and this 
complication is neglected. In order to use figures 24 and 25, a system of two equally 
spaced nonconducting shields is chosen. The overall spacing is taken to be 0.2 and the 
source temperature is taken as 300 K (540° R). The emissivity of all of the surface is 
assumed as 0. 1 while the solar absorptivity is taken as 0.02. The net heat-transfer rate 
to the source can be found from the equation: 



The first heat-transfer term is found from figure 24 and the second from figure 25. Sub- 
stitution of values gives 


Q 


1 = (0. 0545)1. 390X10 5 (0. 02) + (-0. 064) (5 . 67xl0“ 8 )(300) 4 


= 122.6 =38.9 -™L 


m 


(hr) (ft 2 ) 
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This value of 122.6 W/m^ is considerably less than the value of 2734 W/m^ (867 Bin/ 
o 

(hr) (ft )) which would occur if no shields were present. 


Structural Member 

The structural member which is thermally analyzed in this report could be a support 
member of a system for shadow shields. It could also be any of several other structural 
members in an upper stage. The thermal analysis of this member can involve many in- 
dependent variables. This is especially true when the internal emissivity is not zero. 
The following figures present datafor a few cases representative of structural members 
in a space vehicle. In figures 26(a) to (c), the nondimensional total heat-transfer rate is 
given as a function of emissivity. The solid lines are for closed ends while the dashed 
line in each figure is for opened ends and an external emissivity of zero. When the ends 
are closed, the total heat-transfer rate is the conducted heat-transfer rate at the colder 
end plus the radiant energy absorbed by the end disk. When the ends are open, the total 
heat-transfer rate is the conducted heat -transfer rate. When the internal emissivity is 
zero, the internal radiant heat-transfer rate is also zero. For each of these figures, a 

different 1 / d is given and since the diameter of the tube is held constant, the radiation 
i 3 2/ 

to conduction parameter crT^ j kt increases with increasing 1/ d^. It should be noted 
that figure 26(c) has a different ordinate scale from figures 26(a) and (b). 

From these figures, it can be seen that both the internal and external emissivities 
can affect the heat-transfer rate. The heat-transfer rate is more sensitive to changes in 
emissivity at low values of emissivity than at high values. As the external emissivity in- 
creases, more energy is radiated to the surroundings. This, in turn, reduces the tem- 
perature gradient at the colder end of the strut thereby lowering the conducted heat- 
transfer rate. In these cases, the surroundings were taken to be at zero temperature. 
Increasing the internal emissivity causes an increase in the heat-transfer rate at the 
colder end of the strut in a twofold manner. First, as the internal emissivity increases, 
the temperatures along the strut increase. This increases the thermal gradient at the 
colder end which results in an increased conducted heat-transfer rate. Second, for 
closed ends, increasing the internal emissivity increases the amount of radiant energy 
which is absorbed by the end of the strut. 

In the next series of figures, the temperature profiles are given as a function of ex- 
ternal emissivity for constant values of 1/ d^ . In figures 27(a) to (c), the nondimensional 
temperature is plotted as a function of axial position. For all of these figures, the in- 
ternal emissivity is held constant at 0.6 while the l/d f changes for each figure. By 
comparing these figures, it can be seen that as the radiation to conduction parameter in- 
creases, the temperature profile departs more from the straight line which would occur 
if conduction were the only mechanism determining the heat-transfer rate. In each of 
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these figures, increasing the external emissivity causes an increase in the thermal gra- 
dient at the warm end of the strut and a decrease in the thermal gradient at the cold end. 
It should be noted that the thermal conductivity varies greatly between materials. There- 
fore, the values of the radiation to conduction parameter are strongly dependent on the 
strut material used. The values of this parameter which are shown on figure 27 were 
formed using a thermal conductivity representative of fiberglass. 

Generally, it is desired to minimize the heat-transfer rate to the colder end of the 
strut. Returning to figures 26(a) to (c), it can be seen that reducing the internal emis- 
sivity of the strut reduces the heat-transfer rate. If the inside of the tube were filled 
with a very low conductivity material, the effective internal emissivity would be zero and 
the added heat transfer due to the conduction of the filler would be small. Also, if radia- 
tion barriers were placed along the length of the tube, the strut would be divided into a 
series of nearly isothermal compartments and the effective internal emissivity would be 
reduced. 

In appendix D, it is shown that when the internal emissivity is zero and when emis- 
sivity and conductivity are taken to be independent of temperature, the temperature pro- 
file, and, consequently, the heat -transfer rate for a thin-walled tube seeing a zero- 
temperature environment can be expressed by two parameters. One of these parameters 
is the ratio of the boundary temperatures, and the other contains various specifications 
of the strut. The second parameter ( k ) is ere ^ 2 ^- / kt and is dimensionless . 

Figure 28 presents the ratio of the temperature slope at the colder end of the strut to 
that obtained by considering conduction only. This is done for an internal emissivity of 
zero. The data are presented as a function of k and various temperature ratios. The 
ordinate of figure 28 can be directly related to the heat transfer out of the tube. The 
heat-transfer rate out of the tube considering the effects of emissivity is k?rd 0 t(dT/dx) 
while the heat-transfer rate, when only conduction is considered, is kird^Tg - T ^)J l . 
Increasing k decreases the heat-transfer rate. The effects of this parameter become 
more pronounced as the temperature ratio decreases. When the ordinate of this figure 
becomes negative, heat is being transferred into the strut from both ends and is dissipa- 
ted to the surroundings along the length of the tube. 

In figures 29(a) to (c), temperature ratio is plotted as a function of axial position for 
various values of x. In each of these figures, the temperature ratio is held constant; 
and it can be seen that as x increases, the temperature profiles become more dis- 
torted. As x increases, the temperature along the strut becomes lower . If the tem- 
perature ratio is sufficiently low and x is sufficiently high, the temperature along the 
strut can fall below the colder end temperature. Under this circumstance, heat is trans- 
ferred into the strut at both ends. 
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CONCLUSIONS 


In this report, the thermal analyses of both a system of shadow shields and a struc- 
tural member have been carried out. It has been shown that the shadow-shield analysis 
developed herein is capable of accurately predicting the behavior of shadow-shield sys- 
tems. The effects of emissivity and spacing ratio on the heat-transfer rates have been 
shown. This has been done for a system consisting of up to several shadow shields. This 
report shows that the effectiveness of shadow shields is most pronounced when the emis- 
sivity of the surface is low. By targeting shadow shields, large changes in the tempera- 
ture profile of the shields can be achieved. This can be done with only a moderate change 
in the heat-transfer rate and may be advantageous in reducing the heat transferred by 
conduction through the strut if the shields are thermally connected to the support struc- 
ture. In some cases, the heat-transfer rate will be lowered when targeting is applied. 

In addition to protecting a cryogenic propellant, shadow shields may be used to pro- 
tect a payload or other body from incident solar energy. Here again, shadow shields can 
be effective in reducing the heat-transfer rate especially when the emissivities of the sur- 
faces are low. 

The effects of having specularly reflecting surfaces have been examined for the sim- 
plifying case of surfaces with uniform radiosity. It has been shown that with this as- 
sumption specularly reflecting surfaces yield a higher heat-transfer rate than diffusely 
reflecting surfaces. The percentage increase due to specular reflectance increases with 
increasing shield spacing. 

The effects of nonuniform radiosity for diffusely reflecting surfaces have also been 
examined. It has been shown that for reasonable shield spacings the possible error due 
to assuming uniform radiosity surfaces can be significant. This error is often greater 
than the error which would result from treating specularly reflecting surfaces as dif- 
fusely reflecting surfaces at the same spacing. For a more complete understanding of 
the behavior of specular surfaces, an analysis which would account for both specular re- 
flectivity and nonuniform radiosity is necessary. 

The ability of the thermal analysis for a structural member to predict the behavior 
of an actual strut has been examined. By taking advantage of the external emissivity of 
the strut, the rate of heat transfer through the strut can be greatly affected, even to the 
point of becoming a negative value. The thermal interaction between the shadow shields 
and their support members would also be needed for a more complete understanding of 
the entire system. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 3, 1968, 

124-09-05-12-22. 
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APPENDIX A 


SYMBOLS 


A area of a radiating surface 

a cross-sectional area of an element on the strut 

x.i_ 

B,. absorption factor - the fraction of energy which leaves the 1 n surface and is 

l J th 

absorbed by the j surface 

d diameter of the strut 

F view factor between specular surfaces 

f.j diffuse view factor between the i Ln and j surfaces 

H incident energy per unit time and area 

h step size 

K constant used in appendix B 

k thermal conductivity 

L spacing between surfaces in the shield system 

l length of the strut 

m shield number or the number of steps in the solution of the differential equation 

n number of annuli on a shield or the number of elements into which a strut is 

divided. This symbol is also used as a subscript. 

Q heat-transfer rate to a whole surface 

q heat -transfer rate to an element 

# thermal resistance between annuli 

R radius of a shield 

r radial distance 

T absolute temperature 

t shield thickness or thickness of a thin-walled tube 

W emissive power per unit area of an element 

x distance along the strut 

a absorptivity 

solar absorptivity 
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0 


? 

V 

9 

K 

£ 
7 T 


angle of incidence or reflection 

Kronecker delta; equals 1 when i = j; equals 0 if i ^ j 
emissivity 

number of terms used to approximate the finite series 

number of times energy has been reflected off an adjacent surface 

nondimensional temperature along the strut 


radiation to conduction parameter for the strut, 


ejh 2 

O A 
kt 


transformed independent variable in the strut equation 
a constant 


p reflectivity, p = 1 - a 

a Stefan-Boltzmann constant 


<p solar flux 

The following matrices and vectors are used: 

C k-by-k matrix whose elements are the fraction of energy which leaves the q 
surface and is reflected by the p n surface (C) = f qp p p - 5 pq 

D n-by-n matrix whose elements are the sum of the entries in and 'i r which 
refer to the annuli of the shield times the emissivity of the annuli 

(°)pq - ^n+p, n+q e n+p + ^ ^pq e p 

E column vector of size k whose elements are the emissive power of each surface 
<E)p = = «p t £a p 

I identity matrix (I) = 6 pq 

P column vector of size n with entries containing the radiant heat transfer from 

each element of the shield from known temperature surfaces 

Q column vector of size n whose elements contain the heat transfer away from each 

f h 

annulus of the m shield 


U 


u 
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k-by-k matrix whose elements are the absorption factors (U) pi = B pi in the shield 
system 


column vector of U in the shield system (u) p = B p ^ and a vector containing the 

absorption factors divided by the area of the differential element in the strut 
B.. 

system (u)^ = — 

A j 


J 



r 


V k-by-k matrix whose elements contain the negative of fraction of energy absorbed 

by a surface (V) pq = -f pq « q 

v column vector of V in the shadow-shield system (v) = -f q j or j and a vector 

containing the negative of the fraction of energy absorbed by a surface in the 

strut system (vl = -f. — i- 
P JP A 
A P 

Y n-by-n matrix containing the coefficients of the temperatures when the equations 

are linear in temperature 

4 

Z column vector of size n containing (Z)^ = crA^T^ 

0 column vector containing the temperatures of the annuli of the shield or the tem- 

peratures along the strut 

$ column vector of size k whose entries contain the heat-transfer rate to each 
annulus in the system 

X m-by-m matrix whose entries contain the temperature coefficients in the solution 

to the differential equation 

\k k-by-k matrix defined as 'k = U - I 

£2 column vector of size m, the entries of which are the values of the second de- 

rivative in the differential equation 

Subscripts: 

c conductive heat transfer 

1 any surface in the system 

in heat transfer towards the surface 

j particular surface under consideration 

k total number of elements in a system. Equals 2n + 1 for a shield system and 
n + 2 for a strut system . 

n number of annuli on a shield or segments of a strut 

o outside surface of strut 

out heat transfer away from the surface 

p,q dummy subscripts denoting any element in a matrix or vector 
r radiant heat transfer 

s shield 

sr surroundings 
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T total heat -transfer rate through the strut or the total spacing between sources 

t inside surface of the strut 

1, 2 each of the two sources or each end of the strut 

Superscripts: 

XV 11 

c complimentary system composed of the m L and m in + 1 shields and their 

surroundings 

source for which the heat-transfer rate is given has a zero temperature 

* source for which the heat-transfer rate is given is the only source of thermal 

energy in the system 

+ heat -transfer rate is a ratio to the case of uniform radiosity and diffusely reflect- 

ing surfaces 
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APPENDIX B 


DERIVATION OF SHADOW-SHIELD EQUATIONS 

The following derivation for determining the temperature distribution for a system of 
shadow shields follows the work of Gebhart given in reference 10. Heat balances are for- 
mulated which consider radiation and radial conduction. The steady -state temperature 
distribution for each shield is determined from the heat balances, since there is no net 
heat transfer to the elements of the shield or the shield as a whole. The properties of the 
shields are allowed to be functions of the shield temperatures. This results in an itera- 
tive solution. In order to avoid matrices of cumbersome size, the temperature distribu- 
tion for each shield is determined assuming that adjacent shields and surroundings have 
known temperature distributions. As the iterative procedure converges, the assumed and 
calculated temperature distributions become nearly identical. 

Figure 30 gives a schematic of a shield bounded by two other surfaces and the sur- 

i u 

roundings. In the derivation of the equations for the temperature of the m l shield, it is 
assumed that the temperatures of the m^ + 1 and the m^ 1 - 1 shield are known. Each 
of the four surfaces (there are two surfaces to consider for the m L shield) are divided 
into n elements. Because surface properties are assumed constant in the circumferen- 

T Pi 

tial direction, the elements are a series of concentric annuli. Both surfaces of the j 
element of a shield have the same temperature and area, but not necessarily the same 
surface properties. The temperature of both sides of the shield are the same because it 
is assumed that there is no axial temperature gradient. 

The total rate of heat transfer to the j^ 1 element on the m*"* 1 shield is the sum of 
the radiant heat transfer occasioned by the shield seeing the m^ + 1 and m*"* 1 - 1 
shields plus the conducted heat transfer from the j^ - 1 and the j*"* 1 + 1 elements of the 
nr shield to the j element. The rate of heat transfer away from the element is the 
emitted radiation from both surfaces of the element. 

The shield temperatures are determined by constructing two closed systems as is 
illustrated in figure 31. The first system consists of the m L - 1 shield, the m l 
shield, and the surroundings enclosing these two shields. The second system consists of 
the m^ 11 shield, the m^ 11 + 1 shield, and the surroundings enclosing these two shields. 
The m L shield forms the common boundary for the two systems. The temperature of 
the boundary must be the same when viewed from either system. If the enclosures do not 
represent physical surfaces, they are nonreflective. 

There are n elements on each surface of each system. Therefore, there are 
2n + 1 elements in each system. The manner in which elements were assigned sub- 
scripts is also illustrated in figure 31 . Each element in a system is assigned a unique 

f Pi 

subscript; therefore, it is not necessary to distinguish elements on the m shield from 
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those on the other surface with an additional subscript denoting the particular shield. 

The net rate of radiant heat transfer away from the j element is: 

q r,i =W j A r B l) W l A l B 2j W 2 A 2- B 3i W 3 A 3 ’ • • • 'VVk where k = 2n+1 (B« 

W. is the emissive power per unit area of the i in element, and B.. represents the 

1 1J 

fraction of energy which is emitted from the i m surface and is absorbed by the j sur- 
face. It is assumed that the radiant energy is uniformly distributed over each of the ele- 
mental surfaces, and that each of the surfaces emits and reflects energy diffusely. There- 
fore, B^ can be found in terms of the other B's. The following equations give this rela- 
tion. 


B lj = f lj“j + f ll p l B lj + f 12 p 2 B 2j + • • * + f lkPk B kj 
B 2j = f 2j“j + f 21^1 B lj + f 22 p 2 B 2j + • • • + f 2k p k B kj 


B kj = Vj + f kl p l B lj + f k2 p 2 B 2j + • • • + f kk P k B kj 
This set of linear algebraic equations can be rewritten as: 

( p l f ll " ^lj +p 2 f 12 B 2j + • • • +p k f lk B kj = 

p l f 2 1 B 1 j + ( p 2 f 22 " 1 ) B 2j + • • • +p k f 2k B kj = _f 2j a j 

p l f kl B lj +p 2 f k2 B 2j + ' ■ • + ( p k f kk " X ) B kj = “ f kj“j 
These equations can be rewritten in matrix form for each j. This gives: 

Cu = v (B2) 
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where 


( C )qp f qp P p 5 qp 


(u) 


P 



and 


The heat-transfer balance can be taken for each of the other elements in the system. 
The coefficient matrix C is independent of the surface for which the heat balance is taken. 
The entire k set of algebraic equations can be represented in terms of k-by-k matrices. 
The matrix equation for this is: 


CU = V 


(B3) 


Each of the matrices have the following values in their respective elements: 


= B pq 


(v) = -f a 

v 'pq pq q 


and the coefficient matrix C remains the same. 
Solving for U yields : 


U = C -1 V 


(B4) 


The entire process can be repeated for the other system. And an equation similar to 
equation (B4) can be written for this system. 


U c = 



(B4') 


The surface properties and the spacing between surfaces are not necessarily the same in 
each system. Therefore, there is no necessary correspondence between the systems. 

In general, 
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C ^ C c 

u^u c 

v c 

Equation (Bl) can be written for each of the elements in each of the systems. For the 

fh fh 

system containing the rn shield and the m - 1 shield, these equations are: 

-4r, 1 = ( B ll - !) W 1 A 1 + B 21 W 2 A 2 + • • • + B kl W k A k 
_q r, 2 = B 12 W 1 A 1 + ( B 22 " 1 ) W 2 A 2 + • • + B k2 W k A k 


- q r,k = B lk W l A l + B 2k W 2 A 2 + • + ( B kk - l ) W k A k 

These equations can be written in matrix form giving: 

$ = * T E (B5) 

Because the indices are reversed, it is the transpose of the coefficient matrix which is 
used in the matrix equation. $ is a column vector whose elements represent the net 
radiant heat transfer to each of the surfaces in the system. Each of these matrices are 
made up of the following elements: 


<*>pi = <%t - 6 pi or = U - I 

<*>1 * A 

< E >p = W p\ = q£ p T p A P 

i.r J.U 

For the other system involving the rn and nr — 1 shields, it is possible to 
write the heat-balance equations. These equations in matrix form are: 


$ 


c 



(B5’) 
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The m tla shield is common to both systems, and the temperature of each annulus of 

this shield is the same when viewed from either system (A, = A? and T. = T?, where i 

tjh ^ ^ 11 

represents any element of the m shield). However, because the surface properties as 
well as the spacing between surfaces may be different in each system, there is, in gen- 
eral, no correspondence between equations (B5) and (B5') so that 

<f> c t <f> 

* c ^ * 

E c ^ E 

It is desirable to maintain a system of linear equations. Radiation terms in the heat 

4 

balances are linear in T , while conduction terms are linear in T. In this appendix, 
equations are given which are linear in T . The conduction terms are considered known 
during an interation and are recalculated after every iteration. When the conduction 
terms are relatively large compared with the radiation terms, it may be advantageous to 
solve a system of equations that are linear in T. In appendix E, equations, which are 
linear in T, are presented for determining the temperatures of partially conducting 
shields. Either form of the solution should yield the same result since each is iterative 
in nature so that at the final iteration the terms that were assumed known are truly 
known. The heat-transfer rate to the j n element by conduction is 


_Vi 


T. 

1 


*c,J 


3t 


j+l,j 


\ , - T. 

Izl L 


dt is the thermal resistance between the 


,th 


element and an adjacent element 



where i = j ± 1 


(B6) 


(B7) 


In each of the systems considered, only n of the k elements are for the m^ 
shield. It can be seen from figure 31 that for the system containing the m 4 “ and 
m** 1 - 1 shields the elements for the m^ shield are the n^ 11 + 1 to the 2n elements. 
For all of the other 2k - 2n elements in the two systems, temperatures as well as emis- 
sivities and, therefore, reflectivities are considered known. 
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Equations (B5) and (B5’) can be added together and rewritten so as to separate the 
knowns from the unknowns. Our interests are in determining the temperatures of the 2n 
elements common to both systems. The algebraic equations which result from the addi- 
tion of equations (B5) and (B5') are: 


(4+1 + (* C )l - X (4,n+l( J 

V X (* c ) p ,i( ec ) p 

+ (4,n+l(4 

P=1 

p=n+l 


k-1 

n 


+ X (4.n+l( I 

p=n+l 

4 + X (*%>% 

p=i 


(4+2 + (*°) 2 - X (4,n+ 2 ( E )p + X (*%.2( e % + (4,n+ 2 ( E )k 

p=l p=n+l 


k-1 


n 


+ X 

p=n+l p=l 


V (B8) 


Wan + (4n = X (4, 2n( E )p + X W)p, n( E % + (4. 2n( E ) 

p=l p=n+l 


k 


k-1 


n 


+ X (4,2n( E )p + XW)p,„( EC ), 

p=n+l p=l 


It should be noted that in the equations of (B8), because of the order in which the sub- 
scripts were written, the transpose of ^ or was not used. 

For each of the algebraic equations of (B8), the first three terms to the right of the 
equal sign are for elements with known temperatures. Since there is no net heat transfer 
to any annulus of the shield, the two terms on the left side of each equation equal the rate 
of heat transfer away from the annulus by conduction 
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Wn+l+^V-Ic.i 



If the shield does not have a uniform temperature, it is convenient to write equa- 
tion (B8) in matrix form. The algebraic equations become: 

Q = D T Z (B9) 

m 

Equation (B9) is composed of an n-by-n matrix D and two vectors, each containing n 
values. (Q)^ represents the heat transfer away from the i^ annuli of the m^ shield 
due to conduction less the radiant heat transfer to the annuli from all of the surfaces 
whose temperatures are known 



(*)p,n«( E )p- X 
p=n+l 



(Z)j is an entry of a column vector and has the value 
trix containing the entries: 


aA.T?. 
J 1 


D is an n-by-n ma- 



n+j , n+i e n+j + 



The transpose of D can be written in terms of the matrices containing the absorption 
factors. This yields: 


( D )i,j ( U )n+j,n+i " 5 j,i_ c n+j + ( U )j,i 6 j,i_ e ] 


In these equations, e n+ j and e? both refer to the j th annulus of the m th shield, 
but each must be taken from the appropriate system. Aj and T. are the same in both 
systems. 

When one of the sources is solar radiation, the energy which is absorbed on the out- 

J-U 

ward facing surface of the j element is <pAjO!^ where cp is the incident solar flux, 
and ciy is the solar absorptivity. For flat disk shields, Aj equals the area of the annu- 
lus. However, if solar energy is incident upon a curved surface Aj, in this instance 
only, represents the area of the projection of the surface element into a plane perpendicu- 
lar to the solar flux. 

iL 

The heat transfer away from the j Ln annulus when a solar flux is present is: 



^ 1 % - X (* c )p,i( EC )p 

p=n+l 
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Also, when solar radiation is present, the outer shield is not bounded by either surround- 
ings or another surface. Therefore, none of the energy which leaves the surface in the 
direction of the solar flux is returned to any annulus of the shield. Under this circum- 
stance, the matrix D has the entries: 



-6. .e . 

J,i n+j 


+ 



Equation (B9) can be solved for the vector Z to give: 

z = (d t ) *q 

The temperatures of each ai llus can be found from the various values of Z 



(BIO) 


When the temperature of the shield is uniform, equation (B8) can be utilized directly 
to find the shield temperature . Since there is no net heat transfer to the shield as a 
whole, the following equation can be written: 


£ [(*)„« ♦ (* c )i] = 0 
i=l 

The sum of the equations given by equation (B8) can be written as: 

^(*)p,n + i( E )p + 2 
p=l p=n+l 




(Bll) 


S 


k-1 


Z Wi,n»rZMi,i( EC )i 

j=n+l j=l 


(B12) 


All of the quantities on the left side of this equation are known. The shield tempera- 
ture can be factored out of the right side of the equation. This yields: 
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(B13) 


T 


4 

m 



Z (*)p, n«( E )p + Z (*%, i( E % + Wk, n«( E )k 

p=l p=n+l 



When there is solar radiation present, the temperature of an infinitely conducting 
shield can be found from the amount of solar radiation absorbed by the shield as a whole. 
The shield temperature is found from the following equation: 


T 


4 

m 





(B14) 


When the conducted heat transfer to an element (q •) is zero so that the shield is 
either nonconducting or the shield has a uniform temperature, and all properties are in- 
dependent of temperature, the equations for determining the temperature are linear in 
T^. When the equations are linear, the shield temperature and heat-transfer rates can be 
found by summing up the contributions of each source separately. 

If there is no heat transfer between the elements of the shields, entries of the vector 

Q in equation (B9) can be portioned into two parts with each part being the negative of the 

heat transfer to the element due to one of the sources. Since the matrix D is a function 

only of the surface properties of the sources and not their temperatures, the values in the 

vector Z are proportional to the magnitude of the entries of Q. The entries of Z are 
4 4 

crAjT/\ T for every element can be found by summing the contribution due to each 
source. 

Similarly, the heat-transfer rates due to each of the sources appears in the numera- 
tor of equation (B13) in a linear combination. Therefore, if the shield has a uniform tem- 
perature, this temperature can be found by summing up the contributions due to each 
source taken separately. 

When the equations are linear and there is only one source with a finite temperature, 
the shield temperatures are porportional to this temperature. The heat-transfer rates 
are proportional to this temperature raised to the fourth power. In this way, if the source 
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temperature is doubled, the shield temperature will also double while the heat-transfer 
rate will increase by a factor of 16. The additive procedure, when the case resolved it- 
self into a linear problem for a single shield between two sources, is illustrated as fol- 
lows: 


T = t 2 T = t 2 


4 4 4 4 

t: = T* T* = 

s si s s2 

. T s = + T s2 

Q = Q Q = Q 

Q = Q + Q 

T = T. T = 0 

T = T 1 


The procedure can be extended to the situation in which more than one shield is present. 

iu 

The heat-transfer rate into the j tn element is the sum of conduction and radiation 
heat-transfer rates and equals zero 


q c,j +q r 


+ q . = 0 

J+n 4 r,j 


Substituting values from equations (Bl) and (B6) gives: 


M - T .i) , ( T i-1 - T t) jS 

*j+l,j P=1 

So long as the ratio of the radii (r./rj where i = j ± 1) stays constant, the thermal 
resistances are not a function of the radiating area. Each term on the right side of equa- 
tion (B15) is proportional to the radiating area. However, for a linearly scaled system in 
which the spacing ratio (L/R) between surfaces remains constant, a reference area Aj 

may be factored out. Since it is assumed that the projection of all surfaces is circular, 

2 

this reference area is proportional to its lateral radius squared (A^ aeR^). When all prop- 
erties are independent of temperature, it is convenient to express temperatures as a ratio 
to a reference value T 2 . Physically, Aj could be the area of a cold source, while T 2 
could be the temperature of a warm source. Equation (B15) can be nondimensionalized to 
give: 


(*)p,n +j ( E ) P + (*%>%] < B15 > 


(T 
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J 




K is a nondimensional constant determined by the shape of the reference area A^. For 
a plane source 
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APPENDIX C 


HEAT TRANSFER BETWEEN SPECULARLY REFLECTING SURFACES 

The purpose of this section of the report is to present the relations used to determine 
the shield temperature and heat-transfer rates for diffusely emitting, specularly reflect- 
ing surfaces of uniform radiosity. Because of the assumptions of uniform radiosity and 
uniform emittance, there is no thermal gradient in the shields. 

The equations for the heat -transfer rate between two plane specular surfaces having 
uniform radiosity (diffuse emissivity) are given in reference 11. The sources appear 
mirror -like so that the view factor between surfaces is a function of the number of re- 
flections encountered by the bundle of energy. Figure 32 gives an illustration of the way 
the apparent spacing increases due to each reflection for the specular surfaces. The sur- 
faces numbered 1 and 3 are specularly reflecting, while the other two surfaces are simply 
windows to let the radiant energy escape. The apparent spacing for the direct view factor 
from 1 to 3 is the actual spacing between 1 and 3. The apparent spacing for the view fac- 
tor from 1 to 1 after being reflected off 3 is twice the actual spacing. Thus, the energy 
appears to be coming from the dashed surface immediately to the left of surface 3. Fig- 
ure 33 is a schematic of a shield and two sources. The surfaces are specular so that the 
angle of incidence equals the angle of reflection (/3 = j3’) . 

The net heat transfer away from a surface is : 

Q = ae T 4 - aH (Cl) 

A 


where H, the incident energy per unit time and area, is composed of two parts. The 
first part is the energy from the other surface, and the second is the energy from the 
surface under consideration which is incident on itself after it has been reflected off the 
other surface. For surface 1: 


H 1 " ae 3 T 3 |f 1, 3 + F l, 3(1, 3) p l p 3 + F l, 3(1, 3, 1, 3) p l p 3 + •••]+ ae i T i [fl, l(3) p ! 

2 3 2 I 

+ F l, 1(3,1, 3) P 3 P 1 + F l, 1(3, 1,3,1, 3) P 3 P 1 + • ' J 


(C2) 


The view factor from surface 1 to surface 3 is used so as to have only one area (A^) 
in the equation. This is done by the use of the reciprocity theorem which states that: 
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F 1,3 A 1 = F 3,1 A 3 



Rewriting equation (C2) in terms of an infinite series gives: 


H 1 ~ CTe 3 T 3 


E 

?7=o 


VnV 


1 , 3(1^, 3^)^1^3 + a 


’■A Y 

77=1 


1, K 3 7 ?, I 7 !- Y&i 


■M " 1 


(C3) 


The view factor between two disks is given in reference 1 and the view factor is a 
function of the spacing ratio F = g(L/R). As the apparent distance between the shields 
increases due to an increasing number of reflections, the view factors decrease due to the 
widening distances between images. Returning to figure 32 shows that the direct view 
factor between a pair of surfaces (F^ g) is a function of the actual spacing between sur- 
faces Fj g = g(L/R) . The view factor from a surface to itself as seen in the other sur- 
face (Fj 1 ( 3 )) i s a function of twice the actual spacing between the surfaces F^ 

= g(2L/R). The view factor after two reflections from the opposite surface is a function 
of three times the actual spacing F^ g^ g^ = g(3L/R). In general 


F l,3(l ri ,2 r ) = s((2V + 1)L/R) 


and 


F 1 1 ( 3 ^ l 77-1 ) = S(2 V L/R) for 77 > 0 (C4) 

Equation (C3) is evaluated by taking a large number of terms in each of the series and es- 
timating the remainder. The terms in the remainder are estimated by assuming the re- 
mainder forms a geometric series. For a constant view factor, the remainder would be 
exactly a geometric series. The value of the parameter in the geometric series is the 
ratio of the two final terms. Equation (C3) is approximated as: 


H 1 “ CTe 3 T 3 


F i , s ( i 5 + 1 , s 5 +1 ) 


+ oe i T i 


E ( F l,l(J , ,l , ‘VK ') 

> 1=1 


~l\ F l,l(lW,^ 


, F i . 1(1^2. 

F l.l(l W ,l ! ) 


(C5) 
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The value of ? was chosen so that the remainders (the last term in each bracket) 
were less than 0. 1 percent of the first term in the bracket which is the value found by 
summing over £. 

A shield placed between two sources has no net heat transfer to it in steady-state con- 
ditions. This fact can be utilized to determine the shield temperature. Figure 33 gives 
an illustration of a single shield placed between two plane sources. The four-surface sys- 
tem is divided into two two -surface systems, and the heat-transfer rate to the surfaces of 
the shield can be found in terms of the source temperature for each system. The sum of 
the heat transfer to the shield from both systems is zero and the shield temperature is 
determined. 
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APPENDIX D 


DERIVATION OF THE DIFFERENTIAL EQUATION FOR 
A STRUCTURAL MEMBER 

The equation giving the temperature distribution along the strut is derived by con- 
sidering the heat balance on an element of the structural member. Figure 34 gives a 
schematic of a strut. Heat enters the element by conduction from the warmer adjacent 
element and leaves by means of conduction to the colder adjacent element. Heat is trans- 
mitted to or from the element by radiation on both the inside and outside surfaces of the 
element. 

The rate at which heat enters the element due to conduction at the position x on the 
strut is given by: 


^in, c 



kir / ,2 ,2\ dT 

4 \° V dx 


(Dl) 


The rate at which heat leaves the element due to conduction at x + Ax is given by: 


%ut, c 


-ka— i 
dx 


x+Ax 


k7 r (,2 

~ 4 \° 



dT 

dx 


'x+Ax 


(D2) 


Since the environment is assumed black (e = 1), energy which is emitted or re- 

b 1 

fleeted from the outside surface of the strut is absorbed by the surroundings and cannot 
return to any part of the strut. The rate at which heat is absorbed by the element due to 
the surroundings is given by: 


x m 


— (T e A. f 
sr sr sr sr 


4 4 

T a = crA f T^ 1 oi 
-x sr o sr sr-x sr o 


The reciprocity theorm for view factors yields the relation: 


A sr f sr-x = A x f x-sr = " d o Ax f x-sr 


The view factor from the element to the surroundings (f ) is unity. 

X bx 

The rate at which heat is radiated from the outside surface toward the surroundings 

4 

is crire 0 d 0 Ax T . It is assumed that the surfaces are gray; therefore, ^ 0 = a 0 - The 
net rate of radiant heat transfer from the outside surface of the element is therefore: 


37 


III 



%ut,r = ™Vo ax ( T4 - T sr) 


(D3) 


The net rate at which heat is absorbed on the internal surface due to radiation from 
the rest of the strut is found by Gebhart’s solution for problems of radiant interchange in 
an enclosure. This method has been described in appendix B. 

In attempting to determine the temperature distribution along the strut, a differential 
equation is derived. In deriving the equation, the limit is taken as the area for radiant 
interchange approaches zero. The view factor to an infinitesimal area from another sur- 
face approaches zero, though the view factor from the infinitesimal area to the surface is 
not necessarily zero. It is desirable, therefore, to express view factors as going from 
the infinitesimal element to the other surface rather than going from the other surface to 
the infinitesimal area. 

The net radiant heat-transfer rate for the internal surface is found by dividing the 
strut into n + 1 sections. One of these sections is the element under consideration and 
the size of this element approaches zero. Physically, this represents the division of the 
strut into n sections. The amount of radiant energy which leaves the element and is 
later absorbed by the element goes to zero as the size of the element decreases. If j 
represents the element under consideration, this last statement is equivalent to saying 
that Bjj approaches zero in the heat-transfer balances. (The mathematical justification 
for this is shown at a later stage in the development of the equations.) In addition to the 
n physical sections of the strut, there are two additional sections necessary to complete 
the enclosure. These two sections are the end disks of the strut. 

The set of algebraic equations for the n + 2 nonzero values of B.. for the heat- 

th ” 

transfer rate to the j element are: 


B lj = e j f lj + ^l f ll B lj + ^2 f l2 B 2j + ' • • + ^k f lk B kj 
B 2j = e j f 2j + ^l £ 2l B lj + P2 f 22 B 2j + • ■ ■ + ^k f 2k B kj 


(D4) 


B kj = £ j f kj +P l f kl B lj + ^2 f k2 B 2j + * • • + P k f kk B kj 


where k = n + 2. 

These equations can be rewritten as: 
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( p l f ll " ^lj +p 2 f 12 B 2j + • • • +p k f lk B kj - " e j f lj 
p l f 21 B lj + ( p 2 f 22 ” X ) B 2j + • • • +p k f 2k B kj = " e j f 2j 


p l f kl B lj +p 2 f k2 B 2j + • ' • + ( p k f kk - X ) B kj = - e j f kj 


The reciprocity theorem for view factors gives: 


A.f.. = A.f.. 
i ij 3 ]i 


If each side of every equation is divided by A. and the view factors on the right side are 

th J 

transformed so that they are from the j element to another source, one obtains: 
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a j 
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f'i‘ki -r 1 

j j 


A. 

] 


p l f 21 + ( p 2 f 22 ” + • • • + p k f 2k " " 4^ 
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A i 


_fjjk 


(D5) 


These equations can be placed in matrix form 


Cu = v 


(D6) 


where C is a coefficient matrix and is independent of j. The entries of C are: 


(C) pi = P i f pi' 6 pi 
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v is a column vector with the entries (v) = -(e.f. )/A„. u is a column vector which 

_ p j jp P 

contains the values of B^/Aj . u is found from equation (D6) by taking the inverse of C 

u = C -1 v (D7) 

th 

The net heat-transfer rate to the j element due to internal radiation is: 


q 


in, r 


i=l 


B ij®i 


- B..) 
IV 




(D8) 


It should be noted that E contains only entries for the finite area elements, and there- 
fore the emissive power of the element under consideration (j) is not an entry in E. 
Equation (D8) can be expressed in terms of the column vector u to give: 


q in, r A j 


K 

2 < 5 v E >i 


Ll=l 


- (1 - 


“ B ji 


had been included as the k + 1 equation of (D4), the above equation would be: 


q in, r 


- a j 


k 


i=l 


(u) i (E) i - aT?e j 


” T iVi (u| W 


(D9) 


Aj is the internal surface area of the element under consideration (ird^ Ax). In sub- 
sequently forming the differential equation, the right side of equation (D9) is divided by 
Ax and the limit is taken as Ax approaches zero. When this is done, the last term in 
equation (D9) goes to zero. Therefore, the net radiant heat transfer into the internal sur- 
face of the element is: 


q. = 77 d Ax 
in, r l 


K 

X) (H) l (E) i - olfe, 
i=l 


(DIO) 


Equating the rate at which heat enters the element with the rate at which it leaves 
gives: 


q in, c q out, c q out, r q in, r 
Expanding each term and collecting common factors gives: 


(Dll) 
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(d 2 - d 2 ) 

/dT 

dT 

\ a o a J 

ydx 

x _dx 


y ” a x (ir) i A i £ i T i 

x+Ax/ i=l 


Taking the limit as Ax — 0 for both sides of the equation yields: 

k ( d ° ' ^ — = T “( £ A + e o d o) - £ o d o T sr - d t [x (5) i A i 6 i T i 


40 dx 2 


i=l 


(D12) 


In equation (D12), the subscript j has been dropped. However, the entries in u do de- 
pend on the relative position of the differential element and the finite area sections. 

The differential equation is integrated from zero to l where l is the length of the 
strut. By letting £ = x/Z, the range of the independent variable is from zero to one. 
Since dx = Zd£, equation (D12) becomes: 


€ - d ?) 


d 2 T 


4crZ 


dr 


- t4 ( £ A + £ o d o) - £ o d o T sr - d t 




(D13) 


Li=l 


When the thermal conductivity and emissivity are independent of temperature, non- 
dimensional temperature profiles can be obtained. Since both end temperatures 
(T 1 and T 0 ) are known, the differential equation has known boundary values. Dividing 

1 * 4 

each side of the equation by Tg yields: 



Letting 9 = T/T 2 one obtains: 
^ - d ?) 


4aZ 2 T| 


d 2 * 

d| 2 


= 0 VVsr“ d t 


Li=l 


(D14) 


The boundary conditions are 9^ - 1 and 9 ^ = T^/Tg. Since the values of 6>. are func- 
tions of the temperature distribution along the strut, they are determined by 9. 

When the strut is a thin-walled tube, the cross-sectional area is approximately ?rdt 
where t is the thickness of the tube and d represents either the inside or outside diam- 
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eter. For a thin -walled tube equation (D14) becomes: 


kt 





£ o e Jr 


22 (u)^^ 

i=l 


(D15) 


When the temperature of the surroundings goes to zero and the internal emissivity is 
reduced to zero, equation (D15) becomes: 



crZ 2 Tp£ 0 4 
2 o 

kt 


(D16) 


The total heat-transfer rate out of the cold end of the strut is the sum of the conducted 
heat transfer out and the radiant energy absorbed by the end disk. For illustration, con- 
sider the end at which x equals zero. From figure 34 it is seen that the end disk is the 
k^ element in the system. The temperature of the strut at this end is Tj and for 
closed disks T k would probably equal T^. The total heat-transfer rate is: 



+ Q = ka — 
r dx 


k 

+ 

x=0 1=1 


K 

X K - 5 ikH 


(D17) 


Note that since the end disk is a finite area element, Bj^ is not zero. In reference 10, 
it is shown that for gray surfaces. Equation (D17) can be rewritten 

as: 






Introducing nondimensional parameters gives: 


4Qn 


4 2 

« T i d r 


. k ( d o - 0 
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+ e. 


£=0 


X] Ki - s ik>f 


(D18) 


For a constant internal emissivity e k = The absorption factors B k > depend only on 
the internal emissivities and 1 / d t of the tube. 
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APPENDIX E 


SOLUTION OF EQUATIONS 

The purpose of this appendix is to discuss the solution of the equations presented in 
appendix B for the shadow-shield system and the equations presented in appendix D for the 
thermal analysis of a structural member . - 

Shadow-Shield Equations 

In appendix B, equations are given for the shield temperature and concomitant heat- 
transfer rates when the equations are linear in T^. As the conducted heat transfer in- 
creases relative to the radiant heat transfer, it becomes more advantageous to solve the 
equations for the shield temperatures assuming that the equations are linear in T. The 
various terms contributing to the heat balance for an element are shown in figure 30. The 
conduction heat-transfer rate into the element is given by equation (B6) . The equa- 
tions of (B8) give the radiant heat-transfer rate into each element of the shield. In steady 
state, the net heat-transfer rate into each element is zero. Adding equations (B6) and (B8) 

J.U 

for the j element gives: 



Expanding this equation and placing terms involving known temperature sources on the 

J_ V- 

right side yields for the j element: 



Factoring out the temperature of each element of the shield gives: 
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CrA/lf 


T. 

1 


H (*)p, 

LP=1 


n>)p + (*W4 + 21 (Op.iM 


p=n+l 


An equation can be written for each of the n elements of the shield. This results in 
the matrix equation: 


Y0 = P 


(E2) 


Y is an n-by-n matrix containing the entries: 


( Y ). . = (*) . . € . + (**). . 6 ° 
V /!,] \ Vn+j , n+i fc n +3 V /],i ] 


lA ] T r\/ 1 


5 i,j+l + 5 i,j-l 


5-1,3 ^j+lj/ ^i-l,i ^i+l,i 


0 is a column vector containing the temperature for each annulus 


4 - t j 


P is a column vector containing the negative of the radiant heat transfer to each annulus 
from the sources with known temperatures. 


( p )t - - 


2 (*U E ) P *(i.i)c 2 Mp.iMp 

p=l p=n+l 


The temperatures are found by taking the inverse of Y in equation (E2). This gives: 


0 = Y 


(E3) 


The equation for the matrix Y can be expanded to give: 
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When one of the sources is solar radiation, the matrix Y and the vector P become: 



h = - iz (*%>% 

p=n+l 


From the equations presented above and those of appendix B, it can be seen that when 
the equations are linear and only one shield present, the shield temperatures are deter- 
mined explicitly. The equations can be nonlinear due to a finite shield conductivity or 
temperature -dependent surface properties. When the equations are nonlinear or there is 
more than a single shield, an iterative procedure is used. This iterative procedure sim- 
ply sets a new trial value for each annulus of each shield at a value lying between the old 
estimate and the latest calculated value. When the shields are conducting, the new esti- 
mate lies close to the old estimate. The computer program chooses this degree of close- 
ness so as to minimize computer time, but also insuring that the system of equations re- 
mains stable. The tolerance used to determine when the iteration has converged is 
0.0055 K (0.01° R) . The difference between the calculated and estimated temperatures is 
less than this value for each annulus of each shield at convergence. To insure that the 
tolerance is sufficiently fine, the first case in each run is automatically rerun with the 
tolerance halved. 

When the conduction term does not enter into the equations, convergence is fairly 
rapid. With a two-shield system and no conduction, about five iterations are required to 
obtain a solution. The number of iterations increases with the number of shields involved 
going to about 20 iterations for a four-shield system with no conduction. The number of 
iterations required when the shields have a uniform temperature is about the same as 
when there is no conduction between annuli. If there is conduction between annuli, the 
number of iterations increases drastically so that a factor of 10 in the number of itera- 
tions required would not be unexpected. 

The execution time for the program is proportional to the number of iterations and 
the number of shields involved. It is proportional to the square of the number of annuli 
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into which each shield is divided. A two-shield configuration with conduction and 14 annuli 
per surface would take about 2 minutes to execute on an IBM 7094 mod It digital computer. 
If conduction were neglected, the execution time would be less than half a minute. 

In order to obtain heat-transfer rates for surfaces of nonuniform radiosity, each of 
the annuli should represent a differential area. To closely approximate this directly using 
a computer would result in an exceedingly large amount of running time. The results for 
nonuniform radiosity surfaces, which constitute the bulk of the data in this report, are ob- 
tained by graphing the heat -transfer rate as a function of the reciprocal of the number of 
annuli. This procedure is illustrated in figure 35. In this figure, the ratio of the heat- 
transfer rate to the heat-transfer rate for uniform radiosity surfaces (a single annulus per 
surface) is given as a function of the reciprocal of the number of annuli per surface. The 
ratio refers to the heat-transfer rate to the colder source, and is given for two emissivi- 
ties and two spacing ratios. The single shield was spaced evenly between the two sources 
and all surfaces had the same emissivity. The point on the graph corresponding to an in- 
finite number of annuli (1/n = 0) corresponds to the heat -transfer rate for surfaces of non- 
uniform radiosity. The solid line represents points for which data were actually gener- 
ated, and the dashed line represents an extrapolation. It can be seen from this figure that 
the accuracy of the extrapolation decreases with decreasing emissivity. Also, it can be 
seen from this figure that increasing the number of annuli (this results in increased com- 
puter time) reduces the uncertainty in the extrapolation. 


Thermal Equations for Struts 

The differential equation presented in appendix D for the thermal analysis of a struc- 
tural member is a second-order nonlinear equation with known boundary values. A finite 
difference scheme is used to solve the differential equation. In this technique, the length 
of the strut is taken as unity. The strut is divided into several equally spaced intervals, 
each with an axial length h. The second derivative expressed in terms of the dependent 
variable and the step size is 


d T i T i-1 - 2T i + T i + 1 


dr 


(E4) 


Here i represents the particular point along the strut for which the differential equa- 
tion is being approximated. Because the range of the independent variable is from zero to 
one, h = l/(m + 1) . Where m is the number of interior points on the strut at which the 
differential equation is evaluated. The first point along the strut at which the differential 
equation is evaluated is a distance h from the end. Since this is a boundary value prob- 
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lem, the end temperatures (T and T m are known a priori. There are m equations 
of the same form as equation (E4) and these can be written in matrix form as: 


O = X© 


(E5) 


n is a column vector of size m containing the evaluations of the second derivatives. 


(n) p = n 


d 2 T 


d^ 


> P " ( 6 l,P T o + 6 m,p T m+l) 


(E6) 


X is an m-by-m coefficient matrix containing the coefficients of the temperatures, 
entries of this matrix are clustered along the diagonal and have the values: 

^p,q = " 26 p,q + Vq+i + Vq-i 

0 is a column vector containing the temperatures at each node point 


The 


(E7) 


(©) 


q 



(E8) 


Equation (E5) can be solved to yield the temperature distribution along the strut. This 
yields the equation: 


0 = x *0 (E9) 

Since the second derivative is dependent upon the temperature distribution along the 
strut, an iterative procedure is used to determine the correct temperature distribution. 
This iterative procedure is much like the one employed in the solution of the equations for 
the shadow shield. 
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APPENDIX F 


COMPUTER CODE FOR SHADOW-SHIELD CALCULATIONS 

The purpose of this section of the report is to discuss the computer code used to solve 
the shadow-shield equations. The procedure for reading in the necessary data and the re- 
sulting output are discussed. Following this are general descriptions of the various rou- 
tines used to make up the code. Flowcharts are given for these routines. This appendix 
concludes with a listing of the program and a sample output. 


Description of Input 

All of the necessary data enter the program through READ statements in the routine 
INPUT. For every case in a run this routine is called by SHIELD. Whenever view factor 
data for nonplanar sources are read in INPUT is called by TBVUF. All data enter the 
program through the use of NAMELIST statements. NAMELIST is an input/output feature 
of the FORTRAN IV language in which this code was written. The use of this feature is 
described in FORTRAN IV manuals such as reference 12. 

The NAMELIST input/output feature uses names in place of Format numbers in the 
read and write statements. There are two names used in the reading statements of the 
INPUT routine. These two names are INPUT1 and INPUT2. The first of these is used to 
read in all of the data when there is only planar sources or solar radiation. When there 
is one or more nonplanar sources present data are read in for view factors through the 
use of the name INPUT2. 

The necessary data are described in the listing of the routine INPUT. At the begin- 
ning of a run, prior to reading any data, many input variables are assigned certain values. 
This is done to reduce the amount of data which are read in. NAMELIST provides a flex- 
ible means of entering data into core. It is not necessary to specify all of the variables 
which are associated with a NAMELIST name, therefore, only those data necessary to 
execute a case need be provided. There is no provision for storing data for more than 
one case, so that at least one variable must be read in for each case in a multiple case 
run. However, only the data which change need be specified for each case. The program 
continues to search for new data until all of the input has been processed. (This program- 
ming results in a warning diagnostic for the routine SHIELD at compilation.) Because of 
the flexibility of NAMELIST, it is possible to inadvertently omit some of the data neces- 
sary to run a case. This will result in either error messages during execution or erro- 
neous results. The latter situation is the more troublesome of the two. In order to be 
able to check the data, a listing of the input is provided at the beginning of each case. 
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This listing may be more extensive than the input because it consists of all of the data 
which have a bearing on the results. The sample case given at the end of the listing is for 
two nonconducting shields and a planar tank. The outermost shield is exposed to solar ra- 
diation. The variables required are determined from the listing of the routine INPUT. 

The values used for each variable are found in the output. 


Description of Output 

The output from a successfully executed case is written by the routine SHIELD. A de- 
scription of the variables printed out can be found by consulting the listing of this routine. 
In addition to the expected output there may be additional output. Faulty data may result 
in an illegal operation such as division by zero. The execution monitor will then print out 
its own message. Faulty data may also result in the routine FACTOR encountering a sin- 
gular matrix. If this occurs, an error message is printed and the run is terminated. 

As mentioned previously, a large number of iterations is required when the shields 
are partially conducting. If the number of iterations exceeds the value of the variable 
ITMAX, a message is printed along with the last estimates of the shield temperatures. 
When the shields are partially conducting, a solution may become divergent from a given 
starting condition, or may become divergent due to the means used to calculate the solu- 
tion. There are two starting conditions for partially conducting shields and these are 
specificed by the variable IVCOND. These starting conditions are either uniform temper- 
ature shields or nonconducting shields. The solution for partially conducting shields can 
be found by assuming that the equations are linear in T or linear in T . The initial as- 
sumption concerning the linearity of the equations is governed by the value of the variable 
NCALCR. If, starting with the specified values for IVCOND and NCALCR, the solution 
appears to be diverging, one of the variables is changed, and the case is restarted. If, 
after the four possible combinations have been tried, the solution still appears to be di- 
verging, the case is abandoned. Each time a change is made a message is printed. The 
current values of IVCOND and NCALCR are stored in ICOND and NCALC, respectively, 
and these values are printed as part of the message. All of these messages occur from 
the routine SHIELD. 

Whenever a case is run involving one or two nonplanar surfaces, a check is made to 
insure that the necessary view factor data have been read in. This check is for the spac- 
ing ratio between surfaces, the type of surfaces, the ratio of the axial radius to the later- 
al radius for the nonplanar surfaces, and the number of annuli used. If the check fails, 
an error message is printed from the routine TBVUF and the run is terminated. 
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Description of Computer Code 


This section of the appendix discusses the computer program which performs the cal- 
culations of the shadow-shield equations. This program is composed of several routines 
each of which performs a specific function. In the program, one of the sources is re- 
ferred to as a heater and the other source is referred to as a tank. It is convenient to 
describe the computer program using these terms for the sources. The following is a 
brief description of each of the routines used in the program. Most of the routines have 
an accompanying flowchart and following some of the routines is a list of the FORTRAN 
symbols which correspond to symbols used in appendix B. The FORTRAN variables which 
are needed as input are described in the listing of the program and are not repeated here. 
Some variables are used throughout the program while others are limited to one or two 
routines. The FORTRAN symbols which are used throughout the program are given after 
the routine SHIELD and the others are given after the first routine in which they appear. 
The list of FORTRAN symbols is not exhaustive but is limited to the major variables. 
Accompanying the description of the FORTRAN variables is a key. This key describes 
the subscripting of variables and is explained after the discussion of the routines FACTOR 
and INVERT. 

INPUT . - This subroutine reads in the necessary data. All of the input to the pro- 
gram enters through this routine. Figure 36 is a flowchart of this routine. The listing of 
the routine contains a description of the input. 

SHIELD. - This is the main program. One of the functions of this routine is to pro- 
vide a means of controlling the calling sequence to the subroutines which perform the ac- 
tual calculations. Also, the iterative procedure for calculating the shield temperatures 
is done by this routine. With the exception of error messages, all output is generated in 
this routine. Figure 37 is a flowchart of this program. The numbers in connectors or to 
the left of a function block represent statement numbers in the program. 


FORTRAN 

name 

Variable in 
appendix B 

Description 

Key 

AREA 

A 

The area of each annulus of each surface 

1 

AB 

a 

The absorptivity of each annulus of each 
surface 

2 

EMS 

€ 

The emissivity of each annulus of each 
surface 

2 

REF 

P 

The reflectivity of each annulus of each 
surface 

2 

COND 

k 

The thermal conductivity between annuli 
of a shield 

1 
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FORTRAN 

name 


Variable in 
appendix B 


Description 


Key 


RES 

T 

FNNHTT 


FNNTTH 


dt The thermal resistance between annuli of 1 

a shield 

T The temperature of each annulus of each 1 

surface 

f The view factor from the surface of an 3 

annulus facing towards the tank to an 
annulus on a surface facing towards the 
heater. If there are no shields present, 
this is the view factor from an annulus 
on the heater to an annulus on the tank 

f The complement of FNNHTT except that this 3 

is the view factor from an annulus towards 
the tank to an annulus towards the heater 


TBVUF . - This routine calls INPUT to read view -factor data when nonplanar sources 
are present. These view factors were obtained using the computer program given in ref- 
erence 13. This routine is set up so that all the data which are generated for nonplanar 
surfaces may be kept together and read in as a unit. Figure 38 is a flowchart of this sub- 
routine . 

VUFAC1 . - This routine calculates the radii and areas for planar surfaces and sets 
up the call to the routines TBVUF and SUBVF in order to generate the view factors be- 
tween annuli. VUFAC1 calculates the view factors between whole surfaces by summing up 
view factors. Figure 39 is a flowchart of this routine. 

SUBVF. - This subroutine calculates the view factors between planar surfaces. The 
view-factor relations that are given in reference 1 are used in this routine. Figure 40 is 
a flowchart of this routine. 

SCHLC . - This routine calculates the temperatures of the shields assuming that the 
temperatures of adjacent surfaces are known for each shield. Figure 41 is a flowchart of 
this routine. 


FORTRAN 

name 

PHI 

Q 


Variable in 
appendix B 

<P 

Q 


D 


D 


Description 


The solar flux 

The rate of heat transfer to an annulus 
from known temperature sources 

The matrix containing the absorption fac- 
tors times emissivity 
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COEF . - This’ routine determines the absorption factors for each system consisting 
of two surfaces and the surroundings. Figure 42 is a flowchart of this subroutine. 


FORTRAN 

name 

Variable in 
appendix B 

Description 

Key 

A 

C 

The coefficient matrix in the solution of 
the absorption factors 

4 

Y 

V 

The matrix containing the negative of the 
view factors times absorptivity 

4 

AI 

C' 1 

The inverse of A 

4 

B 

u 

The matrix containing the absorption 
factors 

4 


QDOT . - This routine calculates the heat-transfer rates to the tank as well as both 
the adjacent surface and the surroundings between these two surfaces. Figure 43 is a 
flowchart of this routine. 

FACTOR and INVERT . - These two routines are used to invert a matrix. The meth- 
od used is given in reference 14. As a consequence of the way these two routines are 
programmed, four warning diagnostics result at compilation. 

Explanation of keys . - 

Key 

1 These variables are doubly subscripted. The first subscript refers to an annulus 
of a surface. The first element is the innermost annulus. The second subscript 
refers to the particular surface beginning with the heater. Thus AREA (5, 2) is 
the area of the fifth annulus of the first shield. 

2 These variables are triply subscripted. The first two subscripts have the same 
meaning as the two subscripts for variables with a key 1 . The third subscript 
denotes the particular side of the shields. The side designated a 1 faces the 
heater while side 2 faces the tank. Thus AB (5, 2, 2) is the absorptivity of the 
fifth annulus of the first shield on the side feeing towards the tank. 

3 These variables are triply subscripted. The first two subscripts refer to annuli 
on the originating and receiving surfaces, respectively. The third subscript re- 
fers to the pair of surfaces involved. These subscripts are arranged in the same 
fashion as the two previous keys. Thus, FNNTTH (1, 10, 1) refers to the view 
factor from the innermost annulus to the tenth annulus and is from an annulus of 
a shield or source adjacent to the heater to an annulus on the heater itself. 

4 These variables are doubly subscripted and represent matrices. The two sub- 
scripts refer to the rows and columns of each matrix, respectively. 
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$ I BF rc INPUT 


DECK 


SUBRCUTINE INPUT(JR) 

THE PURPOSE OF THIS ROUTINE IS TO READ IN ALL OF THE DATA TO THE 
PROGRAM 

EQUIVALENCE 

I ( NE ( 1 ) * NE I B ) , (NEI21) » NE IT ) , ( NE ( A 1 ) » NE2B ) , (NEI61) ,NE2T), 

2 C NE { 8 1 ) , NE3B ) , ( NE ( 10 1.) , NE 3T ) , ( NE ( 12 I ) f NE4B ) , { NE ( 141 ) , NEAT ) , 

3 ( NE ( 161 ) ,NE5B ) , (NE( 181 ),NE5T) , {.MAC 1 ) ,NA1B) , (NA(21> ,NA1T), 

4 ( N A ( A 1 ) ,NA2B) , (NA161) ,NA2T) , (NA181) ,NA3B) , { NA ( 1 0 1 ) , NA3T ) , 

5 ( N A ( 121 ) ,NA4B ) , ( NA ( 14 1 ) , NA4T ) , ( NA ( 16 1 ) , NA 5B ) , ( NA ( 18 1 ) , NA5T ) 

EQUIVALENCE 

KCVABS(l) » T B ABS 1 ) , (CVABSI5) ,TBABS2> , (C VAB S ( 9 ) ,TBABS3) , 

2 ( CVABS ( 13) » TBABS4 ) , ( CVABS ( 17 ), TBAB S5 ) , (CVABS (21 ) ,TBABS6 ) , 

3 (CVABS( 25 ) ,TBABS 7 ) , ( CVABS( 29 ) , TBABS8 ) , ( C VAB S ( 33 ) , TBABS9 ) , 

4(CVABS( 37) »TBABSL) , (CVEMS(l) » TBEMS1 ) , ( CVEMS ( 5 ) , TBEMS2 ) , 

5 ( CVEMS ( 9 ) ,TBEMS3) , ( CVEMS ( 13 ), TBEMS4 > , ( C VEM S ( 1 7 ) , TBEMS5 ) , 

6(CVEMS(21) » T B EM S 6 ) * < CVEMS ( 25 ), TB EMS7 ) , (CVEMS (29 ) ,TBEMS8 > , 

7(CVEMS( 33) .TBEMS9) , ( CVEMS ( 37 ), TBEMSL ) 

EQUIVALENCE 

l(CVCONDd) , T BCON 1 ) , ( CVCONDI 5 ) , TBC0N2 ) , ( C VCOND ( 9 ) , TBC0N3 ) , 

2(CVC0ND( 1 3 ) , T BC0N4 ) , ( CVCONDI 1 7 ) , TBC0N5 ) 

DATA ( ID(I), 1=1,10) /IH1, 1H2, 1H3, 1H4, 1H5, 1H6, 1H7, 1H8, 1H9, IHL / 
DATA ( DNE( I ) , 1=1, 12 ) /72H NEH = N E 1 B =NE 1 T =NE2B =NE2 T =NE3B =NE 3T = 
1NE4B = N E 4 T =NE5B = N E 5 T = NET = / 

DATAIDNAt I ) , 1 = 1, 12) /72H NAH = N A 1 B =NA l T =NA2B =NA2 T =NA3B =NA3T = 
1NA4B =N A4T = N A 5 B = N A 5 T = NAT = / 

DATA ( 0 T E M ( I ) , I = l , 2 ) /6H TH =,6H TT = / 

DATA DAHS , DATS / 6H AHS = , 6H ATS = / 

COMMON C( 1000) 

EQUIVALENCE 


1 ( C ( 1 ) 

, ABSR 

) 

f 

(C( 11 ) 

, ARES 

) 

» 

( C { 2 1 ) 

♦ 0 

) 


2 { C ( 3 1 ) 

tES 

) 

t 

( C ( 41 ) 

, IVCOND) 

» 

(C (42) 

, ICONOR) 


3 ( C ( A 3 ) 

? UMAX 

) 

» 

(C( 51 ) 

, NAH 

) 

f 

(C ( 71) 

9 NAT 

) 


4(C (91 ) 

,NA 

) 

» 

( C ( 29 1 ) 

, NEH 

) 

» 

(C ( 311 ) 

, NET 

) 


5 ( C ( 3 3 1 ) 

» NE 

) 

f 

( C ( 5 3 l ) 

, NCOND 

) 

» 

(C( 541 ) 

fNCASE 

) 


6(C(542) 

,NEL IPS) 

t 

(C( 543) 

, NR ADS 

) 

t 

(C ( 544) 

t NRNG S 

) 


7(C(545) 

, NS HL DS ) 

t 

{ C ( 546) 

,NCALCR) 

f 

< C C 54 7 ) 

9 N V U F H 

) 


8 (C ( 548 ) 

t N VUFT 

) 

f 

( C ( 549 ) 

, PHI 

) 

» 

(C ( 550 ) 

9 PI 

) 


9 ( C ( 5 5 L ) 

,RS 

) 

f 

( C ( 56 1 ) 

t RS H 

) 

9 

(C ( 562 ) 

9 R ST 

) 


A { C ( 56 3 ) 

, RH 

) 

t 

( C ( 564 ) 

»RT 

) 

» 

(C ( 565 ) 

f SOLAR 

) 


B ( C ( 5 66 ) 

, SIGMA 

) 

t 

(C( 567) 

t TOL 

) 

» 

(C (568 ) 

tWI 

) 


C ( C { 5 7 1 ) 

, SOLARS) 

* 

(C< 591 ) 

f STHICK ) 

f 

( C ( 60 1 ) 

» C V A B S 

) 


0(C (641) 

, CVCOND) 

» 

( C ( 68 1 ) 

, CVEMS 

) 

» 

(C( 721 ) 

9 TH 

) 


E ( C ( /4 1 ) 

.TS 

) 


( C W 51 ) 

» TT 

> 

9 

(C ( 771 ) 

9 REF S 

) 



COMMON / ELIP2 / FN , FL , RMHD , RMTD , XLOR , MXS , ATS(IO) , 
1 AHS ( 10 ) , F(IOO) , LIST , NR 1 


53 



onnoonoonooonono 


INPUT 


EFN SOURCE STATEMENT 


IFN(S) 


DIMENSION ONE ( 12 ) , DNAI12) , DTEMI2) , NF(IO) ,ID(lO) 
DIMENSION CVABS (4,10) , CVEMS(4,10) , CVC0ND(4,5) , NEHI20) , 


1 


NAHI20) , NET ( 20 ) , NATI20) , NEI20 

.10) . NA (20.10) » 

2 


T T ( 20 ) , THC20) , ABSR(IO) » ES(IO) 

, ARES 110) , TSI10I, 

3 


DUO) 





DIMENSION 

TBEMSK4) , 

TBEMS2I 4 ) , 

TBEMS3 (4) 

, TBEMS4 (4) , 

I 


TBEMS5 ( 4 ) , 

TBEMS6C 4 ) , 

TBEMS7 (4) 

. TBEMS8 (4) . 

2 


TBEMS9I4) , 

TBEMSL (4) , 

TBABSI (4) 

, TBABS2 (4 ) . 

3 


TBABS3 ( 4 ) , 

TBABS4I4) , 

TBABS5I4) 

, TBABS6I4) , 

4 


TBABS714) , 

TBABS8C4) , 

TB AB S9 ( 4 ) 

, TBABSL (4 ) , 

5 


TBCON 1(4) , 

TBC0N2 ( 4 ) , 

TBC0N3 (4) 

, TBCCN414) , 

6 


T BCON 5 ( 4 ) 






DIMENSION 

NE1BI20) , NE1T120) , NE2BI20) , NE2TI20) , 

NE3BC20) , 

1 


NE3TI20) , NE4B ( 20 ) , NE4TI20) , NE5BI20) , 

NE5T ( 20 ) , 

2 


NA1BI20) , NA1TI20) , NA2BI20) , NA2T120) , 

NA3BI20) , 

3 


N A3T ( 20 ) , NA4B ( 20 ) , NA4TI20) , NA5BI20) , 

NA5TI20) , 

4 


RS ( 5 ) , STHICK ( 5 ) , NC0NDI5) 




NAMELIST / 

INPUTl / 





I 

ABSft 

# ARES , 

D , 

ES , 

ICONDR 

t 

2 

ITMAX 

, IVCOND , 

NAH f 

NAT , 

NA IB 

9 

3 

NA2B 

, NA3B , 

NA4B , 

NA5B , 

NA1T 

t 

4 

NA2T 

, NA3T , 

NA4T , 

NA5T , 

NCALCR 

9 

5 

NEH 

t NET , 

NELB , 

NE2B , 

NE3B 

9 

6 

NE4B 

, NE5B , 

NE1T f 

NE2T , 

NE3T 

9 

7 

NE4T 

t NE5T , 

NCOND f 

NCASE , 

NELIPS 

f 

8 

NR ADS 

, NRNGS , 

NSHLDS , 

NVUFH f 

NVUFT 

f 

9 

PHI 

t RH , 

RS , 

RSH » 

RST 

9 

A 

rt 

* SIGMA t 

SOLABS , 

SOLAR , 

STHICK 

9 

8 

TBABSI 

, TBABS2 , 

TBABS2 , 

TBABS3 , 

TBABS4 

9 

C 

TBABS5 

, TBABS6 t 

TBABS7 , 

TBABS8 , 

TBABS9 

9 

D 

T8ABSL 

, TBCON 1 , 

TBC0N2 , 

TBC0N3 , 

TBC0N4 

9 

E 

TBC0N5 

, TBEMS1 t 

TBEMS2 , 

TBEMS3 , 

TBEMS4 

9 

F 

TBEMS5 

, TBEMS6 , 

TBEMS7 , 

TBEMS8 , 

TBEMS9 

9 

G 

TBEMSL 

» T H » 

TOL , 

TS , 

TT 

9 

NAMELIST/OUTPUI/ 





I 

ICONDR 

, ITMAX , 

IVCOND , 

NCALCR , 

NCASE 

9 

2 

NELIPS 

, NR ADS , 

NRNGS f 

NSHLDS , 

NVUFH 

9 

3 

NVUFT 

t PHI , 

RH , 

RSH , 

RST 

9 

4 

RT 

, SIGMA , 

SOLABS , 

SOLAR * 

TOL 


NAMELIST / 

INPUT 2 / AHS 

» ATS * F , 

FL , FN , 

LIST , MXS . NR1 , 

2 


RMHD 

, RMTD t XLOR 




NAMELIST / 0UTPU2 / FL , FN , LIST , MXS , NR 1 , RMHD , RMTD ,XLOR 
LIST OF INPUT FOR SHADOW SHIELD PROGRAM 


VARIABLES READ IN THROUGH THE USE OF THE NAMELIST NAME INPUTl 


NAME 

DESIGNATION 

QUANTITY 

CODES 

A8SR 

ABSORPTIVITY OF THE SURROUNDINGS. 
UNLESS INPUT THESE VARIABLES HAVE 
A VALUE OF 1.0 

NSHLDS+ 1 

A 

ARES 

AREA OF THE SURROUNDINGS. UNLESS 
INPUT THESE VARIABLES HAVE A VALUE 
OF 1.0 

NSHLDS+1 

A 

0 

DISTANCE BETWEEN SURFACES. 

NSHLDS+1 

A 

ES 

EMISSIV ITY OF THE SURROUNO I NGS. 

NSHLDS+1 

A 
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I >4 PUT 


EFN 


SOURCE STATEMENT 


IFN(S) 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


UNLESS INPUT THESE VARIABLES HAVE 
A VALUE OF 1.0 

ICCNCR INDICATES THE TYPE OF SHIELDS USED. 1 
0 WHEN THE SHIELDS ARE PARTIALLY 
CONDUCTING. I WHEN THE SHIELDS ARE 
NON— CONCUCT ING. 2 WHEN THE SHIELDS 
HAVE A UNIFORM TEMPERATURE. ICONDR 
IS INITIALLY ASSIGNED A VALUE OF l. 

THIS VARIABLE HAS SIGNIFICANCE ONLY 
IF NRNGS IS GREATER THAN 1, AND 
CANNOT BE ZERO WHEN NRNGS IS 1. 

ITMAX THE MAXIMUM NUMBER OF ITERATIONS 1 

BEFORE THE SOLUTION OF THE 
IS ABANDONED. UNLESS SPECIFIED 
THIS VARIABLE IS ASSIGNED THE 
VALUE 2 CO 

IVCOND 1 WHEN THE SOLUTION FOR PARTIALLY 
CONDUCTING SHIELDS IS STARTED 
FROM A NON CONDUCTING SOLUTION. 

2 IF THE STARTING POINT IS FROM 
THE UNIFORM TEMPERATURE SOLUTION. 


NAH KEY USED TO SPECIFY ABSORPTIVITY NRNGS 

UF ELEMENTS OF HEATER. 

NAT THE EQUIVALENT OF NAH EXCEPT FOR NRNGS 

THE TANK. 

N A 1 6 KEY USED TO SPECIFY ABSORPTIVITY NRNGS 

OF ELEMENTS ON SIDE OF SHIELD 
CLOSEST TO HEATER FACING THE 
HEATER. 

N A I T THE EQUIVALENT OF N A I B EXCEPT NRNGS 

FOR SIDE FACING AWAY FROM HEATER. 

NA2B THE EQUIVALENT OF N A I B OR NAIT NRNGS 


NA3B EXCEPT THAT THESE VARIABLES 

NAA8 REFER TO THE SECOND THROUGH THE 

NA5B FIFTH SHIELDS RESPECTIVELY. 

NA2T 
NA3T 
NA4T 
NA5 T 

NCALCR 0 WHEN THE SOLUTION IS FOUND WTTH I 

THE EQUATIONS LINEAR IN T**4. 

1 WHEN THE EQUATIONS ARE LINEAR IN T. 

NCASE i IF THE LAST CASE IS TO BE RERUN I 

WITH THE TOLERANCE ON THE 
TEMPERATURE DISTRIBUTION HALVED. 

A VALUE OF I IS ASSIGNED FOR THE 
FIRST CASE AND UNLESS A VALUE IS 
INPUT A VALUE OF 2 IS ASSIGNED 
EACH SUCEEDING CASE. 

NCOND KEY DENOTING THE THERMAL CONDUC- NSHLDS 

TIVITY OF EACH OF THE SHIELDS. 

NEH KEY USED TO SPECIFY EMISSIVITY OF NRNGS 

ELEMENTS OF HEATER. 

NET THE EQUIVALENT OF NEH EXCEPT FOR NRNGS 

THE TANK. 

NELIPS THIS VARIABLE HAS MEANING ONLY IF L 


B,C 
B t C 
B.C 

B y C 
B f C 


C 

BfC 

B»C 
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IMPUT 


EFN SOURCE STATEMENT 


IFN(S) 


NE IB 


NE1T 

NE2B 
NE3B 
NE4B 
NE5B 
NE2T 
NE 3 T 
NEAT 
NEST 
NR ACS 


NRNGS 


NSHLCS 

NVUFF 


NVUFT 

PHI 

RH 

RS 

RSH 


RST 

RT 

SIGMA 


VIEW FACTOR DATA FOR A NONPLANAR 
HEATER OR TANK IS READ IN FOR A CASE 
OTHER THAN THE FIRST CASE IN A RUN. 
SET EQUAL TO 0 FOR EACH CASE IN WHICH 
VIEW FACTOR DATA ARE READ IN. CARE 
SHOULD EE TAKEN IF VIEW FACTOR DATA 
ARE READ IN MORE THAN ONCE, SINCE 
SUCEEDING DATA OVERWRITE PRECEEDING 
DATA. 


KEY USED TO SPECIFY FMISSIVITY OF 
ELEMENTS ON SIDE OF SHIELD CLOSEST 
TO HEATER FACING THE HEATER. 

NRNGS 

BtC 

THE EQUIVALENT OF NE I B EXCEPT FOR 
SIDE FACING AWAY FROM HEATER. 

NRNGS 

BfC 

THE EQUIVALENT OF NEIB OR NE I T 
EXCEPT THAT THESE VARIABLES REFER 
TO THE SECOND THROUGH THE FIFTH 
SHIELD RESPECTIVELY. 

NRNGS 

B » C 


I WHEN EQUAL AREA ANNULI ARE CHOSEN 1 
FOR EACH SURFACE. 0 WHEN EQUAL 
RADII ANNULI ARE CHOSEN FOR EACH 
SURFACE. UNLESS A VALUE IS 
SPECIFIED A VALUE OF 1 IS ASSIGNED. 

THIS VARIABLE IS IGNORED FOR SURFACES 
OPPOSITE NON PLANAR SOURCES. 

THE NUMBER OF ANNULI ON EACH SURFACE. I 
NRNGS NORMALLY CANNOT EXCEED 20, BUT 


IT CANNOT EXCEED 10 IF NON PLANAR 
SURFACES ARE INVOLVED. 

THE NUMBER OF SHIELDS. THERE MAY BE I 
UP TO 5 SHIELDS IN THE SYSTEM. 

I WHEN THE HEATER SURFACE IS NON 1 

PLANAR. OTHERWISE 0. A VALUE OF 
0 IS ASSIGNED INITIALLY. 

THE EQUIVALENT OF NVUFH EXCEPT 1 

FOR THE TANK SURFACE. 

SOLAR FLUX PER UNIT TIME AND AREA. 1 

THE RADIUS OF THE HEATER IN THE 
LATERAL DIRECTION. 

THE RADIUS OF THE SHIELDS. NSHLDS 

THE RATIO OF THE RADIUS IN THE AXIAL 1 
DIRECTION TO THE RADIUS IN THE 


LATERAL DIRECTION FOR A NONPLANAR HEATER. 
NEED NOT BE SPECIFIED IF THE HEATER 
IS PLANAR ( NVUFH=0 ). 

THE EQUIVALENT OF RSH EXCEPT 1 

FOR THE TANK. 

THE RADIUS OF THE TANK IN THE 1 

LATERAL DIRECTION. 

STEFAN-BOLTZMAN CONSTANT. THE UNITS I 
OF THIS CONSTANT DETERMINE THE 
UNITS OF ALL OTHER DIMENSIONED 
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INPUT 


EFN 


SOURCE STATEMENT 


IFN(S) 


C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


VARIABLES. UNLESS SPECIFIED 
SIGMA IS IN ENGLISH UNITS AND HAS 
THE VALUE 1.713E-09 BTU/HR/FT**2/R**4 
SOLABS SOLAR ABSORPTIVITY 1 

SOLAR 0.0 WHEN NEITHER SOURCE REPRESENTS 
SOLAR RADIATION. L.O IF ONE SOURCE 
IS SOLAR RADIATION. I N I TALL Y ASSIGNED 


A VALUE OF 0.0 

STHICK THE THICKNESS OF EACH SHIELD NSFLDS 

TBABSL THE CONSTANTS FOR THE POLYNOMIAL 1-4 

GIVING ABSORPTIVITY AS A FUNCTION 
OF TEMPERATURE. THIS SET OF ONE 
TO FOUR CONSTANTS IS USED WHEN AN 
ABSORPTIVITY KEY IS DESIGNATED AS 1 
TBABS2 THE CONSTANTS FOR THE ABSORPTIVITY 1-4 

THROUGH POLYNOMIAL WHEN THE KEYS ARE EACH 

TBABSL DESIGNATED AS 2 THROUGH 10 

TBCONI THE CONSTANTS FOR THE POLYNOMIAL 1-4 

THROUGH GIVING THE THERMAL CONDUCTIVITY AS A EACH 


TBC0N5 FUNCTION OF TEMPERATURE. EACH SET 
OF FROM I TO 4 CONSTANTS IS USED 
DEPENDING ON THE VALUE OF NCOND 
FOR THE SHIELD. 

TBEMSI THE EQUIVALENT OF TBABSI THROUGH 1-4 

THROUGH TBABSL EXCEPT THAT THESE CURVES EACH 

TBEMSL ARE FOR EM I SS I V I T I E S . 

TOL THE MAXIMUM ALLOWABLE DIFFERENCE 

BETWEEN EACH TRIAL AND CALCULATED 
TEMPERATURE OF THE SHIELDS AT THE 
TIME OF CONVERGENCE. ASSIGNED A 


TH 

VALUE OF 0.01 

THE TEMPERATURE FOR 

EACH 

ANNULUS 

NRNGS 

TT 

OF THE HEATER. 

THE TEMPERATURE FOR 

EACH 

ANNULUS 

NRNGS 

TS 

OF THE TANK. 

THE TEMPERATURES FOR 

EACH 

OF THE 

NSHLDS+1 


SURROUNDINGS. UNLESS SPECIFIED 
THESE VARIABLES ARE ASSIGNED A 
VALUE OF 0. 


D 


D 


D 


D 


B 

B 

A 


VARIABLES READ 


I N 


THROUGH THE USE OF 


NAMELIST NAME 


INPUT2 


NAME DESIGNATION QUANTITY CODE 

AHS THE AREA OF EACH ANNULUS ON THE NR 1 E 

HEATER OR THE SURFACE OPPOSITE THE TANK. 

ATS THE AREA OF EACH ANNULUS ON THE TANK NR 1 E 

OR SURFACE OPPOSITE THE HEATER. 

LIST 1 IF ALL THE DATA FOR NONPLANAR 1 

SOURCES ARE TO BE LISTED. 0 IF THIS 
DATA ARE NOT LISTED. INITIALLY 
ASSIGNED A VALUE OF 0 

MX I DESIGNATION FOR TYPE OF SURFACES 1 E 

INVOLVED. 1 IF A NONPLANAR HEATER 
SEES A PLANAR TANK OR SHIELD. 

2 IF A NONPLANAR TANK SEES A PLANAR 
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INPUT 


EFN SOURCE STATEMENT 


IFN(S) 


NR I 
F 

FN 

FL 

RMHD 

RMTD 

XLOR 

CODE 

A 

B 

C 

D 

E 


SHIELD OR HEATER. 3 IF A NONPLANAR 
TANK SEES A NONPLANAR HEATER. 

THE NUMBER OF ANNULI ON EACH SURFACE. I 
VIEW FACTORS FROM ANNULI ON THE NR1**2 

TANK TO ANNULI ON THE HEATER. 

THE SET OF DATA FOR A GIVEN CONFIG- 1 
UR AT ION WHICH IS STORED ON TAPE 1 
CAN BE FOR SEVERAL DIFFERENT NUMBERS 
OF ANNULI. AT THE END OF EACH SET OF 
DATA FN IS SET EQUAL TO I 
SET EQUAL TO 1 AT THE END OF THE 1 

LAST SET OF DATA. THERE IS A MAXIMUM 
OF 30 SETS OF DATA. IF IT IS DESIRED 
TO READ IN ADDITIONAL DATA FOR NONPLANAR 
SOURCES, FL IS SET EQUAL TO 1.0 JLST 
PRIOR TO EACH RETURN TO TBVUF . 

THE RATIO OF THE RADIUS IN THE AXIAL 1 

DIRECTION TO THE RADIUS IN THE LATERAL 
DIRECTION FOR THE NONPLANAR HEATER. 

THE EQUIVALENT OF RMHD EXCEPT FOR 1 

THE TANK. 

THE SPACING RATIO (L/R) BETWEEN THE 1 
NONPLANAR SURFACE AND THE SURFACE NEXT 
TO IT. 


E,F 

G 


E, H 


E,H 


EXPLANATION OF CODES 


WHEN THERE IS SOLAR RADIATION PRESENT ONLY NSHLDS VALUES OF 
THESE VARIABLES NEED BE READ IN. HOWEVER, SINCE THE 
NUMBERING BEGINS WITH THE HEATER, THE VALUE STORED IN THE 
SECOND LOCATION REPRESENTS THE PROPERTY BETWEEN THE FIRST 
AND SECOND SHIELDS. 

IF THE VALUES FOR THESE VARIABLES ARE THE SAME ACROSS A 
SURFACE, THEN DENOTING THE SECOND VALUE AS ZERO CAUSES 
THE PROGRAM TO ASSIGN THE FIRST VALUE TO ALL OF THE VALUES 
FOR THE VARIABLE 

KEYS RELATE THE PROPERTIES OF AN ELEMENT OR SHIELD TO THE 
CONSTANTS GIVING THAT PROPERTY. THUS THE ENTRY NET(3)=4 
INDICATES THAT THE EMISSIVITY OF THE THIRD ANNULUS OF THE 
TANK IS DETERMINED BY THE VALUES IN TBEMS4. IF ALL OF THE 
ANNULI ON A SURFACE HAVE THE SAME KEY FOR A VARIABLE, 
ASSIGNING A VALUE OF 0 TO THE SECOND ANNULUS CAUSES THE 
PROGRAM TO ASSIGN THE VALUE FOR THE FIRST ANNULUS TO ALL OF 
ANNULI ON THE SURFACE. 

THE POLYNOMIALS ARE IN ASCENDING ORDER SO THAT THE FIRST 
COEFFICIENT IS A CONSTANT WHILE THE FOURTH COEFFICIENT IS 
MULTIPLIED BY THE TEMPERATURE CUBED. 

THE VIEW FACTOR DATA CONSISTS OF BETWEEN 1 AND 30 SETS OF 
DATA. EACH SET IS FOR A DIFFERENT SPACING RATIO OR SHAPE 
OF THE NONPLANAR HEATER OF TANK. WITHIN A SET OF DATA 
INFORMATION CAN BE PROVIDED FOR DIFFERING NUMBERS OF ANNULI 
FOR EACH SET OF DATA THERE IS A SUBSET FOR EACH NUMBER OF 
ANNULI. EACH SUBSET CONTAINS THE AREAS OF EACH ANNULUS ON 
THE NONPLANAR SURFACE AND THE SURFACE ADJACENT TO IT, AS 
WELL AS THE VIEW FACTORS BETWEEN ANNULI. 
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INPUT - EFN SOURCE STATEMENT - IFN(S) - 

F THE VIEW FACTORS ARE READ IN THROUGH A ONE DIMENSIONAL 

VECTOR ( F ) » THE VIEW FACTORS ARE STORED SERIALLY. THE VIEW 
FACTORS ARE FROM THE TANK SIDE TO THE HEATER SIDE. THUS 
FOR MX 1=3 THE VALUES OF F ARE VIEW FACTORS FROM ANNULI ON 
THE TANK TO ANNULI ON THE HEATER. THE FIRST VIEW FACTOR IS 
FROM THE INNERMOST ANNULUS ON THE TANK TO THE INNERMOST 
ANNULUS ON THE HEATER. THE SECOND VIEW FACTOR IS FROM THE 
INNERMOST ANNULUS ON THE TANK TO THE SECOND ANNLUUS ON THE 
HEATER. 

G ALL OF THE DATA FOR A CONFIGURATION SHOULD BE READ IN AS 

A BLOCK OF DATA. EACH BLOCK CAN CONTAIN DATA FOR SEVERAL 
DIFFERENT NUMBER OF ANNULI. THUS THE VIEW FACTORS ARE 
STORED IN A TRIPLY SUBSCRIPTED ARRAY WITH THE THIRD INDEX 
GIVING THE NUMBER OF ANNULI ON EACH SURFACE. 

H THESE VARIABLES ARE READ IN SO AS TO DISTINGUISH DIFFERENT 

TYPES OF SURFACES. FOR A HEMISPHERE THEY WOULD HAVE A VALUE 
OF 1.0, AND FOR AN OBLATE SPHEROID THEY WOULD BE LESS THAN 
1.0 


GO TO ( 1,2 ) , JR 

READ IN ALL DATA OTHER THAN THAT FOR THE VIEW FACTORS BETWEEN 
NONPLANAR SURFACES. 

1 RE AD ( 5, I N PUT 1 ) 

WRITEI6, 10 ) 

10 FORMAT ( 42 X , 18 HL I ST ING OF INPUT ) 

WRITEI6, OUTPUi) 

DO 11 1=1,10 
SUM 1 = 0 • 0 
SUM 2 = 0 . 0 
DO 12 J = 1 , 4 

SUM 1 = SUM 1 + CV A 0 S ( J , I ) 

SUM2=SUM2 + CVEM S ( J , I ) 

12 CONTINUE 

IF ( SUM 1 .EQ. 0.0 ) GO TO 14 

WR I TE ( 6 , 13 ) I DC I ) , ( C V A BS ( J , I ) , J= 1 , 4 ) 

13 FORMAT! 3X, 5HT BAB S , A 1 , 4E 1 2. 4) 

14 IF ( SUM2 .EQ. 0.0 ) GO TO 11 

WR I T E ! 6 , 1 5 ) I DC I ) , ( C V E MS { J , I ) , J = 1 , 4 ) 

15 FORMAT! 1 H+ , 66X , 5HT BE MS, A 1 , 4E 12 . 4 ) 

11 CONTINUE 

IF! ICONDR .NE. C .OR. NSHLDS .EQ. 0 ) GO TO 16 
DO 17 1=1,5 
SUM 1 = 0. 0 
DO 18 J = 1 , 4 

SUM 1 = SUM 1 + CVCQND! J, I ) 

18 CONTINUE 

IF ! SUM 1 .EQ. 0.0 ) GO TO 17 

WRITEI6, 19 ) I D( I ) , ( CVCOND! J, I ) , J=l f 4) 

19 FORMAT! 3X, 5HTBC0N, A1,4E12.4) 

17 CONTINUE 

WR I T E < 6 , 3 ) ( STHICK! I), 1=1, NSHLDS) 

3 FORMAT! 18H SHIELD THICKNESS= ,5E12.4> 

WR I T E ( 6 , 6 ) ! NCOND! I ), 1 = 1, NSHLDS) 
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n n 


INPUT 


EFN SOURCE STATEMENT 


I FN ( S ) 


6 FORMATt 19 H CONDUCTIVITY KEYS=,5I5) 

16 WRITE(6,4) ( RS( I I , 1*1, NSHLDSI 
4 FORMAT ( 16H SHIELC RADIUS =,5E12.4) 

NX = NSHL DS + 1 
NS= 1 

IF ( SULAR .NE. 0.0 ) NS = 2 
IF ( NS .GT. NX ) GO TO 21 

WR I T E { 6 » 22 ) ( I , 0(1) , ARES(I) , ABSR(I) , ES(I) , TS(I) , 

1 I = NS, NX ) 

22 FORMAT( //36X, 44HPR0P ERT l ES OF SURROUNDINGS BETWEEN SURFACES 
1 SPACING DISTANCE AREA ABSORPTIVITY EMISSIVITY 

2ERATURE / ( 16, 3 F 1 3 . 3, F 1 5 . 3 , F 1 3. 3 > ) 


21 

DO 

23 K= 1 , 6 



NF (K ) = NRNGS 



IF 

( K .LE. 3 .AND. 

SOLAR .NE 


GO 

TG ( 31,32,33,34, 

35, 36) ,K 

31 

IF 

( NE H ( 2 ) .EO. 0 ) 

NF ( 1 ) = L 


GO 

TO 23 


32 

IF 

( NAH ( 2 ) .EQ. 0 ) 

NF ( 2 ) =1 


GO 

TO 23 


33 

IF 

( TH ( 2 ) .EQ. C.O 

) N F ( 3 ) = 1 


GO 

TO 23 


34 

IF 

( NET ( 2 ) .EQ. 0 ) 

NF ( 4 ) =1 


GO 

TO 23 


35 

IF 

( NAT ( 2 ) .EQ. 0 ) 

NF ( 5)=1 


GO 

TO 23 


36 

IF 

( TT ( 2 ) .EQ. 0.0 

) N F t 6 ) = i 


23 CONTINUE 

IF ( SULAR .NE. 0.0 ) GO TO 27 
N 1 = N F ( 1 ) 

WR I T E ( 6 , 24 ) DNE(l) , ( N EH ( l ) , I = 1 , N 1 ) 

N1-NFC2 > 

24 FORMAT ( 12X,A6,20I3) 

WR I T E ( 6 , 24 ) DNA(L) , ( N AH (I ) , I = 1 , N 1 ) 

N 1 = NF ( 3 ) 

WRIF£(6,25 ) DTEM(l) , ( TH ( I ) , I = 1 , N 1 ) 

25 FORMAT ( 12X * A6» 10F11.3/18X, 10FL1.3) 

27 NX=NSHLDS*2 

IF ( NX .EQ. 0 ) GO TO 28 

DO 26 VI = 1 1 N X 

N1=NRNGS 

IF ( NE ( 2 , N ) .EQ. 0 ) Nl=l 

WR I T E ( 6 , 24 ) DNE ( N + 1 ) , ( NE ( I , N ) , I = 1 , N 1 ) 

N 1=NRNGS 

IF ( N A ( 2 , N ) .EQ. 0 ) Nl = l 

WR I T E ( 6 , 24 ) DNAIN + l) , ( NA(I,N) , 1=1, Nl) 

26 CONTINUE 

28 N 1 =NF ( 4 ) 

WRIFE(6t24) DNE (12) , ( NET { I ) , I = 1 , N 1 ) 

N1=NF ( 5 ) 

WRITE(6,24) D N A ( 12) , ( NAT ( I ) , I = 1 , N 1 ) 

Nl-NF { 6 ) 

WRITE(6,25 ) DTEM ( 2 ) , ( TT ( I ) , I = 1 , N 1 ) 

RETURN 

READ IN VIEW FACTORS AND AREAS FOR NCNPLANAR SURFACES 
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/76H 

TEMP 


I 


INPUT 


EFN SOURCE STATEMENT 


IFN(S) 


C 

2 READ ( 5, INPUT2 ) 

IF I LIST .EO. 0 ) RETURN 
WR ITE ( 6» 5 ) 

5 FORMAT! 1H1 , 42X t 3 CHL 1 ST I NG OF VIEW FACTOR DATA > 

WR IT E ( fa »OUTPU2 ) 

WR I T E ( 6 » 50 ) DAHS t I AHS ( J ) » J= 1 » NR 1 ) 

50 F0RMAT(6X»A6» LOF LI. 4) 

WR l TE ( 6 » 50 ) DATS , I ATS ( J ) , J= 1 , NR I ) 

WR I T E ( fa » 53 ) 

53 FORMAT! I2X, 30H VIEW FACTORS BETWEEN ANNULI ) 

NR2=NR1**2 

WR I T E ! 6 » 51 ) < F(N),N = 1.NR2) 

51 FORMAT! 10E12. 4) 

RETURN 

ENC 
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o o n n 


tIBFTC SHIELD DECK 


DIMENSION ESI 10 ) , ABSR(IO) , REFS(IO) , TS(IO) , ARESI10) 

EQUIVALENCE 

I I C VCOND , CVK ) , (CVEMS , CVX) , (CVABS « CVY) 

DIMENSION RSI9) , STHICK (7 ) , NC0NDI7) , DI8) , 

1 THI2C) , TTI20) , CVCOND(4,10) , CVABS(4,10) , 

2 CVEM S 1 4, 10 ) , NEH120) , NET120) , NAHI20) , NATI20) , 

3 CVXI40) , CVY I 40 ) , CVKI20) , NE(20,10) , NA(20,10) 
DIMENSION NEII20C) , NAK200) 

EQUIVALENCE I NEI , NE ) , I NA1 , NA ) 

COMMON / BULL/ EMS(20,9,2) , ABI20,9,2) , REF(20,9,2) 

COMMCN/ NEW/ FSSHTT 1 6 ) , FSSTTH1 6) , FSNHTT(20,6) , FSNTTH(20,6) f 
L FNNHTTI 20,20,6 ) , FNNTTH ( 20 , 20 , 6 ) , AREAI20.7) , 

2 F NS H T I ( 2 0 , 6 ) , FNSTTH(20,6) , R(20,7) , ARE AT ( 7 ) , 

3 FNHT SRI 20*6) , FNTTSR(20,6) , FSRTNH(20,6) , 

4 FSRTNT (20,6) , FSRSRI6) 

COMMON / EAS / AT , AH , QOA , T(20,10) , F13 , F14 , F23 , F24 
COMMON /CGNDT/ C0ND(20,I0) , RES(20,10) 

COMMON / KSLV / TGI 20, 10) , LOOK , LCOND , ICOND , NCALC , T MAX 
COMMON /T3/ Q0A2 , QSUR , QOAS 
DIMENSION R0RSI2C) 

COMMCN C I 1000 ) 

EQUIVALENCE 


1(C(1) 

y ABSR 

) 

9 

I C l 1 1 ) , ARES 

) 

* 

{ C ( 2 L ) 9 0 

> 


2 ( C ( 3 1 ) 

• ES 

) 

f 

I C I 41 ) , I VCOND) 

» 

( C ( 42 ) 9 I CONDR ) 


3 ( C ( A 3 ) 

t UMAX 

) 

1 

I C 1 5 1 ) , NAH 

) 

9 

( C ( 7 I ) 9 NAT 

) 


4 ( C ( 9 1 ) 

9 N A 

) 

t 

I C I 29 1 ) , NE H 

) 

9 

(C ( 311 ) 9 MET 

) 


5 ( C ( 3 3 1 ) 

9 NE 

) 

t 

I C I 53 1 ) , NCOND 

) 

9 

(C ( 541 ) f NCASE 

) 


6 ( C ( 54 2 ) 

, NEL I PS ) 

* 

1 C I 543 ) »NRADS 

) 

9 

( C ( 544 ) 9 NRNGS 

) 


7 ( C ( 545 ) 

, NS HL DS ) 

» 

(Cl 546 ) ,NCALCR) 

9 

( C ( 547 ) 9 NVUFH 

) 


8 ( C ( 548 ) 

9 NVUFT 

) 

f 

(C (549) f PHI 

) 

9 

( C ( 550 ) , P I 

) 


9 ( C (551 ) 

»RS 

) 

t 

( C ( 56 1 ) 9 RSH 

) 

9 

( C ( 562 ) 9 RST 

) 


A ( C ( 5 6 3 ) 

9 RH 

) 

9 

( C ( 564 ) 9 RT 

) 

9 

( C ( 565 ) 9 SOLAR 

) 


8 ( C ( 566 ) 

t S I GM A 

) 

9 

( C ( 567 ) , TOL 

) 

9 

{ C ( 568 ) 9 WI 

) 


C(C (571 ) 

, SOL ABS ) 

9 

< C ( 59 l ) 9 ST HI CK ) 

9 

(C(601), CVABS 

) 


D ( C ( 64 1 ) 

,CVCOND) 

9 

I C I 68 1 ) , CVEMS 

) 

9 

( C ( 721 ) f TH 

) 


E ( C ( 74 1 ) 

, TS 

) 

9 

I C I 75 1 ) , IT 

} 

9 

CC (771), REFS 

) 



COMMON / ELIP2 / FN , FL , RMHD , RHTD , XLOR , MXS , ATSIIO) , 
I AHSUO) , F I 10, 10 ) , LIST 

LOGICAL PASS 


INITIALIZE DATA AND GIVE INPUT VARIABLES THEIR ASSIGNED VALUES 


00 510 

1 = 1 

»200 




NEK I ) = 

= 0 





NAK I ) = 

= 0 





IF ( I 

-GT 

• 40 ) 

1 GQ 

TO 

510 

CVX( l ) 

= 0 

• 0 




C V Y ( I ) 

= 0 

.0 




IF ( I 

.GT 

• 20 ) 

' GO 

TO 

510 
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non 


SHIELD 


EFN 


SOURCE STATEMENT 


IFN(S) 


CV K ( I ) = 0.0 
NET ( I ) =0 
NEH ( I ) - 0 
NAT ( I ) = 0 
NAH ( I ) =0 
TH( I ) =0 . 0 
TT ( I ) =0 • 0 

IF( I .GT. 10 ) GO TO 510 
D( I ) = 0.0 
TS ( I ) =0 . 0 
ARES ( I ) = 1 . 0 
ESI I > = 1.0 
ABSR ( I ) = 1 • 0 
510 CONTINUE 
TOL = 0.01 
N I =0 
NR ADS = I 
NVUFT =0 
NVUFH - 0 
SOLAR =0.0 
I VCOND= 1 
I CCNCR= 1 
NC AS E = 1 
W I =0 • A 
NEL I PS = 0 
L I ST=0 
I T M A X = 200 
NC ALC R= 1 
PI = 3.14159 
SIGMA= 0. 1713 E— 08 
1000 WRITE (6,100) 

100 FORMAT I 1H1 ) 

READ INPUT DATA 

106 CALL INPUT ( 1 ) 

IF ( NELIPS .ECU 0 .AND. ( NVUFT .NE. 0 .OR. NVUFH . NE . 0 ) ) 

1 CALL TBVUF(O) 

WRITEI6, 107) 

107 F0RMATI42X, 18 ML I STING OF OUTPUT ) 

NE L I PS = 1 

TOLX = TOL 

IF! ICUNDR .NE. C ) ICOND-ICONDR 

IF ( IVCOND .EQ. 0 ) IVC0ND=1 

IF I ICONDR .EQ. 0 ) ICOND=IVCUND 

NCALC=NCALCR 

I C AS E = 0 

SDT=i.0E10 

I T COND= 1 

W=1.0 

LW = 0 

LOCK= l 

WX = Wi 

T EMPH = TH I 1 ) 

TEMPT= TT I l ) 

IF (NI.EQ.O) WRITE (6,100) 
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SHIELD - EFN SOURCE STATEMENT - IFN(S) - 

506 N I = 1 

DO 501 L= 1 » 2 
C 

C STORE DATA FOR EMISSIVITY KEYS AND TEMPERATURES 

C 

DO 501 M= 1 » 5 
J J=2*M + L-2 
DO 501 K= 1*8 

IF ( JJ .EG). 1 ) GO TO 520 
IF { K-7 ) 502 , 517 , 518 

520 GO TO ( 511 , 512 , 513 , 514 , 515 , 516 , 517 , 518 ) , K 

511 NS = TH ( 2 ) 

GO TO 519 

512 NS = TT (2 ) 

GO TO 519 

513 NS = N ET { 2 ) 

GO TO 519 

514 NS = N EH ( 2 ) 

GO TO 519 

515 NS = N AT ( 2 ) 

GO TO 519 

516 NS = N AH ( 2 ) 

GO TO 519 

517 NS = N E ( 2 T J J ) 

GO TG 519 

518 NS = N A ( 2 f J J ) 

519 DO 521 N=2 » 20 

IF ( NS .NE. 0 ) GO TO 502 

GO TO ( 531 , 532 , 533 , 534 , 535 , 536 * 537 , 538 ) , K 

531 TH ( N )=TH ( 1 ) 

GO TO 521 

532 TT ( N ) = T T ( 1 ) 

GO TG 521 

533 NE T ( N ) = NET ( 1 ) 

GO TG 521 

534 NEH ( N ) =NEH ( 1 ) 

GO TG 521 

535 NAT ( N ) =N AT ( 1 ) 

GO TO 521 

536 NAF ( N ) = N AH ( 1 ) 

GO TO 521 

537 NE (Nt JJ )=NE( It JJ ) 

GO FG 521 

538 NA(NtJJ )=NA( 1, JJ ) 

521 CONTINUE 
502 CONTINUE 
501 CONTINUE 

DO 540 N= 1 t NRNGS 
T { N» 1 ) = TH ( N ) 

T ( NfNSHLCS+2 )=TT ( N ) 

540 CONTINUE 

DETERMINE SURFACE PROPERTIES OF HEATER AND TANK 

DO 503 N=l, NRNGS 
J=NEH ( N ) 
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SHIELD 


EFN SOURCE STATEMENT 


IFN(S) 


EMS I N» 1 » 2 ) = CVEMS(L»J) + 1 CVEMSI2.J) + TH ( N ) * ( CVEMS(3,J) + 

1 CVEMS(A,J)*TH(N) ))*TH(N) 

J=NAH(N ) 

AB ( N * I » 2 ) = CVABS(ltJ) + ( CVABS(2,J) + THIN)*! CVABSI3,J) + 

1 CVABS(A,J )*TH( N > ) ) *TH ( N ) 

REFIN, 1,2) = 1.0 — AB ( N, 1,2) 

J=NET ( N ) 

EMSIN.NSHLDS+2, l ) = CVEMS(1,J) + ( CVEMS(2,J) + TTIN)*( CVEMS(3,J) 
1 + CVEMS(A,J)*TT(N) ) )*TTIN) 

J=NA T IN) 

AB(N,NSHLDS+2»1)= CVABS(1,J) + ( CVABS(2,J) + TT(N>*( CVABS(3,J)+ 

1 CVABSI A, J >*TT(N) ) )*TT(N) 

REFIN, NSHLDS+2, 1 )= 1 . 0-AB ( N ,N SHLDS+2 , 1 ) 

503 CONTINUE 

DO 5 A 1 M= 1,10 
REFS(M)=1.0-ABSR (M) 

5A1 CONTINUE 

CALCULATE VIEW FACTORS BETWEEN PAIRS OF SURFACES 
CALL VUFAC1 

IF ( NSHLDS .EO. 0 ) GO TO 10 
IF ( ICONOR .NE. 0 ) GO TO AOO 

CALCULATE THERMAL RESISTANCE BETWEEN ANNULI 

NF = NRNGS - 1 
DO 600 M = l,. NSHLDS 
RESK=SlGMA/STHICK(M)/2.0/PI 
DO 601 N = 1 * NF 

IF(N.EQ.l) RESIN, M)=RESK*ALOG( (R(2,M+1)+R( 1 » M+ 1 ) )/R( l,M+l ) ) 

IF (N.GT . I ) RESIN »M)=RESK*ALOG( (RIN+1,M + 1) + R(N,M+1))/ 

1 (RIN,M+1)+RI N— 1,M+1) ) ) 

601 CONTINUE 
600 CONTINUE 

MAKE INITIAL ESTIMATES OF TEMPERATURE DISTRIBUTIONS 

AOO IF I ICASE .GT. C ) GO TO 20 

IF I SULAR .EQ. 0.0) TEMPX = TEMPH 
IF I SOLAR .NE. C.O) TEMPX = 500.0 
Til, 2) = TEMPX -I TEMPX-TEMPT )/FLOAT (NSHLDS+1 ) 

IF I NSHLDS .EQ. I ) GO TO 20 

NSF=NSHLDS+1 

DO 15 M = 3 » N S F 

T(l,M) = T ( l , M— l ) - (TEMPX -111,21) 

15 CONTINUE 

CALCULATE SURFACE PROPERTIES AND THERMAL CONDUCTIVITIES FOR 
SHIELDS 

20 DO 21 M = 1, NSHLDS 
JK = NCONDIM) 

K = M + 1 

DO 21 L = 1,2 
JJ = 2*M + L - 2 
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SHIELD - EFN SOURCE STATEMENT - IFN(S) - 

DO 21 N = l» 20 

IF ( LOOK .NE. I .AND. ICASE .NE. 0 ) GO TO 17 
T(N,K)=T( 1,K) 

TO (N t K ) = T t N» K ) 

17 IF ( N .GT. NRNGS ) GO TO 21 
J=NE(N, JJ) 

EMS(N,K,L)=CVEMS(lt J )+CVEM$(2, J ) *T ( N , K ) +C VEMS ( 3 , J)*T(N,K)**2 
1 +CVEMS 1 4» J )*T ( N» K ) 

J = NA(N,JJ) 

AB(N,K,L) =CVABS( 1, J )+CVABS(2, J ) *T < N t K ) +CV ABS ( 3 , J)*T(N,K)**2 
1 +CVABS (4, J )*T(N,K)**3 

REF(N,K,L) = 1.0 - A8 ( N » K » L ) 

IFIN. EQ. NRNGS. OR. ICOND. NE.O) GO TO 21 
TMS= ( T ( N »K ) + T(N+1,K) )/2.0 

COND I N» M ) = CVC0N0(ltJK)+CVC0ND(2»JK)* TMS+C VCOND (3»JK)*TMS**2 + 
1 C VCOND ( 4 * JK ) *TMS**3 

21 CONTINUE 
TMAX=0. 0 

CALCULATE SHIELD TEMPERATURES 

201 CALL SHCLC 

IF ( TM AX .NE. 0.0 ) GO TO 7 

202 DO 210 M = 1 i NSHLDS 
K=M+ 1 

DO 210 N * 1, NRNGS 

IF ( ICOND .EQ. 2 ) T(N,K) = T(1,K) 

210 CONTINUE 

TEST TO SEE IF TEMPERATURES HAVE CONVERGED 


SDTX=0 . 0 
P ASS= .T RUE . 

DO 6 M = 1, NSHLDS 
K=M+1 

NF = NRNGS 

I F ( ICOND .EQ. 2 ) NF = 1 
DO 6 N = 1,NF 

TEST = ABS (TO(N» K ) - T(N,K) ) 

IF ( TEST .GT. TOLX ) PASS=. FALSE. 
SDTX=SDTX+TEST 
6 CONTINUE 

IF I PASS .AND. .TRUE. ) GO TO 9 

ICASE=1 

LOOK=LOOK+ 1 

IF I LOOK .GT. ITMAX ) GO TO 14 
IF ( ICOND .NE. C ) GO TO 8 
IF ( SDTX .GT. SCT ) GO TO 13 
SDT-SDTX 
GO TO 8 
13 SDT=1.0E10 
LW=LW+1 

IF I LW .EQ. 8 ) GO TO 7 
ICOND^I CONDL 
W= 1 . 0 
ICASE=0 
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SHIELD - EFN SOURCE STATEMENT - IFN(S) - 

GO TO 400 

8 DO L M= 1 , NSHLDS 
K = M + 1 

DO 2 N= 1 ♦ NRNGS 

IF { ICOND .EQ. 2 .AND. N -GT. I ) GO TO 2 
T(N,K) = TC( N,K > -W*(T0(N,K)-T<N,K) ) 

TO(N,K)=r (N*K) 

2 CONTINUE 
I CONTINUE 
GO TO 20 

7 GO TO ( 22,23,24,25 ), ITCOND 

22 ICON D= I VCQND 
GO TO I L 

23 ICCND=3-IVC0ND 
GO TO II 

24 ICGND=3~IVC0ND 
GO TO II 

25 WRITE(6,101) NCALC , ICOND 

101 FORMAT ( 28H CASE ABANDONED WITH NC ALC= , l 3 , 6 X , 6H I COND= , I 3 ) 

19 M = NS HLDS + 1 

WRITE (6, 109) U T 0 ( N , K ) ,N=1, NRNGS) ,K=2,M) 

109 FORMAT ( 30 H PREVIOUS TEMPERATURE VALUES /(10F12.3) > 

NC AS E=2 
I CAS E=0 
GO TO 10 

11 ITCCND= ITCOND+1 
I C ASE=0 

W= 1.0 

NC ALC= 1 — NCALC 
LW=0 

WRITEI6, 12) NCALC , ICOND 

12 FORMAT ( 44H CASE DIVERGING, CHANGING NCALC AND ICOND ,L2X, 

1 6HNC ALC= ,11, 1 2X , 6H I CGND= , I 1 ) 

GO TC 400 
14 WR I T E { 6 , 18 ) 

18 FORMAT ( 36H EXCESSIVE NUMBER OF ITERATIONS ) 

GO TG 19 

9 IF ( ICOND .EQ. ICQNDR ) GO TO 10 
IF ( LW .EQ. 0 ) GO TO 26 
W^WX/2.0 

I C AS E= 1 
I CCNC= I CON DR 
GO TO 20 

26 I C G N D L = ICOND 
I CCND= I CGNDR 
W= W I 

I C AS E= 1 
GO TO 20 

DETERMINE HEAT TRANSFER RATE FOR TANK 
10 CALL QDOT 

WRITE OUTPUT 

IF ( NSHLDS .NE. 0 ) GO TO 105 
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WR I TE ( 6 » 108 ) SOLAR , QO A , Q0A2 , QOAS 
C 

C LISTING OF OUTPUT VARIABLES 

C 

C SOLAR SAME VARIABLE AS IN INPUT LIST 

C QUA HEAT TRANSFER RATE PER UNIT AREA TC SOURCE DESIGNATED AS 

C THE TANK. 

C • Q0A2 HEAT TRANSFER RATE PER UNIT AREA TO SURFACE OPPOSITE TO 
C THE TANK. 

C QOAS HEAT TRANSFER RATE PER UNIT AREA TO THE SURROUNDINGS 

C BETWEEN THE TANK AND THE OPPOSITE SURFACE. 

C LOOK THE NUMBER OF ITERATIONS 

C TOLX THE MAXIMUM DIFFERENCE BETWEEN EACH ESTIMATED AND 

C CALCULATED TEMPERATURE FOR THE SHIELDS. 

C M THE SHIELO NUMBER. 

C T THE TEMPERATURE OF EACH ANNULUS OF EACH SHIELD. 

C R ( N ) THE RADIAL DISTANCE TO THE OUTER EDGE OF THE NTH ANNULUS. 

C RORS R ( N ) DIVIDED BY THE RADIUS OF THE SHIELD ( RS ). 

C 

L08 FORMAT! IH /12H SOLAR =,F3.1,6X,6H QOA = , E I 2 .4 , 6X , 6H00A2 =, 

1E12.4»6X, 6HQ0AS =,E12.4) 

110 FORMAT ( 1HL ) 

GO TC 106 

105 WR ITE ( 6 » 102 ) LOOK , SOLAR , QOA , TOLX 

102 FORMAT! /// 17H NO. ITERATIONS =,I4,6X,8H SOLAR =,F3.0,10X, 

1 7H QOAT = , E12.4, 10X,6H TOLX=,F7.5) 

WR ITE ( 6 » 1 1 0 ) 

DO 16 M= 1 » NSHLDS 
K=M+ 1 

DO 204 N = 1 » NRNGS 

RORS(N) = R(N»M+1)/R( NRNGS »M+1 ) 

204 CONTINUF 

WR ITE ( 6 f 103 ) M , T < l , K ) , R(1,K) , RORS(l) 

103 F0RM.AH//48H M TEMP R ( N ) R(N)/RS / 

1 I6» 3F12.3) 

IF ( NRNGS .EQ. 1 ) GO TO 16 

WR I T E ( 6 , 104) ( T ( N » K ) , R ( N » K ) , RCRS(N) , N = 2 , NRNGS ) 

104 FORMAT! F18.3,2F12. 3) 

16 CONTINUE 

IF (NCASF .NE. 1 ) GO TO 1000 

NCASE = 2 

TOLX = TOLX/2 .0 

ICASE = 1 

L0CK=2 

I T COND= 1 

WRITE (6,100) 

GO TO 20 
STOP 

ERROR MESSAGE NUMBER 1 
END 
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SIBFTC TBVUF 


DECK 


SUBROUTINE TBVUF { JR) 

THIS ROUTINE READS IN AND STORES DATA FOR VIEW FACTORS AND AREAS 
OF NCNPLANAR SURFACES. 

A MAXIMUM OF 10 ANNULI CAN BE USED WHEN THERE ARE NCNPLANAR 
SURFACES, AND THE NUMBER OF ANNULI USED SHOULD CORRESPOND TO THE 
NUMBER OF ANNULI FOR WHICH DATA HAS BEEN STORED. 

COMMON / ELIPT / M , XOR , MX1 

COMMON/ NEW/ F SSHT T ( 6 ) , FSS TTH ( 6 ) , FSNHTT(20,6) , FSNTTH(20,6) , 

1 FNNHTT ( 20,20,6) , FNNT TH ( 2 0 , 20 , 6 ) , AREA(20,7) , 

2 FNSHTT (20,6) , FNSTTH(20,6) , R(20,7) , AREAT ( 7 ) , 

3 F NHT SR ( 20,6) , FNTTSR(20,6) , FSRTNH(20,6) , 

4 FSRT NT (20,6) , FSRSR(6) 

DIMENSION XT ( 30 ) , MXI30) , RMHS ( 30 ) , RMT S ( 30 ) , NRG(10,30) 

DIMENSION FT ( 100 0 ) , FX(10,10,10) 

DIMENSION ASTORE (10, 10,2) , ASY(200) 

EQUIVALENCE ( FT , FX > 

EQUIVALENCE ( ASY , ASTORE ) 

LOGICAL PASS 
COMMON C( 1000) 

EQUIVALENCE 


1 (C ( 1 ) 

, A B SR 

) 

t 

(C( 11 ) 

, ARES 

) 

» 

(C (21) 

♦ D 

1 


2 ( C ( 3 1 ) 

, E S 

) 

t 

( C ( 4 1 ) 

, IVCOND) 

f 

(C (42) 

, ICONDR) 


3 ( C ( 4 3 ) 

, ITMAX 

) 

, 

( C ( 5 1 ) 

, NAH 

) 

f 

<C( 71) 

, N AT 

) 


4(C (91 ) 

, N A 

) 


(C ( 291 ) 

, NEH 

) 

t 

(C ( 311 ) 

» NET 

) 


5 ( C ( 3 3 1 ) 

f NE 

) 

, 

( C ( 53 1 ) 

, NCOND 

) 


(C ( 541 ) 

, NC A SE 

) 


6 ( C ( 5 42 ) 

* NEL I PS ) 

f 

(C( 543 ) 

, NR ADS 

) 

» 

(C (544) 

, NRNGS 

) 


7(C (545) 

t NSHLDS ) 

t 

( C ( 546 ) 

, NC ALC R ) 

* 

( C ( 547 ) 

, NVUFH 

) 


8(C (548) 

, N V U F T 

) 

t 

(C( 549) 

, PH I 

) 

t 

(C ( 550 ) 

,PI 

) 


9 ( C ( 5 5 1 ) 

,RS 

) 

t 

(C( 561 ) 

, RS H 

) 

, 

(C ( 562 ) 

, R ST 

) 


A ( C ( 563) 

, RH 

) 

* 

( C ( 564 ) 

,RT 

) 

, 

(C ( 565 ) 

, SOLAR 

) 


B(C (566) 

, SIGMA 

) 

t 

( C ( 56 7 ) 

, TOL 

) 

, 

(C (568 ) 

tWI 

) 


C ( C ( 5 7 1 ) 

.SOLABS) 

f 

( C ( 59 1 ) 

, STHICK ) 

, 

( C ( 60 1 > 

, C VA B S 

) 


D ( C ( 6 4 1 ) 

tCVCOND) 

f 

{ C ( 68 1 ) 

» CVEMS 

) 

t 

(C ( 72L ) 

, TH 

) 


E( C ( 741 ) 

,TS 

) 

t 

(C( 751 ) 

t TT 

) 

t 

(C ( 77L ) 

, REF S 

) 



COMMON / ELIP2 / FN , FL , RMHD , RMTD , XLOR , MXS , ATS(IO) , 

1 AHS ( 10) , F ( LOO) , LIST , NR 1 

IF ( JR .NE . 0 ) GO TO 10 
REWIND 1 
1 = 1 

12 FN = 0 • 0 
FL = 0 « 0 

READ VIEW FACTOR DATA. 

CALL INPUT (2) 

CONVERT THE VIEW FACTORS WHICH ARE STORED SERIALLY INTO THE PROPER 
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TBVUF 


EFN SOURCE STATEMENT 


I FN ( S ) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


c 

c 

c 


c 

c 

c 


c 

c 

c 


ARRAY FORM • THE FIRST VIEW FACTOR IS THE ONE BETWEEN THE INNERMOST 
ANNULUS ON THE SURFACE CLOSER TO THE TANK TO THE INNERMOST ANNULUS 
ON THE SURFACES CLOSER TO THE HEATER. THE NEXT VIEW FACTOR IS FROM 
THE INNERMOST ANNULUS ON THE TANK SIDE TO THE SECOND ANNULUS ON 
THE HEATER SIDE. AFTER NRNGS VIEW FACTORS HAVE BEEN READ IN, THE 
VIEW FACTOR FROM THE SECOND ANNULUS ON THE TANK SIDE TO THE 
INNERMOST ANNULUS ON THE HEATER SIDE IS READ IN. THE PROCESS IS 
REPEATED UNTIL ALL THE VIEW FACTORS BETWEEN A PAIR CF SURFACES ARE 
DETERMINED. 

NT=NR1**2 
DO 21 M= 1 , NT 
J= 1 + ( N-D/NR1 
K=N-( J-l )*NR1 
FX { K , J , NR I ) = F ( N ) 

21 CONTINUE 

DO 23 N= I, NR 1 
ASTORE( NtNRl, 1 ) = AHS ( N ) 

ASTOREINtNRl t 2)=ATS(N) 

23 CONTINUE 

NRG ( NR1 , I )=NR1 

IF ( FN .NE. 0.0 ) GO TO 22 

GO TO 12 

22 XT ( I ) =XLOR 
MX ( I ) =MXS 
RMHS ( l ) = RMHD 
RMTS ( I )=RMTD 

STORE VIEW FACTORS AND AREAS ON TAPE. 

W-RITE(l) XL OR , MXS , ( FT ( J ) , J = 1 , 1 000 ) , { A S Y ( J ) , J= 1 , 200 ) 

IF ( FL .NE. 0.0 ) RETURN 
1=1 + 1 
GO TO 12 

CHECK TO SEE IF SPACING AND SURFACES ARE IN TABLES 


10 L = 0 

DO 13 K= 1 , I 
PASS = .TRUE. 

IF ( { MX 1 .EQ. 1 

1 ABS ( RS H-RMHS (K ) > 

IF ( { MX 1 .EQ. 2 

1 ABS { RS T-RMT S ( K ) ) 


.OR. MX 1 .EQ. 3 
.GT. 0.05*RSH ) 
.OR. M X 1 . EQ . 3 
.GT. 0.05*RST ) 


( ABS ( XOR-XT ( K ) ) .GT. 0.05*X0R 


13 


IF 

IF ( NRNGS .NE. 

IF ( PASS .AND. 
CONTINUE 

IF { L .EQ. 0 .OR 
REWIND 1 


) .AND. 

PASS=. FALSE. 

) .AND. 

PASS=. FALSE. 

) PASS=. FALSE. 


NRG ( NRNGS* K ) 
.TRUE. ) L=K 


) PASS=. FALSE, 


NRNGS .GT. 10 ) GO TO 14 


READ VIEW FACTORS FROM TAPE 


DO 15 K = 1 , L 

READU) XL OR , MXS , ( F T ( J ) , J= 1 , 1 COO ) , { A S Y ( J ) , J= 1 , 200 ) 

15 CONTINUE 
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C 

C STORE VIEW FACTORS AND AREAS IN CORE 

C 

AREATIM )=0.0 

AREAT(M+1)=0.0 

DO 17 N=1 * NRNGS 

DO 16 NN* 1* NRNGS 

FNNTTH(NN,N,M) = FXCNNtNt NRNGS ) 

16 CONTINUE 

AREA ( N» M ) = ASTORE { N t NRNGS t 1 ) 

AREA (N» M + l )-ASTORE( Nt NRNGS t 2 J 
AREAT(M)=AREAT(M ) + AREAIN.M) 

AREAT (M + l ) = AREATIM + ll + AREA! N » M+ 1 ) 

17 CONTINUE 
C 

C READJUST RADII OF PLANAR SURFACES TO AGREE WITH INPUT DATA 

C 

I F ( NSHLDS .EQ. 0 ) GO TO 18 

IF { M .EQ. 1 ) R(N, 2)=SQRT(AREA(N*2) /PI I 

IF { M. EQ. NSHLDS+1 ) R ( N , M ) =SQR T < ARE A ( N » M > / P I) 

18 RETURN 

14 WR ITE ( 6 t 20 ) 

20 FORMAT ( 42H NO SPACING ENTRY IN TAOLE ) 

ASTOP=SQRT (-2.0 ) 

RETURN 

END 
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$IBFTC VUFAC1 DECK 


SUBROUTINE VUFACI 

THIS ROUTINE DETERMINES VIEW FACTORS BETWEEN ANNULI AND BETWEN 
WHOLE SURFACES. 

DIMENSION ESI 10 ) , ABSR(IO) , REFS(IO) , TS(IO) , ARESUO) 

COMMON / ELIPT / M , XOR , MX1 
DIMENSION RS(IO) , Dill) 

COMMON/NEW/ FSSHTT ( 6 ) , FSSTTHI6) , F SNHTT (20*6) , FSNTTH(20,6) , 

1 FNNHTT (20,20,6) , FNNTTH ( 20 , 20 , 6 ) , ARFA(20,7) , 

2 FNSHTT (20,6) , FNSTTHI20.6) , R(20,7) , AREAT ( 7 ) , 

3 FNHT SRI 20,6) , FNTT SR { 20, 6 ) , FSRTNH(20,6) , 

4 FSRTNT ( 20» 6 ) , FSRSRI6) 

COMMON / EAS / AT , AH , QOA , T(20,10) , F13 , F14 f F23 , F24 

COMMON CllOOO) 

EQUIVALENCE 

llC(l) , ABSR ) , (C(ll) * ARES ) , ( C ( 2 L ) ,D ) , 

2 1 C ( 3 1 ) ,ES ) , ( C ( 41 ) * I VCOND ) , C C ( 42 ) .ICONDR) , 

3 ( C ( 4 3 ) , I TM AX ) , (C(51) ^ NAH ) , ( C ( 7 1 ) > NAT ) , 

4ICI91) , NA ) , ( C ( 29 1 ) » NEH ) , (C(311),NET ) , 

5(0(331), NE ) , ( C ( 53 l ) , NCOND ) , (C ( 54 1 ) , NC A SE ) , 

6(C (542) .NELIPS) , ( C ( 543 ) , NR ADS ) , ( C ( 544 ) , NRNGS ) , 

7(C<545) .NSHLDS) , ( C ( 546 ) , NC ALCR ) , (C I 547 ) ,NVUFH ) , 

8 ( C ( 548 ) , NVUFT ) , (C(549),PHI ) , (C(550),PI ) , 

9(C(551),RS ) , ( C 1 56 1 ) , RSH ) , (C(562),RST ) , 

A ( C ( 563 ) , RH ) , ( C ( 564 ) , RT ) , ( C ( 565 ) , SOL AR ) , 

B(C(566> , SIGMA ) , (C(567),T0L ) , (C(568),WI ) , 

CICI571) ,SOLABS) , ( C ( 59 1 ) , STHICK ) , (C ( 601 ) ,CVABS ) , 

D(C (641) ,CVCQND) , ( C ( 68 1 ) ,CVEMS ) , (C(721),TH ) , 

E(C (741) ,TS ) , ( C ( 75 1 ) , TT ) , (C(771),REFS ) 

DETERMINE THE RACII OF EACH ANNULI OF EACH SURFACE 

FN=NRNGS 

IF ( NR ADS .NE. C ) FN= SQRT ( FN ) 

NSP=NSHL0S+2 
DO 101 J = 1 , NS P 
SUMA=0.0 

DO 102 N=l, NRNGS 
FN I = N 

IF ( NR AOS .NE. 0 ) FN I =SQRT ( FN I ) 

IF ( J .NE. 1 ) GO TO 120 

IF ( SOLAR .NE. 0.0 ) GO TO 101 

R(N, J) = RH=f‘FNI/FN 

IF ( NVUFH .EQ. 0 ) GO TO 121 

GO TO 130 

120 IF ( J .EQ. NSP ) GO TO 122 
R(N, J)=RS( J-l >*FNI/FN 

IF ( J .NE.NSP-1 .OR. J .NE. 2 ) GO TO 121 

IF ( J .EQ. NSP— l .AND. NVUFT .EQ. 0 ) GO TO 121 
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IF ( J . EQ . 2 .AND. NVUFH .EQ.O ) GO TO 121 
GO TO 130 

121 IF ( N .EQ. 1 ) AREA! 1, J )=PI*R( l f J)**2 

IF ( N .NE. 1 ) A R E A ( N ? J )=PI*( R(N,J)**2 - R(N~1,J)**2 ) 
GO TO 130 

122 R(N*J)=RT*FNI/FN 

IF ( NVUFT • EQ. C ) GO TO 121 

130 SUMA= SUM A + ARE A ( N » J ) 

102 CONTINUE 

CALCULATE THE AREA OF THE PLANAR SURFACES 

IF { SOLAR • NE* 0.0 .OR. J .NE. 1 ) GO TO 131 

AH=SUMA 
AREAT ( 1 ) = A H 
GO TO 101 

131 AREAT { J ) =SUM A 

IF { J .EQ. NSP ) AT = AR E AT ( J ) 

101 CUNTINUE 

MF = NSHL DS + 1 
DO 103 M= 1 » MF 

IF ( M .EQ. 1 .AND. SOLAR .NE. 0.0 ) GO TO 103 
DETERMINE IF NON PLANAR SURFACES ARE INVOLVED. 



I F 

( 

M .NE. 

1 . 

AND 

. M 

• NE 

. MF 

) GO 

TO 

300 


IF 

( 

NSHLDS 

.NE 

. 0 

) GO 

TO 

301 





IF 

( 

NVUFT 

.EQ. 

C 

.AND. 

NVUFH 

.EQ. 

0 ) 

GO TO 


IF 

( 

NVUFT 

.NE. 

C 

. AND . 

NVUFH 

.NE. 

0 ) 

MX 1 = 3 


IF 

( 

NVUFH 

.NE. 

0 

.AND. 

NVUFT 

.EQ. 

0 ) 

MX1 = 1 


IF 

( 

NVUFT 

.NE. 

c 

-AND. 

NVUFH 

. CO. 

0 ) 

M X 1 =2 


GO 

TO 

302 









301 

IF 

( 

M .NE. 

1 ) 

GO 

TO 303 






IF 

( 

NVUFH 

..EQ. 

c 

) GO 

TO 

300 





MX 1 

= 1 











GO 

TO 

302 









303 

IF 

( 

NVUFT 

.EQ. 

G 

) GO 

TO 

300 





MX 1-2 

EXTRACT DATA FOR VIEW FACTORS FROM STORAGE. 

302 IF ( M .EQ. 1 ) MT = L 
IF ( M .EQ. MF) M T =MF 
XOR= C ( M T )/R(NRNGS,MT+l) 

CALL TBVUF(l) 

DO 10 5 N= 1 f NRNGS 

SUMH 1=0 • 0 

DO 132 NN= l » NRNGS 

FNNHTTl NtNNtM)=FNNTTH(NN f N,M ) *AREA < NN f M+ 1) /AREA (N, M) 
SUMH 1= SUMH 1 + FNNHT T ( N » NN , M ) 

132 CONTINUE 

FNSHTT ( N»M )=SUMH 1 

FSNTTHt Nr M ) = FNSHTT (N,M) *AREA (N»M )/AREAT(M+l ) 

FNHTSRI N,M)= 1 . 0- FN S HT T ( N , M ) 
FSRTNH(NtM)=FNHTSR(N*M)*AREA(N*M)/ARES(M) 

105 CONTINUE 


1 
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CALCULATE THE VIEW FACTORS FROM AN ANNULUS TO A WHOLE SHIELD. 

112 DO 133 N=1 1 NRNGS 
SUMT 1=0 . 0 
DO 134 NN= 1 f NRNGS 
SUMT 1 = SUMT 1 + FNNTTH ( N> NN» M ) 

134 CONTINUE 

FNSTTH{N»M)=SUMT 1 

FSNHTTI N»M ) = FNST T H ( N f M ) * ARE A ( N , M+l )/AREAT(M ) 

FNTT SR ( N » M ) = 1.0 - FNSTTH ( N # M ) 

FSRTNT ( N» M ) = FN T T SR ( N , M ) * ARE A { N T M + 1 ) / ARES ( M ) 

133 CONTINUE 
GO TC 304 

CALCULATE THE VIEW FACTORS BETWEEN PLANAR ANNULI. 

300 DO 106 N=l, NRNGS 

IF ( N .EQ. 1 ) GO TO 107 

CALL SUBVF ( 2 , R ( NRNGS *M+i ) f 0.0* R<N-1«M)« R(N,M), D ( M ) ) 

FSNTTHtN* M ) = F 14 

FNSHTH N f M)=FSNTTH(N t M)*AREAT(M+l ) /AREA ( N f M ) 

GO TO 108 

107 CALL SUBVF ( 1 * R ( NRNGS » M + 1 ) f 0.0 t R(l f M) , 0.0 , D ( M ) ) 

FSNTTH ( l f M)=F13 

FNSHTT ( 1 1 M ) =FSNT TH ( 1 » M ) * AR E A T { M + 1 ) /AREA! 1, M) 

108 FNHTSR( N,M )=1.C - FNSHT T ( N »M ) 

FSRTNH(N,M) = FNHTSR ( N, M )* AREA ( N f M }/ ARES ( M ) 

DO 109 NN= 1 » NRNGS 

IF ( N .EO. 1 ) GO TO 110 

IF ( NN .EQ. 1 ) GO TO 111 

CALL SUBVF ( 4 , R ( NN- 1 f M+ 1 ) T R(NN»M+1) , R ( N- 1 1 M ) f R{N,M) , D ( M ) ) 
FNNTTH (NNfN f M)=F24 

FNNHTT ( N f NN f M ) = F NNTT H ( NN , N , M )*AREA (NN f M + l ) /AREA IN, M ) 

GO TC 109 

111 CALL SUBVF ( 2 , R(1,M+1) , 0.0 , R ( N- 1 t M ) f R(N,M) f DIM) ) 

FNNTTH ( 1, N f M)=F14 

FNNHTT ( N y 1 t M ) =FN NT TH ( 1 * N , M )* ARE A ( 1 , M+ 1 ) / ARE A ( N , M ) 

GO TO 109 

110 IF ( NN .EQ. 1 ) GO TO 113 

CALL SUBVF ( 3 , R(NN-l f M+l) , R(NN,M+1) f R(l f M) t 0.0 » D ( M ) ) 

FNNTTH ( NN ? 1 » M ) = F23 

FNNHTT ( 1,NN,M )= FNNTTH (NNf 1 # M)*AREA(NN t M+l) /AREA (1»M) 

GO TO 109 

113 CALL SUBVF ( 1 t R(l f M+l) , 0.0 , R(1,M) , 0.0 T D ( M ) ) 

FNNT TH ( 1, 1,M)=F13 

FNNHTT ( 1 » 1 » M ) =FNNTTH ( 1 , 1 , M ) * ARE A ( 1 , M+ 1 > / ARE A ( 1 f M ) 

109 CONTINUE 
106 CONTINUE 

GO TO 112 
C 

C CALCULATE VIFW FACTORS BETWEEN WHOLE SURFACES. 

C 

304 SUMSR 1 = 0 . 0 
SUMS R2 = 0 • 0 
SUMT 3=0.0 
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SUMH3=0 . 0 

DO 114 N= 1 » NRNGS 

SUMSR 1= SUMSR1 + FSR TNH ( N» M ) 

SUMSR2=SUMSR2 + FSRTNT(N.M) 

SUMT3=SUMT3 + FSNTTH(N»M) 

SUMH3=SUMH3 + FSNHTT(N»M) 

114 CONTINUE 

FSRSR(M) = 1.0- SUMSR 1-SUMSR2 
FSSTTHI M)=SUMT3 
FSSHTT ( M ) = SUMH3 
103 CONTINUE 
RETURN 
END 


75 


L 


non non non non 


$ I BFTC SUB VF 


DECK 


SUBKUUr I ME SUBVF ( I J , R 1 A , R2 A, R3A , R4A, ALA ) 

COMMON / EAS / AT . AH , QOA , T(20»LO> , F13 , F14 , F23 , F24 

DOUBLE PRECISION Rl , R2, R3, R4 ,AL ,C 1 , C2 , C3,C4 ,C5, C6,C 7, C8 , C9 

R 1 = R 1 A 

R2=R2 A 

R3=R3A 

R4=R4 A 

AL=AL A 

GO TO ( 1,2, 3,4), IJ 

CALCULATE THE VIEW FACTOR BETWEEN TWO DISKS 

1 Cl = (1.+(R3**2+AL**2)/R1**2)**2-4.*(R3**2)/R1**2 
F I 3 = ,5*( l. + (R3**2 + AL**2)/RI**2-SQRT(CI ) ) 

RETURN 

CALCULATE THE VIEW FACTOR BETWEEN A DISK AND AN ANNULUS 

2 C2 = (L.+(R4**2+AL**2)/R1**2)**2-4.*(R4**2)/RL**2 

C3 = ( 1 . + ( R3**2 + AL**2 )/R 1**2 )*#2-4 ( R 3**2 ) /R 1**2 

F 1 4 = . 5* ( (R4**2-R3**2)/R1**2-SQRT(C2)+SQRT(C3) ) 

RETURN 

CALCULATE THE VIEW FACTOR BETWEEN AN ANNULUS AND A DISK 

3 C4=(R2**2+R3**2+AL**2 ) **2- ( 2 .*R2*R3 ) **2 

C5 = (R1**2+R3**2+AL**2)**2-(2.*R1*R3)**2 

F23 = ( .5/ (R2**2-Rl**2) ) * ( R2**2-SQRT ( C4 ) +SQR T ( C 5 )-R 1**2 ) 

RETURN 

CALCULATE THE VIEW FACTOR BETWEEN TWO ANNULI 


4 


C6 = (R2**2+R3**2+AL**2 ) **2- (2.*R2*R3)**2 
C 7 = (R2**2+R4**2+AL**2 ) **2- ( 2 . *R2*R4 ) **2 

C8 = (Rl**2+R4**2+AL**2)**2-(2.*R1*R4)**2 

C9 = (R1**2+R3**2+AL**2)**2-(2.*R1*R3)**2 

F24 = { .5/<R2**2-Rl**2) )*( SORT ( C 6 ) -SORT ( C 7 )+SQRT (C8 I-SQRT (C9 ) ) 

RETURN 

END 
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$IBFTC SHCLC DECK 


SUBROUTINE SHCLC 

THIS ROUTINE CALCULATES THE RADIAL TEMPERATURE DISTRIBUTION OF A 
SHIELD ASSUMING THAT THE TEMPERATURES OF THE ADJACENT SURFACES ARE 
KNOWN 

DIMENSION ESI 10 ) , ABSRI10) , REFS(IO) , TS(IO) , ARESI10) 

COMMON / T2 / A(41,4l> , Q<41) 

EQUIVALENCE ( A( 1),B ) 

DIMENSION B( 41» 4 l ) 

DIMENSION 0(20,20 , 01(20,20) , DUMMY(20) 

COMMON / CONDT/ C0ND(20,10) , RES(20,10) 

COMMON/ NEW/ FSSHTTI 6 ) , FSSTTHI 6 ) , FSNHTT(20,6> , FSNTTH(20,6) , 

1 FNNHTTI 20,20,6) , FNN T TH I 20, 20 , 6 ) , AREA(20,7) , 

2 FNSHTT (20,6) , FNSTTH(20,6) , R(20,7) , AREAT ( 7 ) , 

3 FNHTSR ( 20, 6 ) , FNTTSR(20,6) , FSRTNH(20,6) , 

4 FSRTNT ( 20,6 ) , FSRSRI 6 ) 

COMMON / BULL/ EMS(20,9,2) , AB(20,9,2) , REF(20,9,2> 

COMMON / EAS / AT , AH , QOA , T(20,10) , F13 , F14 , F23 , F24 
COMMON / KSLV / 10(20,10) , LOOK , LCOND , ICOND , NCALC , TMAX 

COMMON C(IOOO) 

EQUIVALENCE 


1 ( C C 1 ) 

, ABSR 

) 

V 

( C C LI) 

9 ARES 

) 

f 

( C ( 2 1 ) 

9 DIS 

) 


2 ( C ( 3 1 ) 

fES 

) 

F 

( C ( 4 1 ) 

, IVCOND) 

* 

(C (42) 

, ICONOR) 


3<C (43) 

, I TMAX 

) 

9 

( C ( 5 1 ) 

9 N AH 

) 

t 

( C C 71 ) 

9 NAT 

) 


4 ( C ( 9 1 ) 

,NA 

) 

t 

( C ( 29 1 ) 

, NEH 

) 

t 

(C (311 ) 

9 NET 

) 


5(C (331) 

, NE 

) 

t 

( C ( 5 1 1 ) 

9 NC QND 

) 

t 

(C ( 541 ) 

9 NC ASE 

) 


6(C (542) 

t NELIPS) 

f 

( C ( 543 ) 

t NR AOS 

) 

* 

( C ( 544 ) 

9 NRNGS 

) 


7 ( C ( 54 5 ) 

* NSHLDS ) 

F 

( C ( 546 ) 

tNCALCR) 

9 

( C ( 54 7 ) 

t N VUFH 

) 


8 ( C ( 5 4 3 ) 

* NVUFT 

) 

f 

(C( 549) 

f PHI 

) 

9 

( C ( 5 50 ) 

fPI 

) 


9 ( C ( 5 5 1 ) 

t RS 

) 

f 

( C ( 56 1 ) 

9 RS H 

) 

1 

( C ( 562 ) 

, RST 

) 


A ( C ( 5 6 3 ) 

,RH 

) 

t 

( C ( 564 ) 

t RT 

) 

t 

(C( 565) 

9 SOL AR 

) 


8 ( C ( 566 ) 

9 S I GM A 

) 

9 

(C( 567) 

9 TOL 

) 

f 

(C ( 568) 

fWI 

) 


C ( C ( 5 7 1 ) 

9 SQL A BS ) 

f 

( C ( 59 1 ) 

9 STHICK ) 

F 

(C (601 ) 

9 C VA8 S 

) 


D(C (641) 

y C VCOND ) 

f 

( C ( 681 ) 

9 CVEM S 

) 

t 

(C ( 721 ) 

f TH 

) 


E(C( 741) 

, TS 

) 

9 

(C( 751 ) 

9 TT 

) 

F 

( C ( 77 1 ) 

fREFS 

) 



DO 1 M= I »NSHLDS 
SUMQT=0.0 
SUMQX=0 . 0 

IF ( M .NE. 1 .OR. SOLAR .EQ. 0.0 ) GO TO 15 

CALCULATE HEAT TRANSFER DUE TO SOLAR FLUX. 

DO 17 N= 1 » NRNGS 

IF ( N .EQ. 1 ) G( 1)=PHI*PI*R( 1,2)**2*S0LABS/SIGMA 

IF ( N .NE. 1 ) Q(N)=PHl*PI*(R(N,2)**2-R(N-l»2)**2)*S0LABS/SIGMA 

IF ( ICOND .EQ. 2 ) GO TO 7 

DO 18 NN=1, NRNGS 

0 ( N , NN ) =0 . 0 

IF ( N. EQ. NN ) DIN, NN ) =-EM S ( NN ,2 , 1 ) 
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18 CONTINUE 
GO TO 17 

7 SUMQX=SUMQX-EMS( NN»2« 1)*AREA(NN,2 > 

17 CONTINUE 
GO TO 16 
C 

C CALCULATE HEAT TRANSFER FROM SURFACE CLOSER TO HEATER. 

C 

15 CALL COEF(M,B) 

IF ( ICOND .EQ. 2 ) GO TO 11 

DO 30 N= 1 , NRNGS 

SUMQ=0. 0 

L=N+NRNGS 

DO 21 NN=1, NRNGS 

SUMQ- SUMQ + 8(NN,L ) *T ( NN , M ) **4*EMS ( NN , M , 2 > *ARE A ( NN , M ) 

K=NN+NRNGS 

DK=0 • 0 

IF ( N .EQ. NN ) DK = 1.0 
D ( N , NN ) = ( B(K,L >-DK ) *EMS ( NN, M+l , 1 ) 

21 CONTINUE 

Q(N)=SUMQ + B(2*NRNGS+1,L)*TS<M)**4*ES(M)*ARES<M) 

30 CONTINUE 
GO TO 16 

11 SUMQ T=0 • 0 
SUMQX=0 . 0 

DO 12 N= 1 » NRNGS 
L=N+NRNGS 
SUMQ=0. 0 
SUM 1 = 0. 0 

DO 13 NN=1, NRNGS 
K=NN+NRNGS 

SUM 1 = SUM 1 + B(K,L)*EMS(NN,M+1,1)*AREA(NN,M+1) 

SUMQ = SUMQ + B<NN,L)*EMS(NN,M f 2)*AREA<NN,M)*T{NN,M)**4 
13 CONTINUE 

SUMQX=SUMQX+SUM1 

Q < N ) = SUMQ + B( 2*NRNGS+1?L )*TS(M)**4*ES(M)*ARES (M) 
SUMQT=SUMQT-Q(N> 

12 CONTINUE 
C 

C CALCULATE HEAT TRANSFER FROM SURFACES CLOSER TO TANK. 

C 

16 CALL COEF ( M+l , A ) 

IF ( ICOND .EQ. 2 ) GO TO 2 

IF ( NCALC .EQ. 0 .OR. ICOND .NE. 0 ) GO TO 3 

DO 4 N= 1 f NRNGS 

SUMQ=0 . 0 

DO 5 NN = 1 y NRNGS 

K=NN+NRNG$ 

SUMQ = SUMQ + A(K,N)*EMS(NN,M+2, 1 ) *T { NN , M+2 > **4* ARE A < NN , M+2 ) 
DK=0 . 0 

IF { N .EQ. NN ) DK= 1.0 

D ( N» NN ) = ( ( A ( NN » N ) - DK ) *EM S ( NN , M+ 1 , 2 J + 0(N,NN) )* 

1 ARE A (NN, M+l )*T(NN, M+l )**3 

IF ( N .NE. NN ) GO TO 6 

IF ( N.NE. 1 ) D ( N , NN ) =D ( N , NN ) - COND ( N- 1 , M ) /RE S (N- 1 , M ) 

IF ( N .NE. NRNGS ) D ( N , NN ) = D ( N , N N ) - COND < N , M ) /RES ( N, M ) 
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GO TO 5 

6 IF ( N.EQ. NN+1 ) D(N,NN> = D<N,NN) + COND ( N- 1 , M ) /RES ( N- l , M ) 

IF ( N .EQ. NN-1 ) D(N,NN)=D(N,NN) + COND < N , M ) /RES ( N ,M > 

5 CONTINUE 

Q(N) = - Q ( N ) - SUMQ - A ( 2*NRNGS+1 »N)*TS(M+1)**4*ESCM+1)*ARES(M+1) 
4 CONTINUE 
GO TO 9 

3 DU 8 N=l, NRNGS 
SUMG=0.0 

DO 10 NN=1, NRNGS 
K=NN+NRNGS 

SUMQ = SUMQ + A(K,N)*T(NN,M+2)**4*EMS(NN,M+2,l)*AREAINN,M+2) 

DK=0.0 

IF ( N .EQ. NN ) DK= l . 0 

D ( N , NN ) = ( A ( NN « N ) - DK ) *EMS ( NN , M+ 1 , 2 ) + D(N,NN) 

IF ( ICOND .NE. C ) GO TO 10 
IF ( N .NE. NN ) GO TO 28 

IF { N .NE. 1 ) C(N,NN)=D(N,NN)-CONO(N-L, MI/TIN, M+L)**3/RES<N-1,M> 
I F ( N .NE. NRNGS) 0 ( N , NN ) =D ( N, NN )-COND ( N f M ) /T ( N , M+ 1 ) **3/RES ( N , M ) 

GO TO 10 

28 I F ( N .EQ. NN+l)D(N,NN)=0(N,NN)*C0ND(N-l,M)/T(NN,M+l)**3/RES(N-l,M) 
I F ( N .EQ. NN-1) D ( N, NN ) =D( N, NN ) +COND (N,M)/T(NN»M+1)**3/RES(N,M) 

10 CONTINUE 

Q(N) = -Q(N) - SUMQ - A(2*NRNGSd,N)*TS(M+l)**4*ES(M+l)*ARES(M+l) 

8 CONTINUE 

9 IF ( NRNGS -NE. 1 ) GO TO 14 

Did. 1) = 1. 0/0(1, 1) 

GO TO 19 

14 CALL FACTOR ( D , DUMMY , NRNGS , 20 ) 

CALL I N V E R T ( D , DUMMY , NRNGS , 20 , DI ) 

19 DO 22 N=l, NRNGS 
SUM = 0 .0 

DO 23 NN=1, NRNGS 

SUM = SUM + D I ( N , NN ) *Q( NN ) 

23 CONTINUE 

IF ( NCALC .EQ. C .OR. ICOND .NE. 0 ) GO TO 24 
IF ( ABS ( SUM ) .GT. L.0E04 ) TMAX=1.0 
IF ( TMAX .NE. 0.0 ) RETURN 
IF ( SUM .LT. 0.0 ) SUM=0.8*T(N,M+1) 

T ( N » M + 1 )=SUM 
GO TO 22 

24 IF ( SUM .GT. L.0E15 .OR. SUM .LT. 0.0 ) TMAX=1.0 
IF ( TMAX .NE. 0.0 ) RETURN 

TIN.M+l) = { SUM/AREA1N.M+1) )**0.25 
22 CONTINUE 
GO TC 1 

CALCULATE SHIELD TFMPER ATURE FOR INFINITELY CONDUCTING SHIELD. 

2 DO 25 N=l, NRNGS 
SUM2=0. 0 
SUMQ=0. 0 

DO 26 NN* 1 , NRNGS 
K= NN +NRNGS 

SUM2 * SUM2 + A( NN.N )*EMS( NN .M+1,2 )*AREA (NN.M+l ) 

SUMQ= SUMQ + A{K,N)*nNN,M+2)**4*EMS(NN,Mt2,l)*AREA(NN,M4-2) 
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26 CONTINUE 

Q(N) = SUMQ + A( 2*NRNGS-H,N)*TS ( M+l )**4*ES (M+l ) *ARES(M+1 ) 

SUMGX = SUMQX + SUM2 - ( EMS(N f M+l,I) + EMS(N,M+lt2) ) *ARE A ( N f M + L ) 
SUMQT = SUMQT -Q ( N ) 

25 CONT IN J E 

T ( I » M + 1 ) = ( SUMQT /SUMQX) **0.25 
DO 27 N= I » NRNGS 
T ( N» M + l ) = T ( 1,M+I ) 

27 CONTINUE 
I CONTINUE 

RETURN 

END 
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SUBROUTINE COEF ( M M , B ) 

THIS SUBROUTINE DETERMINES THE VALUE OF THE COEFFICIENTS ( B ) 
USED TO FIND THE HEAT TRANSFER. 

DIMENSION ESI 10 } , ABSR(IO) , REFS(IO) , TS(LO) , ARES(IO) 

COMMON/NEW/ FSSHTT ( 6 ) , FSSTTHI 6 ) , FSNHTT(20,6) , FSNTTH(20,6) , 

1 FNNHTT(20,20,6) , FNNTTH ( 20 , 20 , 6 ) , AREA(20,7) , 

2 FNSHTT(20,6) , FNSTTH(20,6) , R(20 f 7) , AREAT ( 7 ) , 

3 FNHT SR ( 20 » 6 ) , FNTTSR(20,6) , FSRTNH(20,6) , 

4 FSRTNT ( 20* 6 ) , F SR SR ( 6 ) 

CGMMCN / BULL/ EMS(20,9,2> , AB(20,9,2) , REF(20t9,2) 

DIMENSION A ( 4 1» 4 1 ) , Y(41,4i> , A I I 4 1 , 4 1 ) , DUMMY ( 4 L ) 

DIMENSION B ( 4 If 4 1 ) 

COMMON C( 1000 ) 

EQUIVALENCE 


1 ( C ( 1 > 

, ABSR 

) 

t 

(C( 11 ) f AR ES 

) 

9 

(C (21) 

tO 

) 


2 ( C C 3 1 ) 

f ES 

) 

» 

(CI41) fIVCOND) 

9 

(C ( A2) 

, ICONDR) 


3(C (A3) 

t ITMAX 

) 

? 

( C ( 51 ) f NAH 

) 

9 

( C ( 7 1 ) 

* NAT 

) 


A { C ( 9 1 ) 

t NA 

) 

9 

(C ( 291 ) , NE H 

) 

9 

( C ( 3 1 1 ) 

? NET 

) 


t>(C(331) 

, NE 

) 

9 

( C ( 5 3 1 ) , NCOND 

) 

9 

(C ( 5A1 ) 

,NCASE 

) 


6 ( C ( 15 A 2 ) 

,NELIPS) 

9 

( C ( 543 > f NR ADS 

) 

9 

( C ( 5 AA ) 

» NRNGS 

) 


7 ( C ( 5 A 5 ) 

,NSHLDS) 

9 

(C( 5A6), .\CALCR) 

9 

( C ( 5A7 ) 

f N VUFH 

) 


8(C(5Atf ) 

t N VUF T 

) 

9 

( C ( 5A 9 ) » PH I 

) 

9 

( C ( 550 ) 

,PI 

) 


9{C (551 ) 

t RS 

) 

9 

(C( 56 1 ) t R S H 

) 

9 

( C ( 5 62 ) 

? RST 

) 


A ( C ( 56 3 ) 

, RH 

) 

9 

( C ( 56 A ) » RT 

) 

9 

( C ( 56 5 ) 

, SOLAR 

) 


B ( C ( 5 66 ) 

, SIGMA 

) 

9 

( C ( 56 7 ) t TOL 

) 

9 

( C ( 568 ) 

,WI 

) 


C ( C ( 5 7 1 ) 

,SOLABS) 

9 

( C ( 59 1 ) » ST H I C K ) 

9 

(C (601 ) 

,CVABS 

) 


D ( C ( 6 A 1 ) 

,CVCOND) 

9 

( C ( 68 1 ) , CVEMS 

) 

9 

(C ( 721 ) 

,TH 

) 


E ( C ( 7 A 1 ) 

*TS 

) 

9 

(C( 75 1 ) f TT 

) 

9 

(C ( 771 ) 

? REF S 

) 



M=MM 
ND I M = 4 1 
NS=NRNGS+ I 
NF = 2 *NRNGS 
DO 4 N= 1 f NRNGS 

CALCULATE THE VALUES OF THE MATRICES A AND Y BETWEEN ONE SURFACE 
AND ITSELF 

DO I NN= I f NRNGS 
Y( Nf NN) =0.0 

IF ( N .EQ. NN ) GO TO 2 
A ( N , NN ) =0 . 0 
GO TO 1 

2 A(N,NN)=— 1.0 
1 CONTINUE 

CALCULATE THE VALUES OF THE MATRICES A AND Y BETWEEN ONE SURFACE 
AND THE OTHER SURFACE 
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DO 3 NN=NS » NF 
N1=NN-NRNGS 

A(N»NNI*REF(Nl,M+l, l ) *FNNHTT ( N, N 1 , M ) 

Y { N i NN ) = — AB ( Nl» N + 1 » 1 ) *FNNHT T ( N ,N 1 ,M ) 

3 CONTINUE 

4 CONTINUE 

DO 5 N=NS»NF 
C 

C CALCULATE THE VALUES OF THE MATRICES A AND Y BETWEEN THE OTHER 

C SURFACE AND THE ORIGINAL SURFACE 

C 

NI=N-NRNGS 
DO 6 NN* 1 • NRNGS 

A(N,NN)=REF(NN,M,2)*FNNTTH(N1,NN,M) 

Y(N,NN)= -ABINN, R,2 > *FNNTTH { N1 , NN , M ) 

6 CONTINUE 
C 

C CALCULATE THE VALUES OF THE MATRICES A AND Y BETWEEN THE OTHER 

C SURFACE AND ITSELF 

C 

DO 7 NN = NS » NF 
Y(N»NN)=0.0 

I F ( N .EQ. NN ) GO TO 8 
A ( N » NN ) =0 . 0 
GO TO 7 

8 A I N» NN) =-I .0 

7 CONTINUE 

5 CONTINUE 
C 

C CALCULATE THE VALUES OF THE MATRICES A AND Y BETWEEN THE SURFFACES 

C AND THE SURROUNDINGS. 

C 

DO 9 N= 1 » NF 

IF ( N .GT. NRNGS) GO TO 10 
A ( NF + 1 1 N ) = REF(N»M»2)*FSRTNH(N»M) 

Y ( NF+ L, N ) = -AB(N,M,2)*FSRTNH(N,M ) 

GO TO 9 

10 N1=N-NRNGS 

A(NF+l,N)=REF(Nl,M+l, 1 ) *FSRTNT I N 1 , M ) 

Y(NF + 1,N)= -AB(N1,M+1, l ) *F SR TNT I N l , M ) 

9 CONTINUE 

DO 11 N=1 t NF 

IF ( N .GT. NRNGS ) GO TO 12 
A(N,NF+1 ) = REFS(M ) *FNHTSR ( N »M ) 

Y(N,NF+1) = - ABSR(M)*FNHTSR(N,M) 

GO TO 11 
12 NUN— NRNGS 

A(N»NF+l!=REFS(M)*FNTTSR(NltM) 

Y ( N * NF+ 1 ) = - ABSR(M)*FNTTSR(N1»M) 

11 CONTINUE 

A!NF+1,NF+1)= REFS (M )*FSRSR( M) - 1.0 
Y ( NF+ 1 » NF+ 1 ) = — ABSR(M)*FSRSR(M) 

NFZ=NF+1 

C 

C DETERMINE THE INVERSE OF THE MATRIX A 

C 
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CALL FACTOR ( A , DUMMY f NF + 1 , NDIM ) 

CALL INVERT ( A , DUMMY , NF+1 , NDIM , AI ) 

COMPUTE THE VALUES OF THE ABSORPTION COEFFICIENTS (B) 

DO 13 I=1,NFZ 

DO 13 J ~ 1 , N F Z 

SUM=0 • D 

DO 14 K=1,NFZ 

SUM= SUM + AI CI«K)*Y(K V J) 

14 CONTINUE 
B ( I f J ) = SUM 
13 CONTINUE 
RETURN 
END 


(B=AI*Y) 
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noon 


$ I BFTC QDOT 


DECK 


SUBROUTINE QDOT 

THIS SUBROUTINE DETERMINES THE HEAT TRANSFER TO M , M+ 1 , AND 
SURROUNDING SURFACES 

DIMENSION ES { 10 ) , ABSR(IO) , REFS(IO) , TS(10) , ARES(IO) 

COMMON / EAS / AT , AH , QOA , TC2Q,10) , F13 , F14 , F23 , F24 
COMMON/ NEW/ FSSHTTC6) , FSSTTHI6) , FSNHTT(20»6) , FSNTTH(20,6) , 

1 FNNH TT (20,20,6) , FNNTTH ( 20, 20 , 6 ) , AREA(20,7) f 

2 FNSHTT (20,6 ) , FNSTTH(20,6) , R(20,7) , AREAT ( 7 ) , 

3 FNHTSR (20,6) , FNTTSR(20,6) , FSRTNH (20*6) , 

4 FSRT N T ( 20 » 6 ) , FSRSR(6) 

COMMON / BULL/ EMS(20,9,2) , AB(20,9,2) , REF(20,9,2) 

COMMON /T3/ Q0A2 , QSUR , QOAS 
COMMON / T2 / 8(41,41) , Q(41) 

COMMON C( 1000) 

EQUIVALENCE 

1 ( C I I ) ,ABSR ) , C C C 1 1 ) ,ARES ) , (C(21) ,0 ) , 

2 ( C ( 3 1 ) , ES ) , ( C C 41 ) ,IVC0ND) , ( C C 42 ) ,ICONDR) , 

3 ( C ( 43 ) ,ITMAX ) , (C(51) , NAH ) , (C(71) , NAT ) f 

4 ( C ( 9 1 ) , NA ) , ( C ( 29 1 ) , NEH ) , (C(311),NET ) , 

5 ( C ( 3 3 1 ) , NE ) , ( C ( 53 1 ) , NCOND ) , ( C ( 541 ) , NC ASE ) , 

6(C (542) ,NELIPS) , ( C ( 543 ) , NR ADS ) , ( C ( 544 ) , NRNG S ) , 

7 (C (545) , NSHLDS ) , ( C ( 546 ) , NC ALCR ) , ( C ( 547 ) , NVUFH ) , 

7 ( C ( 545 ) , NSHL DS ) , ( C ( 546 ) , NS TE P ) , ( C ( 54 7 ) , NVUFH ) , 

8 ( C ( 548) ,NVUFT ) , (C(549),PHI ) , (C(550),PI ) , 

9 ( C ( 55 1 ) , RS ) , ( C ( 56 1 ) , RS H ) , (C(562 ),RST ) , 

A(C(563) ,RH ) , ( C ( 56 4 ) , RT ) , ( C ( 565 ) , SOL AR ) , 

8(C(566) , SIGMA ) , (C(567),T0L ) , (C(568),WI ) , 

C(C ( 571) tSOLABS) , ( C ( 59 1 ) , STHI CK ) , ( C ( 601 ) , C VABS ) , 

D(C (641) tCVCOND) , ( C ( 68 1 ) , CVEMS ) , (C(72i),TH ) , 

E ( C ( 7 4 1 ) *TS ) , ( C ( 751 ) , TT ) , (C(77i),REFS ) 

M=NSHLDS+ 1 

IF ( M .NE. 1 .OR. SOLAR .EQ. 0.0 ) GO TO 9 
C 

C CALCULATE THE HEAT TRANSFER DUE TO SOLAR FLUX. 

C 

QSUM=0. 0 

DO 8 N= 1 , NRNGS 

IF ( N .EQ. 1 ) Q(N)=PHI*PI*R{ 1 , 2 ) **2* SGLAB S 

IF ( N .NE. 1 ) C(N)=PHI*PI* (R(N,2 )**2 - R ( N— 1 , 2 ) **2 ) *S0L ABS 
Q ( N ) - Q ( N ) /SIGMA - EMS ( N , 2 , 1 ) *ARE A ( N , 2 ) *T ( N , 2 ) **4 
QSUM=QSUM + Q (N ) 

8 CONTINUE 

Q0A=QSUM/AREAT(2 )*SIGMA 

QO AS=0. 0 

QSUR=0.0 

Q0A2=0 . 0 

RETURN 

9 CALL COE F ( M , B ) 
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ooo ooo o o o o 


QDOT 


EFN SOURCE STATEMENT 


IFN(S) 


NS=NRNGS + 1 

NF=2*NRNGS 

NZ=NF+1 

CALCULATE THE HEAT TRANSFER TO THE TANK. Q ( J+NRNGS ) IS THE HEAT 
TRANSFERED TO THE J TH ELEMENT ON THE TANK 

QSUM=0.0 
DO 13 N=NS,NF 
SBE=0.0 
DO 14 NN= 1 » NF 

I F ( NN .GT. NRNGS ) GO TO 15 

SBE= SBE + B(NN,N>*T(NN,M)**4*EMS(NN,M,2)*AREA(NN,M) 

GO TO 14 

15 N 1 =NN— NRNGS 

SBE = SBE + B(N:N,N)*T(N1 ,M + 1 )**4*EMS ( Nl , M+ l , l ) * ARE A ( N1 , M+ 1 > 

14 CONTINUE 

SBE = SBE + B(NF+ltN)*TS(M) **4*ES (M)*ARES(M) 

N1=N-NRNGS 

Q ( N ) = SBE - T(N1»M+1 )**4*EMS(Nl,M+l, 1 )*AREA(N1 f M+l ) 
QSUM=QSUM+Q(N) 

13 CONTINUE 

QOA=QSUM/AREAT(M+l ) *S I GMA 

CALCULATE THE HEAT TRANSFER TO THE SURFACE ADJACENT TO THE TANK. 

QSUM=0 . 0 

NF=2*NRNGS 

DO 10 ,N=1, NRNGS 

SBE=0.0 

DO 11 NN= l , NF 

IF ( NN.GT. NRNGS ) GO TO 12 

SBE = SBE + B(NN,N)*T(NN,M)**4*EMS(NN,M,2)*AREA(NN,M) 

GO TO 11 
12 N1=NN— NRNGS 

SBE= SBE + B ( NN» N ) * T ( N 1 » M+ 1 ) **4#CMS ( N l » M+ 1» l ) * ARE A ( N 1 » M+ 1 ) 

11 CONTINUE 

SBE=SBE + B( NF+1 ,N ) *TS( M )**4*ES( M ) *ARES( M) 

Q(N)= SBE -T(N»M )**4*EMS(N,M,2)*AREA(N,M) 

QSUM=QSUM + Q ( N ) 

10 CONTINUE 

CALCULATE THE HEAT TRANSFER TO THE SURROUNDINGS. 

ADVD= ARE AT (M ) 

IF ( M .EQ. 1 .AND. SOLAR .NE. 0.0 ) ADVD=AREAT (M+l ) 

Q0A2=QSUM/ ADVD*S IGMA 
QSUR=0. 0 
DO 17 N=lt NF 

IF ( N .GT. NRNGS ) GO TO 16 

QSUR = QSUR + B(N,NF+l)*T(NtM)**4*EMS(NtM,2)*AREA(N,M) 

GO TO 17 

16 Nl=N-NRNGS 

QSUR=QSUR+B{ N,NF+ 1 ) *T (N1 ,M+1 ) 4<=EMS ( N l » M+ 1 , l ) * ARE A { N l , M+ 1 ) 

17 CONTINUE 

QSUR = QSUR + ( B(NF+l,NF+l) - 1.0 ) *TS ( M ) **4*E S ( M ) * ARE S ( M ) 
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UU I 


tt- N 


SUUKCt STATEMENT 


I FH( S ) 


Wi>UK =w iUrC* Si o*‘1A 
wu AS = w SUK / aklS( i'1 ) 
Kt TU kn 
ENU 
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ooooooooono 


$1 BFTC FACTOR DECK 


SUBROUTINE F ACTOR ( A, ROW , ORDE R, D IM ) 

FACTOR THE MATRIX A INTO TWO TRIANGULAR MATRICES ONE OF WHICH 
HAS A UNITARY DIAGONAL. 

A THE MATRIX OF SIZE ORDER WHICH IS TO BE INVERTED. 

D THE INVERSE OF A. ( USED IN INVERT ) 

ORDER THE LENGTH OF EACH SIDE OF THE MATRICES A AND D. 

DIM THE DIMENSION OF THE MATRICES A AND D AND THE VECTOR ROW 

IN THE CALLING PROGRAM. 

ROW A DUMMY VARIABLE OF LENGTH DIM TRANSFERED FROM THE CALLING 

PROGRAM. 

INTEGER ONE, TWO, ROW, ORDER, DIM 
DIMENSION A( DIM, DIM I ,R3W(DIM ) 

DOUBLE PRECISION T,ZERO 
REAL MAX 

DATA ZERO/O. ODD/ 

N=ORDER 

ASSIGN 150 TO ONE 
ASSIGN 160 TO TWC 
DO 95 1 = 1, N 
I P 1= I + 1 
I M 1= I- l 

GO TO ONE, { 10,150 
10 DO 30 L= I , N 
T=ZERQ 

DO 20 K= 1 , I M 1 
20 T= T + A ( K » l ) *A ( L» K ) 

30 A(L,I)=A(L«I )-T 
40 MAX=0 .0 

DO 50 L= I » N 

IF (ABS(A(L,I ) l.LE.ABS(MAX) ) GO TO 50 
N I = L 

MAX= A ( L , I ) 

50 CONTINUE 

IF (MAX. EQ. 0.0 )G0 TO 140 
ROW ( I ) = N I 

IF ( I .EQ.NI )G0 TO 180 
DO 60 K=l,N 
S= A ( I , K ) 

A ( I,K) = A(NI,K) 

60 A ( N I , K ) =S 
65 GO TO TWO, (70,160 
70 DO 90 L= I P 1 » N 
T=ZERO 

DO BO K=l» IM1 
80 T=T+ A ( I »K)*A(K»L ) 

90 A ( I , L ) = ( A ( I , L )-T ) /MAX 
95 CONTINUE 
140 ORCER= I M 1 
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FACTOR 


EFN SOURCE STATEMENT 


I FN ( S ) 


WR I T E ( 6 f 141 ) 

141 FORMAT ( / 53H SUBROUTINE FACTOR HAS ENCOUNTERED A SINGULAR MATRIX.) 
CALL EXIT 

150 ASSIGN 10 TO ONE 
GO TO 40 

160 DO 170 L = 2»N 
170 A(lfL)-A(lfL)/MAX 
ASSIGN 70 TO TWO 
GO TO 95 

180 I F ( I.EQ.NJRETURN 
GO TO 65 
END 

ERROR MESSAGE NUMBER 1 

ERROR MESSAGE NUMBER 2 

ERROR MESSAGE NUMBER 3 
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non 


$ I BFTC INVERT DECK 


SUBROUTINE I N VER T ( A , ROW , OR DER, D I M , D ) 

DETERMINE THE INVERSE OF A FROM THE TWO TRIANGULAR FACTORS OF A 

INTEGER ROW, ONE, CIM, ORDER 
DIMENSION At DIM, DIM) , RO W ( D I M ) , D ( D I M , D I M ) 

DOUBLE PRECISION T,ZERO 
DATA l ERO/O.ODO/ 

N=CRDER 

ASSIGN 60 TO ONE 
I = N 

5 IP1=I+1 
IMl=I-i 
DO 30 M=l, I 
L=IP1-M 
T= ZE RQ 

GO TO ONE, ( LO ,60 ) 

10 LP1=L+1 

DO 20 K = L P 1 , N 
20 T=T+A(K,L)*D< I,K) 

30 D( I ,L)= ( FLOAT (L/ I )-T )/A( L,L ) 

I F ( I . EQ • L ) GO TO 70 
DO 5 0 M = 1 , I M 1 
L= I -M 
T = Z E R 0 
L P 1 = L + 1 
DC 40 K = L P 1 , N 
40 T=T+A{L,K)*D(K, I ) 

50 0 ( L , I ) =— T 
I = I M 1 
GO TC 5 

60 ASSIGN 10 TO ONE 
GO TC 30 
70 I=N~1 
80 K= ROW ( I ) 

DO 90 L= 1 , N 
S= D ( L , I ) 

D ( L , I )=DtL,K) 

90 D ( L , K ) = S 

I F ( I • EQ . 1 ) R £ TURN 
1 = 1-1 
GO TO 80 
END 

ERROR MESSAGE NUMBER 1 



listing of input 


*OUTPUl 




I CONDR= 

1, ITMAX = 

2C0» 

I VCONO= 

NEL I PS= 

0, NRADS = 

l. 

NRNGS = 

NVUFT = 

0, PHI 

4 .4200000E 02, 

RH 

RT = 

1 . OOOOOOOF CO, SIGMA = 

1.7130000E-C9, 

SOLAB S= 

S END 





TBABSl 0.2000E-01 0. 0. 0. 

SHIELD RADIUS = 0.100CE 01 O.iOOOE 01 


1, NCALCR= 


1* 

NCASE 

= 

It 

20, NSHLOS= 


2, 

NVUFH 

- 

0, 

0 . , RS H = 

0. 

» 

RST 

= 

0. , 

A • 9 99 9 9 99 E- 02 , SOLAR = 

l.OOOOOOOE 

00, 

TOL 

= 

1 • 0000 00OE-02 , 

TB6MSL 0.2000E-01 0. 

0. 



0. 



SPACING 

2 

3 


DISTANCE 
0.100 
0. LOO 
NELB = 1 

N A 1 B = 1 

NE1T = l 
NA IT = l 
NE2R = 1 

NA2B = l 
NE2T = 1 

NA2T = l 
NET = 1 

NAT = 1 

TT = 


AREA 
L.000 
1 • 000 


PROPERTIES OF 
ABSORPTIVITY 
1.000 
1.000 


SURROUNDINGS BETWEEN SURFACES 
EMISSIVITY TEMPERATURE 
1.000 0 . 

l.COO 0. 


37.000 


LISTING OF OUTPUT 


NO. ITERATIONS = 4 


SOLAR 


QOAT * 0.L188E 00 


TOLX=0.OlCOO 


M 

l 


M 

2 


TEMP 

RtN) 

RtN) /R S 

772.855 

0.224 

0.224 

772.292 

0.316 

0.316 

771.660 

0.387 

0.307 

771.002 

0.447 

0.447 

770.325 

0.500 

0.500 

769.630 

0.540 

0.548 

768.918 

0.592 

0.592 

768.188 

0.632 

0.632 

767. A40 

0.671 

0.671 

766.672 

0. 707 

0. 707 

765.882 

0. 742 

0.742 

765.069 

0.775 

0.775 

764.231 

0.006 

0.806 

763.363 

0.837 

0.837 

762.462 

0.866 

0.866 

761.525 

0.894 

0.894 

760.548 

0.922 

0.922 

759.530 

0.949 

0.949 

758.486 

0.9 75 

0.975 

757.460 

1.000 

l.OOC 


TEMP 

RtN) 

RtN) /RS 

443.449 

0.224 

0 . 22 A 

440.095 

0.316 

0. 316 

436.257 

0. 387 

0.387 

432.164 

0.447 

0.AA7 

427.846 

0.500 

0.500 

423.298 

0.540 

0.5A8 

418.505 

0.592 

0.592 

413.443 

0.632 

0.632 

408.081 

0.671 

0.671 

402.378 

0.707 

0.707 

396.284 

0. 742 

0.7A2 

389.737 

0. 775 

0,775 

382.650 

0.006 

0,806 

374.913 

0.837 

0.837 

366.367 

0.866 

0.866 

356.787 

0.894 

0.894 

345.835 

0.922 

0,922 

333.038 

0.949 

0.949 

317.852 

0.975 

0.975 

300.015 

1.000 

1,000 
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0.1188E 00 


T0LX=0. 00500 


NO. ITERATIONS * 2 


SOLAR 


COAT * 


M 

i 


M 

2 


TEMP 

R ( N) 

R<N)/RS 

772.855 

0.224 

0.224 

772.292 

0.316 

0.316 

771.660 

0.387 

0.387 

771.002 

0.447 

0.447 

770.325 

0.500 

0.500 

769.630 

0.548 

0.548 

768.918 

0.592 

0.592 

768.108 

0.632 

0.632 

767.440 

0.671 

0.671 

766.672 

0.707 

0.707 

765.882 

0.742 

0.742 

765.069 

0.775 

0.775 

764.231 

0.806 

0.806 

763.363 

0 . 837 

0.837 

762 .462 

0.866 

0-866 

761.525 

0.894 

0.894 

760.548 

0.922 

0.922 

759.530 

0.949 

0.949 

758.486 

0.975 

0.975 

757.460 

1.000 

1.000 


TEMP 

RIN) 

R(N)/RS 

443.449 

0.224 

0.224 

440.095 

0.316 

0.316 

436.257 

0.387 

0.387 

432.164 

0.447 

0-447 

427.846 

0.500 

0.500 

423.298 

0.548 

0.548 

418.505 

0.592 

0-592 

413.443 

0.632 

0.632 

408.081 

0.671 

0.671 

402.370 

0.707 

0.707 

396.284 

0.742 

0.742 

389.737 

0.775 

0.775 

382.650 

0.806 

0.806 

374.913 

0.837 

0.837 

366.367 

0.866 

0.866 

356.787 

0.894 

0.894 

345.835 

0.922 

0.922 

333.038 

0.949 

0.949 

3L7.852 

0.975 

0.975 

300.015 

1.000 

1.000 
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Solar radiation 



Figure 1. - Schematic drawing of two possible shadow shield configurations. 



Figure 2. - Representation of hollow 
structural member. Inside surface 
diameter of strut, d t ; outside sur- 
face diameter of strut, d Q ; strut 
length, l; temperatures at ends of 
strut, T} and T 2 ; distance along 
strut, x. 



Heat-transfer rate, Qi/A 1( W/m 



o 


Targeted 

e ~ q 94 ^ section 



.2 .4 .6 .8 1.0 

Radius ratio, r/R 


Figure 5. - Comparison of analytic and experimental 
temperature distributions for targeted shield 
system. Spacing ratio, Lf/R = 0.47; two shields. 


Temperature ratio, (T -T^/O^ - T^) 



Dimensionless axial distance, xll 


Figure 6 . - Comparison of analytic and experimental strut 
temperature distributions. 
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Emissivity, € 

Figure 11. - Heat -transfer rate as function of emissivity for 
different surface assumptions. No shields. 


Heat -transfer rate, Qf/oA^ 


Surface Radiosity 


Diffuse Nonuniform, nonconducting shield 

Specular Uniform 
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Figure 12. - Heat-transfer rate as function of emissivity for different 
surface assumptions. One shield. 





Heat-transfer rate, Q^/aA^ 



10' 2 10' 1 10° 
Emissivity, e 


Figure 14. - Heat-transfer rate as function of emissivity and 
spacing ratio for two oblate spheroids. No shields. 


Heat -transfer rate, QjjoA^ 


ffl 



Figure 15. - Heat-transfer rate as function of emissivity and 
spacing ratio for two oblate spheroids. One shield. 





Heat-transfer rate, -Qi/oAjTj 





Heat-transfer rate 






(a) Uniform temperature shields. (b) Nonconducting shields. 

Figure 20. - Effect of targeting on heat-transfer rate for uniform temperature and 
nonconducting shields. 


T^0 | ^ 


t 2 *o| ^ 


1 2>° C^ 
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' rU Lz^~~J |Targeted 

| area 
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Targeted 


_ 1 
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0 .2 .4 .6 .8 1.0 


0 .2 .4 .6 .8 1.0 

Radius ratio, r/R 


0 .2 .4 .6 .8 1.0 


(a) No shields targeted; heat- 
transfer rate, Qi/oAiTo, 

0. 64xl0~ 4 . 


(b) One shield targeted- heat- 
transfer rate, Qf/oAiT?, 

0. 70x10 . 


(c) Two shields targeted; heat- 
transfer rate, Qf/aAiTo, 

0. 78x10" . 


Figure 21. - Effect of targeting on radial temperature distributions for nonconducting shields. Spacing 
ratio, 0.3; emissivity of nontargeted surfaces, 0.03. 
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Ratio of heat-transfer rates to cold source, Q/Q k= 
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Ratio of cold source heat transfer to 
that for evenly spaced shields 


O 




(a) One shield. 



0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 


Dimensionless spacing, L/Lj 

(b) Two shields. (c) Three shields. (d) Two shields; variable 

spacing. 

Figure 23. - Effect of variation in shield position on heat-transfer rate. Spacing ratio, 0.2. 



(e) Three shields; variable 
spacing. 




Spacing ratio, Ly/R 

Figure 25. - Heat -transfer rate to source in absence of solar flux as 
function of spacing ratio and emissivity. Two shields. 




Total heat-transfer rate, 



(a) Ratio of strut length to 
diameter, 5 ; radiation to 
conduction parameter, 
125. 



(b) Ratio of strut length to 
diameter, 10, radiation 
to conduction parameter, 
500. 


External 

emissivity, 



(c) Ratio of strut length to 
diameter, 20, radiation 
to conduction parameter, 
2000. 


Figure 26. - Total heat-transfer rate through strut as function of internal and external emissivity. 
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Ratio of conducted heat transfer, 





.2 .4 .6 .8 1.0 


(a) Ratio of strut length to 
diameter, 5; radiation to 
conduction parameter, 
125. 


0 .2 .4 .6 .8 1.0 

Dimensionless axial distance, x/l 

(b) Ratio of strut length to 
diameter, 10; radiation 
to conduction parameter, 
500. 


0 .2 .4 .6 .8 1.0 


(c) Ratio of strut length to 
diameter, 20; radiation 
to conduction parameter, 
2000 . 


Figure 27. - Effect of external emissivity and length to diameter ratio on strut temperature distribution. End 
temperature ratio, 10; internal emissivity, 0.6. 


End temperature 
ratio, 

T 2 /Ti 


I I I i I 1 1 1 1 t I I I I 1 1 1 1 I l LHjJjj 1 I I i I 1 1 jj 

.2 . 4 . 6 1 2 4 6 10 20 40 60 100 200 400 600 1000 

Radiation to conduction parameter, k , 0 € 0 T 2 l 2 /kt 

Figure 28. - Effect of radiation to conduction parameter on conducted heat-transfer rate for constant end 
temperature ratios. Internal emissivity, 0. 


Temperature ratio, (T - Ti)/(To - Ti) 



Dimensionless axial distance, x/i 


(a) End temperature (b) End temperature (c) End temperature 

ratio, 3. ratio, 5. ratio, 15. 

Figure 29. - Temperature distributions for constant values of radiation to conduction 
parameter and end temperature ratios. Internal emissivity, 0. 



Figure 30. - Schematic drawing of heat balance for element of shield. 
Surface being considered, j; conducted heat transfer rate, q c ; 
radiant heat transfer rate, q r 
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Figure 31. - Schematic drawing of element numbering used in heat 
balances. Number of elements on each surface of each system, 
n, 5; total number of elements in system, k, 2n + 1. 
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Figure 32. - Apparent spacing between surfaces due to specular refiectivity. 
Surfaces 1 and 3 are specularly reflecting and 2 and 4 are windows. 
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Figure 33. - Schematic drawing of 
reflected energy for specular 
surfaces. 
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Figure 34. - Sources of heat transfer to element of strut. Inside diameter 
of strut, d j ; outside diameter of strut, d 0 ; conducted heat-transfer rate 
toward element, q in c ; radiant heat-transfer rate toward surface, qj n r ; 
heat-transfer rate from surroundings, qj n sr ; conducted heat-transfer 
rate away from element, q out c ; radiant he’at-transfer rate away from 
surface, q ou t r ; temperature’ at ends of strut, Tj and I 2 ; temperature 
of end disks, f T n+i and T n+2 ; distance along strut, x. 
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Figure 35. - Ratio of heat-transfer rate to that for single 
annulus as function of reciprocal of number of annuli. 
One shield. 







Figure 37. - Flow chart of routine SHIELD. 
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Figure 38. - Flow chart of subroutine TBVUF. 
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Figure 39. - Flow chart of subroutine VUFAC1. 
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Figure 40. - Flow chart of subroutine SUBVF. 
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Figure 41. - Flow chart of subroutine SHCLC. 
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Figure 42. - Flow chart of subroutine COEF. 
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Figure 43. - Flow chart of subroutine QDOT. 
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